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ABSTRACT 



In this thesis we treat the “small cell” problem encountered when building an attri- 
tion rate generator for large-scale manpower flow models, specifically for the USMC 
Officer Corps. Such models have a large number of low-inventory (i.e. small) personnel 
cells. This presents a dilemma: on one hand we want to preserve as much fidelity as 
possible in our work by preserving a great deal of detail in each cell; on the other hand 
our statistical estimation techniques require larger cell sample sizes than intrinsically 
occur cell-by-cell in actual sample data. Our approach to producing stable attrition rates 
for such cells involves two efforts: (i) the aggregation of cells into groups that exhibit 
homogeneity of attrition behavior, and (ii) the development of “shrinkage” estimation 
techniques for use in the individual groups. A heuristic algorithm is developed and 
tested to treat the aggregation problem. Empirical Bayes methods are developed to serve 
the multi-cell estimation requirements needed to preserve the fidelity. Cross validation 
techniques are used to verify these methods. 

The present work builds upon the results of previous studies; we integrate what was 
learned into a coherent package that is ready for use. 
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THESIS DISCLAIMER 

The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er- 
rors, they cannot be considered validated. Any application of these programs without 
additional verification is at the risk of the user. 



IV 



TABLE OF CONTENTS 

I. INTRODUCTION 1 

A. GENERAL 1 

B. BACKGROUND 2 

C. ORGANIZATION 3 

D. SMALL CELL PROBLEM 4 

E. DATA BASE 5 

II. CELL AGGREGATION 7 

A. GENERAL 7 

B. BACKGROUND 7 

C. EXPANSION METHOD 9 

D. AGGREGATION METHOD 14 

III. ESTIMATION METHODS 16 

A. GENERAL 16 

B. EMPIRICAL BAYES 17 

1. Transformed Scale 17 

2. Original Scale 20 

C. EFRON-MORRIS METHOD 23 

D. VECTOR METHOD 25 

IV. CROSS VALIDATION 28 

A. GENERAL 28 

B. MEASURES OF EFFECTIVENESS 28 

1. .Mean Absolute Deviation 28 

2. Chi Square 29 

3. .Mean Squared Error 30 

4. Vector Method .MOEs 31 

C. TEST CASES 31 

D. RESULTS 35 



V 



V. CONCLUSIONS AND RECOMMENDATIONS 45 

A. CONCLUSIONS 45 

B. RECOMMENDATIONS 46 

APPENDIX A. AGGREGATION ALGORITHMS 47 

A. HEURISTIC ALGORITHM 47 

B. INTEGER LINEAR PROGRAM 47 

APPENDIX B. COMPUTER PROGRAMS 49 

A. GENERAL 49 

B. MAIN PROGRAM AND AGGREGATION SUBROUTINES 51 

C. ESTIMATION SUBROUTINES 73 

D. VECTOR METHOD SUBROUTINE 88 

E. EXEC PROGRAM 94 

F. SAMPLE DATA FILE 95 

APPENDIX C. SAMPLE OUTPUT ? 98 

A. GENERAL 98 

B. SAMPLE OUTPUT (TEST CASES 1-30) 98 

C. SAMPLE OUTPUT (VECTOR TEST CASES) 101 

REFERENCES 102 



INITIAL DISTRIBUTION LIST 



104 



LIST OF TABLES 



Table 1. MOS GROUPS 11 

Table 2. YCS EXPANSION BOUNDS 12 

Table 3. TEST CASES FOR METHODS 1-5 33 

Table 4. TEST CASES FOR VECTOR METHOD 34 

Table 5. SUMMARY OF RESULTS (CASES 1-10) 40 

Table 6. SUMMARY OF RESULTS (CASES 11-20) 41 

Table 7. SUMMARY OF RESULTS (CASES 21-30) 42 

Table 8. SUMMARY OF RESULTS (CASES 7 AND 9 EXPANDED) 43 

Table 9. SUMMARY OF RESULTS (VECTOR METHOD) 44 



vii 



I. INTRODUCTION 



A. GENERAL 

The Officer Planning and Utility System (OPUS), a comprehensive and fully inte- 
grated manpower management system, is currently being implemented by the U.S. Ma- 
rine Corps (Decision System Associates, 1986). This system contains a set of 
computer-based manpower planning models and is used by the Officer Plans Section 
(MPP-30), Headquarters, U.S. Marine Corps, to produce several manpower planning 
documents. The system must be able to accurately predict personnel attrition, i.e., offi- 
cers leaving the service for purposes such as resignation, retirement, discharge, disability, 
or release. The forecasting of attrition is accomplished by the Marine Corps Officer 
Rate Projector (MCORP), developed by the Navy Personnel Research and Development 
Center (XPRDC), San Diego, California (XPRDC, 1985). 

The attrition rate generator developed by XPRDC calculates empirical attrition 
rates using historical data with user-defined weights and threshold parameters (Siegel, 
1983). This subjective input makes the current generator susceptible to unintentional 
misuse. 

In support of .MCORP, Professor Robert R. Read of the Xaval Postgraduate School 
has been working on the “small cell” problem: applying multiparameter statistical esti- 
mation schemes to estimating attrition when there is low personnel inventoiy, or small 
cells, which generally exhibit unstable empirical rates. 

A comment on terminology is in order. By attrition rate generator we mean meth- 
odology for estimating attrition probabilities for the various cells. The expression “em- 
pirical rates” refers to the ratio of leavers to inventory' for each cell, unmodified by any 
information contained in “neighboring” cells. In contrast to this, the expression “em- 
pirical Bayes” refers to Bayes estimators whose prior parameters are estimated from 
data. 

Accurate forecasting of losses is extremely important to the manpower planner. 
Overestimating losses causes excess accessions, promotion delays, underutilized person- 
nel and increased costs, while underestimation causes personnel shortages and decreased 
readiness. The problem is compounded in that all but a few accessions must start at the 
bottom, i.e.. Second Lieutenant, and work their way up to the higher ranks only after 
many years of service. For example, if a shortage of Lieutenant Colonels arises, it can 
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only be remedied by promoting more Majors, which has a rippling efTect down the rank 
structure. 

B. BACKGROUND 

There have been seven Master's theses over the past four years which have studied 
various aspects of the attrition estimation problem. A concise summary of these works 
is given by Read (NFS Report NPS55-88-006, 1988, pp. 16-23). These studies can be 
grouped into three general areas: shrinkage methods, cell aggregation and peripheral 
studies. 

The application of a shrinkage method begins by identifying a number of personnel 
inventory cells, followed by the development of the empirical rates for individual cells 
and a weighted grand mean of these empirical rates. The final estimate for a cell is a 
convex combination of its empirical rate and the grand mean. There are numerous 
methods for accomplishing this, several of which have been applied in previous studies. 

Tucker was the first to investigate the application of these methods to attrition es- 
timation. He compared traditional estimators to the James-Stein estimator and the 
minimax estimator for a few selected grades and occupational fields. His results gave 
strong support to the James-Stein estimator; minimax was discarded as being too con- 
servative for small cell use. However, there remained pockets of instability for which 
goodness-of-fit tests failed. (Tucker, 1985) 

Following Tucker was Robinson, who introduced the Efron-Morris limited trans- 
lation shrinkage alternative to augment the James-Stein estimator. These methods were 
evaluated with a broader set of test cases. Robinson was able to confirm Tucker's re- 
sults, but the limited translation option failed to provide any consistent relief in the un- 
stable areas. (Robinson, 1986) 

The final application of shrinkage methods to estimating officer attrition rates was 
undertaken by Dickinson. He applied the previously used methods and an empirical 
Bayes estimator to a new and refined data base. Improved results were obtained, but 
the instability remained. Dickinson also performed some exploratory side studies dealing 
with the Freman-Tukey transform and the use of empirical Bayes methods that allow 
non-uniform shrinkage, both of which provided the impetus for the present study. 
(Dickinson, 1988) 

These three studies used ad hoc methods to deal with the second general area of 
study-cell aggregation. Aggregation of cells with low personnel inventory into sets of 
cells, often of larger inventory, is required when applying these shrinkage methods. The 
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desire is to use cells which exhibit similar attrition behavior. Two previous studies have 
investigated this area. 

Amin Elseramegy used the Classification and Regression Trees (CART) program, 
which at the time was newly acquired by the Naval Postgraduate School, in an attempt 
to form aggregates of cells that exhibited homogeneous attrition behavior. Several dif- 
ficulties in using this program were encountered, e.g., because of insufficient memor}' 
allocation he found it necessary' to partition the data base into nine sets and apply 
CART to each. The resulting aggregations were generally unusable. (Amin Elseramegy, 
1985) 

Major breakthroughs in cell aggregation were made by Larsen. He applied a hi- 
erarchical clustering algorithm to the new data base. The resulting rules for building 
aggregates are well defined and especially viable from an intuitive point of view. 
Larsen's work provides the framework for the cell aggregation method developed in 
Chapter II of this thesis. (Larsen, 1987) 

The remaining two theses of the seven were peripheral studies which applied alter- 
nate methods to attrition estimation. Hogan attempted multi-year forecasting using 
exponential smoothing; the smoothing constants were rather unusual and extreme and 
his results inconsistent (Hogan, 1986). Yacin applied logistic regression in the attempt 
to develop an attrition rate scheme; the only new results were the identification of some 
areas that exhibited similar attrition behavior (Yacin, 1987). 

This thesis is the first to integrate the two main areas of study. Whereas previous 
studies of shrinkage methods have used ad hoc aggregation schemes, we now combine 
the implementation of a defensible aggregation method with empirical Bayes estimators. 
Moreover, these are applied to a larger and more refined data base. The results have 
been quite promising in that we have achieved greater stability in attrition rate esti- 
mation; we have defined guidelines for a heuristically appealing aggregation scheme; and 
we have acquired an increased understanding of the data base and developed more effi- 
cient ways to use it. 

C. ORGANIZATION 

The remainder of this introductory chapter provides a more detailed description of 
the small cell problem and the data base. The aggregation problem is discussed and the 
proposed aggregation method is presented in Chapter II. 

The shrinkage estimation methods, generally classified as empirical Bayes type esti- 
mators, are described in Chapter III. Several variations are presented to allow 
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comparison and to gain further insight into their performance. Testing of these methods 
is important but for practical purposes must be carried out using samphng methods. 
The rationale used to select test cases, the cross validation techniques, and the measures 
of effectiveness used to evaluate the results are discussed in Chapter IV. A discussion 
of the results of the cross validation is also included in this chapter. 

Finally, conclusions and recommendations based on these results are contained in 
Chapter V. 

D. SMALL CELL PROBLEM 

Marine Corps officers can be classified and thus partitioned by several attributes. 
The major partitioning of officers is by grade, years commissioned service (YCS), and 
military occupational specialty (MOS). Grade describes the position an officer holds in 
the service. The numbers of officers in each higher grade have a pyramid structure, i.e., 
there are more officers in the lower grades than in the higher grades. YCS is the total 
number of years served since becoming a commissioned officer. There is a strong cor- 
relation between grade and YCS since an officer generally moves up the grade structure 
as he gains in YCS. .MOS is a four-digit code identifying the specific skill for which a 
.Vlarine is trained. .MOS need not remain constant over an officer's career, although 
most changes in .MOS occur in the early years of commissioned service. An officer has 
a single primary .VI OS, however as he develops new job skills he may be assigned one 
or more additional .MOSs. 

For many purposes, partitioning by grade, YCS and MOS is sufficient. However in 
some applications additional refinement by service component, commissioning source, 
sex, race or education level may be necessary. Service component consists of three cat- 
egories; regular officers, reserve officers, and reserves who have augmented to become 
regulars. It is strongly correlated with commissioning source, i.e., an officer receives a 
regular or reserve commission depending upon the commissioning source. Both affect 
an officer's initial service obligation, which is generally three to five years (except 
aviators, whose obligation is dependent upon the amount of flight training). Officers 
who receive a reserve commission normally serve three to four years active duty (except 
aviators), by the end of which they must have either augmented into the regular force 
or are then separated from active service. 

These cross-classifications may be viewed as breaking the officer population into a 
multidimensional array, with each specific intersection of the classifications called a cell. 
The total number of possible cells is quite large, on the order of 10‘. .Vlany of these cells 
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are structurally infeasible in that no officer could possibly fit the cell characteristics, e.g., 
there are no Majors with two years commissioned service. The total officer inventory 
of approximately 20,000 officers is partitioned by the remaining feasible cells; some cells 
have inventory' as large as 50, however most have less than five. An officer's character- 
istics are dynamic, i.e., as an officer moves through the grade, YCS and MOS structure 
he moves from one cell to another. As a result, the inventory of the feasible cells is also 
dynamic and fluctuates between zero and low inventory (less than five) over time. 

These sparsely populated cells have very unstable empirical attrition rates. For ex- 
ample, a cell whose inventory' is two officers of which one leaves the service during a 
given time period yields a 50% empirical attrition rate, whereas a cell whose inventory 
is one officer who remains in service during the same time period yields a 0% attrition 
rate. It is obvious that neither of these empirical estimations provides a usable attrition 
rate. Furthermore, these two rates could change dramatically during the next time pe- 
riod, typifying their instability. 

Even when more modern estimation techniques (e.g. shrinkage) are applied, these 
small cells can still create statistical instability, thereby producing intolerably variable 
attrition rate estimates. The problem then is how to deal with these low inventory' cells, 
or “small cells” in order to achieve stability. 

E. DATA BASE 

In this thesis we benefit from a refined data base compiled by NPRDC and made 
available to the Naval Postgraduate School in 1987. This data base, used by Larsen in 
his aggregation work (Larsen, 1987), was not available for the previous estimation 
studies at NPS. 

The new data base provides more detailed information about the officer population. 
The grade structure now allows separation of Limited Duty Officers (LDO) as well as 
Warrant Officers (WO) from unrestricted officers. Officers who have failed selection to 
the next higher grade can also be identified. YCS is listed instead of length of service 
(LOS), which became ambiguous when dealing with officers who have prior enlisted 
service. MOS can now be broken out completely into 236 .MOSs or summarized by the 
39 occupational fields. Service component and commissioning source are both new 
categories. Other new categories that are not considered here are education level, race, 
additional MOSs, and military schools completed. Larsen gives a complete description 
of the classifications (Larsen, 1987, pp. 66-82). 
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The data base also allows attrition to be broken out by retirement, release, dis- 
charge, resignation, etc., but for our immediate purpose we are only concerned with the 
total number of losses for any reason. 

This refinement of the data presents a dichotomy: we can now break the data into 
more definitive cells to search for homogeneous attrition behavior and stability in esti- 
mation, but this leads to an even greater number of low inventory cells. 

The new data base contains ten years of inventory' and attrition data from the period 
1977-1986, a significant improvement from the previous seven year data base covering 
the period 1977-1983. The inventory data is now obtained from quarterly vice yearly 
snap-shots of the officer population. The attrition data is annualized, i.e., the attrition 
count for a cell reflects the number of personnel who leave the service at any time during 
the year. Attritions are credited to the cell which the officer occupies at the time he 
leaves. 

Two problems arise from this quarterly versus annualized data. First, it is possible 
for a cell to record zero inventory' via the snap-shots, yet be credited with one or more 
attritions. To avoid this situation, the cell inventory used in all calculations is defined 
to be the larger of the inventorv and the attrition count. This ensures that the inventorv 
for a cell is as least as large as its recorded number of leavers. (This override occurs in- 
frequently; a more sophisticated treatment would require significant model enhance- 
ment.) Second, to use the inventory' and attrition data together we must divide the 
inventory' data by four. This poses a philosophical problem when invoking a binomial 
model: the sample size may not be integral. However, for our application the usual 
mean and variance formulas are usable and can still serve in the interpolative sense. 
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II. CELL AGGREGATION 



A. GENERAL 

The aggregation problem takes on new meaning with the use of shrinkage estima- 
tors. Originally, aggregation had only one concern: how to pool cells together into a 
single cell in order to meet a user-defined minimum inventory threshold. This single 
aggregated cell was then used to determine the attrition rate estimate for the original, 
unaggregated cell. In this way an estimated value for a cell is obtained by using the 
grand mean for many cells. 

The empirical Bayes multiparameter estimation techniques provide a way to com- 
promise, using both the stability of a grand mean and the specific information of an in- 
dividual cell. Now we pool cells together and obtain a number of cells that meet the 
user-defined minimum inventoiy threshold. It is important to note that we should be 
able to use a lower inventory threshold with empirical Bayes, thus retaining individual 
cell behavior to a greater extent. It is also important to use cells with homogeneous 
attrition behavior in the aggregation process. 

B. BACKGROUND 

The aggregation method currently used by .MCORP is called the Small Cell Override 
Methodology (N'PRDC, 1985, Appendix H). It is used to solve the original aggregation 
problem, i.e., if a cell is below the user-defined threshold, then cells are adjoined to the 
original cell until the threshold is met. The process for selecting cells for adjunction is 
rather crude, and large-scale with only a few levels (prior to using the entire officer 
corps). The attrition rate estimate for the original cell is the empirical rate from this 
aggregated cell. 

To begin the process, the user defines a cell for which an attrition rate estimate is 
required by grade, YCS and .MOS. The user also defines the minimum cell inventoiy 
threshold (and other parameters which are not relevant here). If the cell he identifies 
meets the threshold, no aggregation is required and the empirical attrition rate is deter- 
mined. If the cell is below the threshold, additional cells must be added until the 
threshold is met. 

This search for additional cells occurs by expanding by YCS and MOS, with grade 
remaining fixed throughout. Expanding in this sense means changing the YCS or MOS 
parameter to identify the additional cells to be added to the original cell. Initially, the 
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single cell is expanded by YCS. For example, if the original cell's grade/ YCS/M OS was 
Capt/7/0802, the cells identified by Capt/6/0802, Capt/8/0802, etc., are added sequen- 
tially until the threshold is met. This YCS expansion has an upper bound at the 20 YCS 
point; an obvious boundary' for attrition behavior due to retirement eligibility. If the 
original cell's YCS is above 20, then 20 would serve as the lower YCS bound. 

If the threshold is not met after the YCS expansion, the override method starts over 
with the original cell and expands by MOS. Each MOS belongs to one of nine MOS 
groups which are defined along traditional Marine Corps functional areas, e.g., all heli- 
copter pilot MOSs are grouped together as are all combat support MOSs. MOS ex- 
pansion adds those cells identified by the MOSs in the same MOS group as the MOS 
of the original cell for the original YCS and grade. If the threshold is not met, all the 
MOSs in the MOS group are expanded by YCS in the same manner as the YCS expan- 
sion discussed previously. 

If this MOS group and YCS expansion is unsuccessful, the override method starts 
over with the original cell and expands by all .MOSs for the original grade and YCS. If 
necessary, all the MOSs are expanded by YCS as before. 

Cell aggregation using this expansion method can potentially include all MOSs and 
YCS bounded only at the 20 year point. The desire to aggregate using cells with ho- 
mogeneous attrition behavior is obviously compromised. Larsen provides a more com- 
prehensive description of the current method (Larsen, 1987, pp. 16-22). 

Larsen examined attrition behavior in the MOS and YCS structure. He applied a 
hierarchical clustering algorithm in an attempt to find MOSs and YCSs that displayed 
homogeneous attrition behavior. He confirmed the belief that YCS is an important 
factor. The YCS expansion bounds he proposed reflect points at which officers reach 
the end of their initial service obligation as well as when they are eligible for retirement, 
which makes them especially viable from an intuitive point of view. Larsen also found 
that some MOSs did not cluster strictly by functional areas. This was especially signif- 
icant in the aviation community. Whereas the previous data base allowed aviators to 
be considered only as one occupational field, the refined MOS information was able to 
identify six distinct homogeneous groups of aviators. 

Larsen uses these results to define more refined .MOS groups and YCS boundaries. 
To avoid the giant expansion leap from MOS group to all MOSs, he proposed a hier- 
archy of small .MOS groups, large .MOS groups and major .MOS groups developed by 
observing which MOSs tend to exhibit similar attrition behavior. Homogeneity is 
greatest within the .MOS group, and becomes successively worse as we move to the large 
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MOS group and then the major MOS group. Each MOS is assigned to a small MOS 
group. Small MOS groups combine to make a large MOS group, and large MOS groups 
combine to make a major MOS group. 

Each small MOS group is assigned a set of YCS expansion bounds. Due to the 
different attrition behavior of the small MOS groups with respect to YCS, three different 
sets of YCS expansion bounds are proposed. 

Initial expansion is by YCS within the specified boundaries, with grade and MOS 
held constant. If more expansion is required, we retain this aggregated cell and expand 
by small MOS group for the original grade and YCS. If the aggregated cell is still below 
the threshold, the MOSs in the small MOS group are expanded by YCS. Subsequent 
expansion to large MOS group and YCS, and major MOS group and YCS is accom- 
plished until the threshold cell inventor}' is met. 

Unlike the current expansion method, expansion using Larsen's proposed method 
will not cross defined MOS groups or YCS bounds to ever include all MOSs and YCSs 
bounded only at the 20 year point. Larsen provides a more detailed description of his 
recommended expansion rules (Larsen, 1987, pp. 45-61). 

C. EXPANSION METHOD 

We now address the methods used to obtain the cells required for use with empirical 
Bayes estimation techniques. Expansion continues to mean finding more cells to be 
used, however we no longer simply add these cells to the original cell to form a single 
aggregated cell. The cells identified by the expansion process are now aggregated to- 
gether to produce a number of cells. After the discussion of the expansion process in 
this section, an actual aggregation scheme is introduced in the next section. 

To begin the estimation process, the manpower planner defines a specific cell by 
grade, YCS and MOS. The attributes service component and commissioning source are 
also included as possible cell descriptors for the purposes of this study. All other de- 
scriptors listed in the data base--sex, education level, additional MOSs, race and military 
schools— are ignored. Loss types are considered as a combined total, i.e., in this study 
we do not discriminate among the various types of losses. The first three user-defined 
descriptors— grade, YCS and MOS— are single-value inputs. The last two descriptors, 
service component and commissioning source, can be single values, or either one of them 
can be treated as a vector of values for each single cell. This vector is collapsed (total 
the components) during the aggregation process, i.e., all records which meet any of the 
vector's values are included in the same cell. As in the previously described expansion 



9 



methods, only YCS and MOS change during expansion, the remaining cell descriptors 
remain constant. 

To use shrinkage techniques, the amount of expansion required not ‘only depends 
upon the minimum cell inventor}’ threshold but also upon a new input parameter: the 
threshold number of cells. These two parameters are denoted: 

1. To - cell inventory threshold. The minimum average inventory for a cell obtained 
by averaging the cell inventory over the ten years of data. 

2. Ko - threshold number of cells. The minimum number of aggregated cells whose 
inventory exceeds Tq. These aggregated cells are the input cells for the empirical 
Bayes techniques. 

For example, if To = 5.0 and Ko = 10, the expansion algorithm continues until at 
least ten aggregated cells, each with average inventory 5.0 or larger are obtained. Since 
we are concerned primarily with the small cell problem the values of To and Ko used are 
selected to range from five to 30. It is also presumed that To is less than or equal to 
Ko . These threshold values can certainly exceed 30 for other applications, however the 
resulting cells are not considered small and their attrition behavior most likely would not 
be as unstable, therefore not requiring special attention. 

Prior to explaining the expansion process, we first define the MOS groups and YCS 
bounds. We have adopted much of Larsen's work in this area; many of the changes are 
minor but are necessary for implementation purposes. 

The general idea of a hierarchy of MOS groups is repeated, as shown in Table 1. 
Each MOS belongs to a small MOS group, a large MOS group and a major MOS group. 
Listed are 14 small MOS groups, which combine to make six large MOS groups, which 
combine to make four major MOS groups. For example, small MOS groups one and 
two form large MOS group one, and small MOS groups three through six form large 
MOS group two. Large MOS groups one and two, which collectively contain small 
•MOS groups one through six, make up major MOS group one. Major MOS group one 
contains only ground MOSs, and major MOS group two contains only aviation MOSs. 
Major MOS groups three and four are special cases as discussed below. 

A subjective decision was made to keep the ground MOSs in groups defined along 
the more traditional functional areas. This is reflected in small MOS groups one through 
six. For estimation purposes it is advantageous if the cell inventories are not too vari- 
able in size (Carter and Rolph, 1974, p.882). It is also desirable to avoid having too 
many MOSs in each small MOS group. This allows the expansion to occur more grad- 
ually, and is especially important for small values of To and A',,. As a result, MOS 0302 
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Table 1. MOS GROUPS 



Group Name 


MOSs 


Small 

MOS 

Group 


Large 

MOS 

Group 


Major 

MOS 

Group 


Combat 


0302 


1 


1 


1 


Combat Support 


0802 1302 1802 1803 


2 


Combat Service 1 


0180 0202 2502 2602 


3 


2 

o 


Combat Service 2 


3415 4002 4302 5803 


4 


Combat Logistics 


0402 3002 3060 3502 6002 


5 


Air Control 


7204 7208 7210 7320 


6 


Fixed Wing Pilots 


7501 7511 7522 7542 7543 
7545 7576 


7 


3 


2 


F-18 Pilots 


7521 7523 


8 


Rotary Wing Pilots + 


7556 7557 7562 7564 7565 
7566 7587 


9 


4 


Naval Flight Officers + 


7508 7509 7563 7581 7583 
7584 7585 7586 7588 


10 


Basic Ground 


0101 0201 0301 0401 0801 
1301 1801 2501 2601 3001 
3401 3501 4001 4301 4401 
5801 6001 7201 7301 9901 


11 


5 


3 


Student Aviators 


7580 7597 7598 7599 


12 


Basic Pilots 


7500 7510 7520 7540 7550 
7560 7575 


13 


Lawxers 


4402 


14 


6 


4 



(infantry) is placed alone in a small MOS group. This MOS contains approximately 
15% of the total officer population, and therefore its respective cells normally contain 
large inventory. The MOSs in small MOS group two also contain fairly large inventor}', 
therefore are grouped together and their first expansion is with MOS 0302. The 
remaining ground MOSs in small MOS groups three through six have similarly small 
inventory. 

The aviation small MOS groups (seven through ten) remain relatively unchanged 
from Larsen's recommendations. MOS 7564 (CH-53 pilot), was removed from a ground 
MOS group and added to small MOS group nine, which reflects its functional area. 
MOSs 7551 (C-9 pilot), 7552 (TC-4C pilot), 7555 (UC-12B pilot) and 7559 (CT-39 pilot) 
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were deleted since they are not primary MOSs. MOS 7530 (basic pilot VMFA (F-4)) 
was deleted since it is not a current MOS. (MCO P1200.7G, 1988) 

Officers who have not acquired sufiicient schooling or field experience to qualify for 
a primary MOS listed in small MOS groups one through ten are gathered together as 
basic officers or students in small MOS groups 11-13. These officers are generally sec- 
ond lieutenants or junior first lieutenants with three or fewer YCS. They are disregarded 
for the remainder of the study because their attrition rates are extremely low; probably 
because none of the officers in these groups have reached the end of their initial 
obligations. 

MOS 4402 (lawyers) is considered a special case and is not addressed in this study. 

All MOSs listed in Table 1 are primary MOSs for unrestricted officers as listed in 
the current Military Occupational Specialties Manual (.MCO P1200.7G, 1988). It would 
be a logical and relatively simple extension of this table to create additional groups 
containing LDO and WO MOSs. These grades are not considered in this study and 
therefore their respective .VIOSs are excluded from the table. 

Several of these seemingly ad hoc decisions to alter Larsen's recommended MOS 
groups are due to the YCS expansion bounds shown in Table 2. Every effort was made 
to group MOSs with similar YCS expansion bounds to allow for feasible implementation 
of the expansion algorithm. This is especially applicable when expanding to large and 
major MOS groups. 



Table 2. YCS EXPANSION BOUNDS 



MOS Group 


Small 

MOS 

Groups 


Bounded YCS Groups 


Fixed Wing Pilots, F-IS Pilots, 
Lawyers 


7, 8, 14 


(1-6, 8-19) (7) (20-25) (26) 


Rotary Wing Pilots, Naval Flight 
Officers 


9, 10 


(1-5,8-19) (6,7) (20-25) (26) 


All Others 


1-6, 11-13 


(1-3,6-19) (4,5) (20-25) (26) 



The YCS expansion bounds reflect the maximum expansion allowed from the initial 
YCS defined by the user. For example, if the original cell's grade/ YCS/.M OS is 
Capt/9/7501, we see from Table 1 that this MOS belongs to small .MOS group seven. 
Thus its YCS expansion bounds are listed on the first line of Table 2. The value of nine 



12 



for YCS falls in the first YCS range, thus we could expand using all YCSs from one 
through 19, excluding seven. If the YCS for this original cell had been seven, no YCS 
expansion would be allowed. 

These YCS expansion bounds are used with the MOS groups to define the addi- 
tional cells which can be used with the original cell to obtain the required number of 
cells, K^, each with minimum average inventory, Tq. The expansion stages are: 

1. Stage 1 - Locate the small MOS group which contains the user-defmed MOS. The 
initial cells are those specified by the MOSs in this group for the user-defined YCS, 
grade, service component and commissioning source (grade, service component and 
commissioning source remain fixed throughout the expansion process and thus are 
not repeated). These cells are aggregated to obtain cells with average inventory 
greater than or equal to T^. After aggregation, if the number of cells is greater than 
Kfi, stop, otherwise go to Stage 2. 

2. Stage 2 - Expand by incrementing YCS (YCS-1, YCS+1, YCS-2, YCS + 2, etc.) 
within the YCS bounds listed in Table 2 for all .MOSs in the small MOS group. 
After each YCS increment, aggregate the cells to obtain cells with average inven- 
tory greater than or equal to Tg. After aggregation, if the number of cells is greater 
than Kf„ stop, otherwise continue to increment by YCS. If the YCS bounds are 
reached before obtaining enough aggregated cells, retain the cells identified in 
Stages 1 and 2 and go to Stage 3. 

3. Stage 3 - Expand to the large MOS group for the single user-defined YCS. Ag- 
gregate the cells to obtain cells with average inventoiy greater than or equal to Tq. 
After aggregation, if the number of cells is greater than K,,, stop, otherwise go to 
Stage 4. 

4. Stage 4 - Expand by incrementing YCS for the large .MOS group. After each YCS 
increment, aggregate the cells to obtain cells with average inventory greater than 
or equal to T^. After aggregation, if the number of cells is greater than K,,, stop, 
otherwise continue to increment by YCS. If the YCS bounds are reached before 
obtaining enough aggregated cells, retain the cells identified in Stages 1 through 4 
and go to Stage 5. 

5. Stage 5 - Expand to the major .MOS group for the single user-defined YCS. Ag- 
gregate the cells to obtain cells with average inventoiy greater than or equal to Tq. 
After aggregation, if the number of cells is greater than K^, stop, otherwise go to 
Stage 6. 

6. Stage 6 - Expand by incrementing YCS for the major .MOS group. After each YCS 
increment, aggregate the cells to obtain cells with average inventoiy greater than 
or equal to Tg. After aggregation, if the number of cells is greater than K^, stop, 
otherwise continue to increment by Y'CS. If the YCS bounds are reached before 
obtaining enough aggregated cells, stop. No more expansion is allowed. Inform 
the user that the thresholds are unattainable. Do not cross any major MOS group 
or YCS bounds. 

Tw'o points about the expansion process are emphasized. First, we retain the cells 
identified by all previous stages as we progress to the next stage. As stated before, the 
degree of homogeneity decreases as we move from small to large to major .MOS groups. 
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Thus we want to locate as many cells from the small MOS group as possible before we 
expand to the large MOS group, and then locate as many cells from the large MOS 
group as possible before expanding to the major MOS group. The YCS expansion for 
each group may be different, e.g., the small MOS group may be expanded by all YCSs 
within the given YCS range, but the large MOS group may only be expanded by a few 
YCSs before the thresholds are attained. 

The second point is that, when aggregating cells, any aggregation that was per- 
formed previously is discarded and all cells currently identified are pooled and made 
available for aggregation. This affords the aggregation algorithm greater flexibility and 
could create more aggregated cells than if the aggregated cells from previous stages were 
left intact, thereby keeping the amount of expansion to a minimum. 

The aviation small MOS groups (seven through ten) make up the only major MOS 
group (two) that contains different YCS bounds, i.e., small MOS groups seven and eight 
have different YCS expansion with regard to year six than do small MOS groups nine 
and ten. To implement the expansion algorithm in a computer program, this difference 
is overcome by using the YCS bounds for the original user-defined MOS. For example, 
suppose MOS 7501 from small MOS group seven is the original MOS. If MOS expan- 
sion continues into the major MOS group, the MOSs in small MOS groups nine and ten 
would follow small MOS group seven's YCS expansion bounds. 

In summar}', this method of grouping MOSs should provide greater homogeneity 
among cells which are used in estimating attrition rates. Unlike the current method, 
ground and aviation MOSs are never used together. The YCS bounds provide a logical 
and effective way to treat periods of different attrition behavior. However, the greater 
the expansion the less homogeneous the cells become, which should be kept foremost in 
mind when setting the threshold parameters. 

D. AGGREGATION METHOD 

While the expansion steps are being undertaken in order to achieve the threshold 
levels specified by the user, those cells with inventory less than must be gathered up 
into larger, aggregated cells whose combined inventory exceeds Tq. In order to limit the 
expansion to as few additional MOSs and YCSs as possible, we desire to maximize the 
number of aggregated cells obtained at any stage of the expansion. 

The term maximization suggests the possible use of linear programming (LP). While 
an LP would ensure maximization, this would not be a trivial problem to solve, i.e., the 
LP relaxation would almost certainly fractionate cells, using their inventory' in more than 
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one aggregated cell. This is not allowed since a cell may be assigned intact to only one 
aggregated cell. Thus an integer LP would be required which would typically contain 
500 or more integer variables. This method would not be expedient in terms of computer 
usage, especially considering the potential number of integer LPs that may have to be 
solved for a single estimation cycle. 

While we are trying to maximize the number of aggregated cells, it would be satis- 
factory^ to obtain close to the maximum if we could preclude the expense in computer 
time required by an integer LP. For this reason, a heuristic “greedy” algorithm was de- 
veloped. Complete descriptions of the heuristic algorithm and the LP formulation are 
contained in Appendix A. The performance of this heuristic is discussed along with the 
results of the empirical Bayes methods in Chapter IV. 
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III. ESTIMATION METHODS 



A. GENERAL 

Once the cell aggregation phase is completed, we begin the attrition rate estimation 
process. The following notation is used to define the cell data 

K = number of cells ... 

T = number of years of data. ^ ' 

Then for / = 1, ..., K and / = 1, ..., T 

(t) = inventorv of cell / in vear t 

' ( 2 ) 

Yj (/) = number of attritions in cell i in year i. 

The cell data is assumed to be independent binomial, i.e., T, (/) ~ Bin{N^{i), p,) . A 
success is defined to be an attrition, i.e., an officer from that cell leaves the service during 
the year. The empirical attrition rate for cell i is given by the Maximum Likelihood 
Estimator (MLE) 



E 



A / 




/ 



( 3 ) 



This estimate of p works well for cells with large inventory, but not those with small in- 
ventory, which is most often the case in our application. 

The MLE has been shown to be dominated by shrinkage methods for K>2 (Carter 
and Rolph, 1974; Efron and .Morris, 1975; Casella, 1985). These methods find a grand 
mean or central attrition rate for the group of cells and a shrinkage factor for each cell. 
Previous theses have primarily used a common shrinkage factor for all cells; w'e now' al- 
low this shrinkage factor to vary from cell to cell. Each cell's MLE is shrunk towards 
the central rate by its shrinkage factor. In this w^ay, attrition information from one cell 
“spills over” into other cells. 

The shrinkage methods are developed under the theoretical assumption that the data 
is normally distributed. Most of the previous studies using empirical Bayes methods 
have used independent normal data with constant variance (Efron and Morris, 1972, 
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1973, 1975; Dickinson, 1988). Some applications have used binomial data, using a 
transformation to make it behave more like normal data (Carter and Rolph, 1974; Efron 
and Morris, 1975). In Carter and Rolph's estimation of fire alarm probabilities, trans- 
formation of the binomial data did not have a large effect on the results (Carter and 
Rolph, 1974). Our application allows us to investigate the impact of the transformation 
when applied with more extreme values of p. 

Six variations of the empirical Bayes method are applied to the attrition rate esti- 
mation problem. The first four are similar in that they use the same iterative procedure 
to compute the amount of shrinkage for each cell. Of these four, two are on the trans- 
formed scale and two on the original scale. Each scale includes two methods of com- 
puting the cell variance: one method where the variance is time dependent and the other 
where it is time independent. The two variance calculations, if they produce like results, 
provide supporting evidence for the assumption that the data is independent and iden- 
tically distributed over time. This assumption is certainly questionable, since an officer 
who remains in a given MOS will move through the YCS and grade cell structure in a 
predictable manner. As a result, variance that is constant in time (time independent) 
may not perform as well as one that allows for time variation. The fifth method uses a 
different iterative procedure to determine the amount of shrinkage and is addressed 
separately in paragraph III.C.. The final method breaks the cell data into its vector 
components (service component or commissioning source) before shrinkage techniques 
are applied and is addressed in paragraph III.D.. 

B. EMPIRICAL BAYES 
1. Transformed Scale 

We begin our application of empirical Bayes methods on the transformed scale 
in an effort to overcome some of the weaknesses in our assumptions. The transforma- 
tion we use is the Freeman-Tukey transform, a modification of the basic arcsin trans- 
formation for binomial data. Its purpose is to stabilize the variance at one and make 
the data behave more like normal random variables. The form used is 



Xt{t) = yV'A^(/) + .5 






2F,(/) 

A'/ (0+1 



-I- arcsim 2 



>/(0+i 

a;w + i 




(4) 
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Now, let 



X (/') 

XTt{t) = , for t= (5) 

7a',W + -5 

except when N, (r) = 0 (no inventory in year /), in which case XT, (r) does not exist and 
we reduce T, by one. The time average of the transformed values for cell / becomes 

XTB, = ( 6 ) 

h V 

We now need to compute the variance of these time averages. Two methods 
are used; the first calculation is time dependent, i.e., the variance changes over time, the 
second is time independent. 

The transform stabilizes the variance at one for large values of n and non- 
extreme values of p. These requirements on n and p are often violated in our application, 
therefore we have many combinations of n and p for which the variance is less than one. 
Dickinson was able to discover an interpolative formula which provides a good approx- 
imation for the variance of the transformed values, X, (r), for small values of n and p, and 
K>3 (Dickinson, 1988, pp.8-11). This variance is given by 

VariX.iO) = min{l,V(Xim ( 7 ) 

where V(X, (i)) is found by solving 

V{X, (0) = u {X, (/) + q'’* {X, (/) + C - 1)'’^ (8) 

with 

C = ViV,(/)-F.5(-^) (9) 

and 

a = 1.6835 6, = -.8934 = .9881 . (10) 

Equation (8) obviously breaks down if >T,(/)-t-C<l. When this occurs, we set 
X, (i) + C = 1.001 and continue. The effect is to use a small but positive variance. The 
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value of one in Equation (7) dominates for about X, (/) + C ^ 2.2. The variance of the 
time average is then 



VariXTB^ = -^Y.^ar{XT,{i)) = 



r; 



VariXiii)) 
Bi'i (/) + .5 



( 11 ) 



The second method of computing the variance is the more familiar one. Con- 
tinuing from Equation (6), the variance of the transformed values is given by 



( 12 ) 

^ r 

The variance of their average is therefore 

Var{XTBi) = Var{XT,) . (13) 



Regardless of which variance calculation we use, the same iterative algorithm is 
used to determine the empirical Bayes estimate for each cell. This estimate, XEB,, is 
found by solving 



XEBi 



A 

A + Var{XTBi) 



XTBi + 



VariXTBi) 

A -t- Var(XTBi) 



XBB 



(14) 



where XEB/, XTB^ and Var{XTB) are cell specific, XBB is the (weighted) grand mean or 
central attrition rate, and A is the variance of the prior distribution of the cell means. 
These latter two values must be estimated simultaneously using the following iterative 
algorithm. 

We initialize the algorithm with A = 0 and store the previous value of A by 



^0 



^A. 



Now compute the (w-eighted) grand mean, XBB. Let 



(15) 



1 

A + Var{XTBi) 



(16) 
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and 



Vi 



K 




7=1 



Then 



K 

XBB = Y^yiXTBi 
1=1 



The updated value of A becomes 



A 



A - 



K 

A" - 1 - ^ OLiiXTBi - XBBf 

£=1 

K 

^ o^iiXTBi - XBBf 

i=l 



(17) 



(18) 



(19) 



If < 0, set = 0 and exit. This represents the case when there is 100% shrinkage to- 
ward the grand mean. If .4 > 0, then check \ A — Aq \ < e. (e.g., c = .0001). If false, return 
to Equation (15) for another iteration. If true, the iterations have converged. Exit with 
the current values of A and XBB for use in Equation ( 14) to solve for the XEB,. 

Close study of Equation (14) shows that the amount of shrinkage changes from 
cell to cell since the variance terms are generally not equal. Specifically, cells with higher 
variance are shrunk more than those with lower variance. In addition, if A is small the 
shrinkage is greater towards XBB. As A -* oo, the shrinkage is minimal and the indi- 
vidual cell means dominate. 

Once the XEB, are determined, these values must be transformed back to the 
original scale. We use 



A- = y{l-f sin(A£i5,)}. 



( 20 ) 



2. Original Scale 

We return to the assumption of binomial data for original scale calculations. 
As in the transformed scale, two methods to calculate the variance are used. We begin 
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with 



A Yi (0 

xm = Am = ^ . 



( 21 ) 



As before, if A, (r) = 0 (no inventory in year t), XT, (t) does not exist and we reduce T, 
by one. This leads to the time average for cell / as 

= ^EAM. (22) 

' / ' r 

The variance calculation which is time dependent, i.e., changes over time, is 

given by 



Var(XTBi) 



-^ 2 ] Var(xr,(,)) 

h t 




Pirn -pm 
N;{t) 



(23) 



We return to Equation (15) with these variance values to perform the iterative 
algorithm for finding the empirical Bayes estimate, XEB,, given by Equation (14). Since 
we are already in the original scale, the transformation given in Equation (20) is ignored, 
U.,p, = XEB, 

A problem arises while performing the iterations if a cell has Y,{t) = 0 'it (zero 
attrition for every year). In this case, the variance given by Equation (23) equals zero. 
When this value is used in Equation (16), the formula for a, becomes undefined. We 
resolve this problem using the Laplace Law of Succession. Assume that 
T,(r)-5//7(A^(r),A)andlet 






and 



<h 



A/ it) - Y; (/) 

A)W+1 



(24) 



be the estimates as prescribed by this law, i.e., Bayes estimator using uniform prior. 
Then 



Var 



Yj{t) 

a; (0 



Pi (h 



Niit) 



/ y,(/) + i Y AK0-r,(r) \ 

\ A/(r)+l j\ A;(r)+1 / 



(25) 
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If Y, it) = 0, then 




A/ (0+1 



A^(0+1 ) 



Niit) 



1 

(A^/(0+l)' ■ 



(26) 



This value is used as the summand in Equation (23) whenever Y, (0 = 0 (zero attrition 
in any year). 

For comparison purposes we again compute an alternate variance which is time 
independent. Continuing from Equation (22), let 



Pi 



/ 

t 



(27) 



The alternate variance is given by 



VariXTB,) = ^ 7;. (,)) = Z ^ ' P*’ 

The problem with cells that have T, (0 = 0 V r (zero attrition for every year) also 
occurs here, since the variance given by Equation (28) would equal zero. Using the same 
concept as before, we obtain the formula 



VariXTB;) 




(29) 



However in this case, this variance formula is necessary only if all years have zero 
attrition. 

As before, we return to Equation (15) to perform the iterative algorithm for 
finding the empirical Bayes estimate, XEB,, given by Equation (14). 
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C. EFRON-MORRIS METHOD 

This method is a modification of the iterative algorithm used to estimate A and 
XBB given by Efron and Morris (Efron and Morris, 1973, pp. 127-129). It differs from 
the method given by Equations (14) through (19) in that it allows the variance of the 
prior. A, to change from cell to cell. It also gives greater weight to the cells with low 
variance, and reduces to the James-Stein estimator when the cell variances are constant. 

Only one scenario for this method is considered, corresponding to the initial trans- 
formed scale, time dependent variance case. Thus, Equations (4) through (11) are re- 
peated, and we begin from the point where we are entering the iterative algorithm. To 
simplify the following equations, let D, = Var{XTB^ as given by Equation (11). 

We initialize the algorithm with A, = 0 and SP, = 0 (previous values of S) for 
i=l,...,K. Let 




(30) 



and 



K 



(31) 




Then 



K 




(32) 



and 



5,. = {XTBi-Xf. 



(33) 



Now set i = 1 and let 




( 34 ) 
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and 



SDi 



y I — . 

Ml 



(35) 



We then use the Newton-Raphson iteration method to solve 



4- = 



iSi-Wi) + {A, + DfSN) 
2 + {Ai+DfSDi 



= ^(4)- 



(36) 



First set A Pi*- A, for i=\,...,K (previous values of A^. The updated value for Ai 
becomes 



^i 



Aj - g(Ai) 

1 -m) ■ 



(37) 



If Ai<0, set A, = 0, let / = /+!, and return to Equation (34). If /1,>0, then test 
IAi — AP,j > e. If true, return to Equation (36). If false, let /=/ + 1 and return to 
Equation (34). 

In either case after incrementing /, if / = 1, we exit and test | S, — SPi | < e V /. 

If false, we update SP, *- Si and return to Equation (30) with updated values of /!,. If 
true, the iterations have converged and we must finalize the estimators, XEMf Let 






M (4 +^/ 



(38) 



and 



Bi 




Di 

Ai + Di 



(39) 



If 5, > 1 , set 5, = 1, or if P, < 0, set 5, = 0. Then 

XEMi = X+{1- Bi){XTBi-^. (40) 

This equation is comparable to Equation (14), which was used to determine the trans- 
formed scale estimates, XEB„ using the previous iterative algorithm. The quantity P, is 
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the amount of shri^a^^^oward the grand mean, X. The corresponding quantity in 
Equation (14) is ^ + Var(XTB) ' 

To obtain p„ the XEM, must be transformed back to the original scale using XEM, 
in place of XEB^ in Equation (20). 

D. VECTOR METHOD 

This method is similar to the previous methods in the sense that it uses the aggre- 
gated cells produced to meet the defined threshold levels. However, prior to the esti- 
mation process, we now partition each aggregated cell by either service component or 
commissioning source, thus obtaining cells whose elements are vectors. The procedure 
given by Efron and Morris, modified to compensate for the assumed variance of the time 
averages, provides the framework for this method (Efron and Morris, 1972, pp.341-344). 

The separation by service component or commissioning source requires us to define 
a third index: the components of the vector. Let 

P = number of service components/commissioning sources. (41) 

Then for / = \,...,K,j= l,...,Pand/= 1,..., T 

Nij (r) = inventory’ of cell / and vector component j in year i 

(0 = number of attritions in cell / and vector component j in year t . 

Before we had K cells with scalar information, but now we need z. Kx P matrix where 

p p 

^/V,,.(r) = A)(r) and y,-;(t) = k;-(t). (43) 

;=i ;=i 

A requirement for this method is that K> P + 2, for reasons that will soon become 
obvious. 

We begin by defining X^ (r) as the transformed value for (r) and 7,^ (r) as given 
by Equation (4). Continuing in similar manner as the previous transformed scale 
methods, let 



XTijit) = 









for t = 1, ..., Ti 



ij 



(44) 



except when N/j (t) = 0, in which case XT,j (r) does not exist and we reduce T^j by one. 
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The time averages of the transformed values are then 



XTBij = -r^yXTijit). (45) 

hj V 

-f 

Here we obtain a vector of grand mean values, with each of the P grand means defined 
by 



K 

XBBj = -^Yj^TBij . (46) 

i=i 

The transformed scale estimate, , is then found by solving 

5ji = XBBj+{Ip-{K- P-2)S~^]{XTBji-XBBj) (47) 

where Ip is the identity matrix of order P, and S“‘ is found as discussed below. Reversal 
of the i j index in this and subsequent equations simply means the transpose of the 
Ky. P matrix to a P x A' matrix. 

To solve for 5"‘, we begin by defining 

S = XjiXji (48) 

where 

Xj, = {XTBji-XBBj)Jv;; . (49) 

The Fy, matrix is the modification required by our application. (The multiplication in 
Equation (49) is element-wise as opposed to normal matrix multiplication.) The Efron 
and Morris method was developed under the assumption that XTB,j~ A'(0,;,1) , 
whereas we are using 



provided that XTij (/) has variance of one. Therefore 
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1 



(51) 



^IJ = 






Nij{i) + .5 



We use the requirement that the P x P matrix resulting from the operations within 
the brackets in Equation (47) must be nonnegative definite to solve for S'* without 
having to actually compute its inverse. We proceed by doing an eigenanalysis of S, 
which is seen be Equation (48) to be a real symmetric matrix. We form the diagonal 
matrix E, which has the eigenvalues, ej, as its diagonal elements, and the matrix E, which 
has the corresponding eigenvectors as its columns. For any ej<{K— P— 2), we replace 
it with the value {K — P —2) . The eigenanalysis provides us with the solution to 



or Sr = FE. 

Post-multiplying by F’^, we obtain 

S = FEF^ 

since F is ortho-normal and therefore F F^ = Ip. We then have 

5"’ = (F E F^)~’ = F E~’ F^ 



(52) 



(53) 



(54) 



which is easily solved since E“‘ is found by replacing the diagonal elements of E by their 
reciprocal. This solution for S’* is then used in Equation (47) to solve for the trans- 
formed scale estimates, . To obtain the attrition rate estimate for a cell, py, , we use 
the inversion formula given by Equation (20) with in place of X£S,. 
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IV. CROSS VALIDATION 



A. GENERAL 

The six estimation methods discussed in Chapter III are evaluated using cross vali- 
dation of the data base. This consists of successively holding out one year's data while 
the other nine years are used to estimate that year's attrition rates. Three measures of 
effectiveness (MOEs) are used to evaluate the validity of our assumptions and the per- 
formance of the estimation methods. Two of these are original scale MOEs-mean ab- 
solute deviation (MAD) and chi square statistic. The third is a transformed scale 
MOE--mean squared error (MSE). Test cases are chosen as input. The results of the 
cross validation are then discussed. 

B. MEASURES OF EFFECTIVENESS 
1. Mean Absolute Deviation 

The MAD is probably the most useful MOE to the manpower planner. Our 
version is augmented to display overestimation and underestimation information. Along 
with the MAD we observe the magnitude of our errors in both directions, which is es- 
pecially important since the cost of overestimating is generally not the same as the cost 
of underestimating. While it does not provide a specific value or standard to gauge the 
performance of our estimation methods, it does provide vejp^ useful insight into tenden- 
cies to consistently underestimate or overestimate. 

For comparison of the estimation methods, we desire a MAD measure that does 
not depend upon cell inventories, yet still displays the overage/underage information. 
For these reasons, we use the attrition rate estimates, p,, as opposed to the estimated 
number of attritions, (p, • A' (0) (''’here t = validation year), in our MAD calculations. 
For those estimates obtained in the transformed scale, the XEBi are inverted back to the 
original scale using Equation (20) prior to calculating this MOE. 

We define the empirical attrition rate for cell i in validation year t as 



Pi 



Yiit) 

AiW 



(55) 



except when A', (r) = 0 (no inventoiy’ in cell / for the validation year). In this case we do 
not compute the cell's deviation from the estimated attrition rate since it would 
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artificially create an overage situation. Therefore, we reduce K by one and continue with 
the remaining cells (the reduced value of K is then used in the following formulas). 

The MAD measures generated for each validation year are 

Ku 

— 77 - = fraction of cells with underage (56) 

K 

where K„ is the number of cells which have underage, 

2 (Pi-Pif 

= fraction of MAD due to underage (57) 



and 



MAD = Ift - ftl . (58) 

/ 

We also calculate the average MAD over the validation years. Here we use a 
weighted average, since the number of cells may have been different in some validation 
years, i.e., a reduced value of K was used in these years. The weighted average takes the 
form 



Y,k,mad, 

Avg MAD = — 

1 




(59) 



where K, is the (possibly reduced) number of cells used in validation year i. 

2. Chi Square 

The chi square test is used as an indicator of how well the binomial model serves 
as a description of the attrition process. The test statistic is 



= Z 

i 






( 60 ) 
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where t is the validation year. As with the MAD calculations, if yV, (/) = 0 we reduce K 
by one and continue. Additionally, if p, = 0 or 1, the denominator equals zero and the 
summand is undefined. The same course of action is used if this occurs—reduce K by 
one and continue. Those estimates obtained in the transformed scale are inverted back 
to the original scale prior to using Equation (60). 

This MOE can be used as a gauge. The chi square statistic given by Equation 
(60) has expected value K and variance 2K. We are looking for a value that is less 
than two standard deviations to the right of the mean, or 

< K-\-2j2K. (61) 

A weighted average chi square is computed in the same manner as the weighted 
average MAD in Equation (59). However, if the number of cells and thus the degrees 
of freedom, K, are different over the validation years a problem arises in determining the 
degrees of freedom for the weighted average. We solve this dilemma by assuming that 
the weighted average chi square has the original value of K degrees of freedom. 

3. Mean Squared Error 

The MSE is used to check the validity of our theoretical basis. It is the average 
squared deviation of the estimated rate from the actual rate, both rates on the trans- 
formed scale. The actual rate used is the transformed validation year data. The MSE 
is defined as 



L(5,m) = ^ ((5; - /i;)^ (62) 

where 

Si = XEBi 

Hi = XTi (r) (/ = validation year) . 

Again, if iV, (r) = 0, we reduce K by one and continue. A weighted average MSE is also 
computed similar to Equation (59). 

The MSE also has a standard to gauge our model. Using Equations (5) and (12) 
we can compute a baseline variance for any given validation cell. The MSE for that cell, 
when compared to the baseline value, provides a figure to gauge the value of using 
shrinkage estimators instead of the cell averages, AT5, . There is considerable variability 
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in these ratios, ranging from 20% to 100%, but 80% appears to be a fair median figure. 
For example, for cell variances computed from Equation (12) running about 0.15, the 
MSE hovers about 0.12. 

4. Vector Method MOEs 

The MOEs discussed above require slight modification before being applied to 
the vector method described in paragraph III.D.. Recall that this method uses K cells 
with service component or commissioning source broken out into a vector of length P. 
An attrition rate estimate, is obtained for each of the Kx. P matrix components. 
Thus now we have KP estimates which are compared to the corresponding empirical 
rates for the validation year. Equations (55) through (63) are modified by replacing all 
i subscripts with ij, replacing all summations over / by double summations over / and J, 
and replacing all instances of K by the product KP. 

C. TEST CASES 

The selection of test cases takes on great importance since they provide the foun- 
dation for comparison of these methods. It would be impossible to test every permuta- 
tion of input parameters; therefore we seek a representative fraction of these which 
would give an accurate account of the performance of our aggregation and estimation 
methods. Because we are using a different data base from previous theses on estimation 
methods, no attempt to duplicate their test cases was made. 

An approach based upon Latin square experimental design principles was used to 
select 30 test cases for the first five estimation methods. The test cases for the vector 
estimation method are addressed later. In determining the test cases, we randomized 
when possible and intervened to force pairings only when necessary'. To begin, we se- 
lected values for the input parameters- r„, K^, grade, YCS and MOS. Service compo- 
nent and commissioning source are ignored for these test cases, i.e., all classifications 
of both are accepted. 

To ensure proper representation from small MOS groups one through ten, one MOS 
from each group was randomly selected: 0302, 1802, 2502, 4002, 3060, 7204, 7545, 7523, 
7557 and 7563. Since YCS and grade are strongly correlated, these parameters were se- 
lected jointly. To ensure each YCS range within the bounded YCS groups was repres- 
ented along with a fair representation of grades, four grade/YCS pairs were selected: 
lLt/4 YCS, Capt/7 YCS, LtCol/20 YCS and LtCol(failed select)/26 YCS. The two 
threshold parameters were also selected jointly, resulting in ten pairs (TJKq): 30.0/30, 
20.0/30, 20.0/20, 10.0/30, 10.0/20, 10.0/10, 5.0/30, 5.0/20, 5.0/10 and 5.0/5. 
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With these choices in place, it was necessary to combine them to define the actual 
test cases. It was decided to limit the grade/ YCS pairs for these cases to lLt/4 YCS and 
Capt/7 YCS due to the large values for the first five threshold pairs. With ten MOSs 
specified, we sought ten test cases. Thus, each of the first five threshold pairs was listed 
twice. Each of the four aviation MOSs was randomly assigned to one of the five 
threshold pairs; the six ground MOSs were then randomly assigned to the remaining 
pairs. The two grade/YCS pairs were then randomly assigned within a set of common 
threshold pairs, e.g., for the two cases with jK^ of 30.0/30, one was randomly assigned 
lLt/4 YCS, the other was then assigned Capt/7 YCS. 

All four grade/YCS pairs would be used with the five remaining threshold pairs. 
Thus 20 more test cases were generated, with each of the five threshold pairs listed four 
times. Each MOS was randomly assigned to two distinct threshold pairs, ensuring that 
the large and major MOS groups were evenly spread throughout the pairs. The four 
grade/YCS pairs were assigned in random order to each set of common threshold pairs, 
ensuring that they were evenly spread across large and major MOS groups. The 30 test 
cases are summarized in Table 3. 

The input parameters for six vector test cases were selected from the 30 test cases: 
Nos. 2, 3, 6, 10, 11 and 20. A small number of vector test cases was initially chosen to 
investigate the possible advantages of the vector method. If this method appeared to 
be favorable, then further testing would be conducted. 

The six test cases contain a cross section of the input parameters. They include 
three ground and three aviation MOSs, and use each of the first three grade/YCS pairs 
twice. The grade/YCS pair of LtCol(FS)/26 YCS was not used because of its extremely 
low inventory numbers, which when broken out into a vector would have been of little 
exploratory use. These cases also include six different \K^ pairs. 

Each of the vector test cases is used twice: first with service component and then 
with commissioning source as the vector component. All three service components-- 
regular, augmented regular and reserve-were used as vector components. Rather than 
use all 15 commissioning sources (these 15 are listed below Table 4) as vector compo- 
nents (many of them would contain little or no inventory') five commissioning sources 
for the ground test cases and five for the aviation test cases were chosen. These five were 
determined to be the sources which contain the largest percentage of inventoiy for the 
respective ground or aviation MOSs. The specific commissioning sources selected along 
with the other vector test case input parameters are summarized in Table 4. 
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Table 3. TEST CASES FOR METHODS 1-5 



No. 


To 


Ac 


MOS 


S : L : M 


YCS 


Grade 


1 


30.0 


30 


0302 


1:1:1 


4 


ILt 


2 


30.0 


30 


7523 


8:3:2 


7 


Capt 


3 


20.0 


30 


3060 


5:2:1 


7 


Capt 


4 


20.0 


30 


7563 


10 : 4 : 2 


4 


ILt 


5 


20.0 


20 


2502 


3:2:1 


7 


Capt 


6 


20.0 


20 


7557 


9:4:2 


4 


ILt 


7 


10.0 


30 


7204 


6:2:1 


4 


ILt 


8 


10.0 


30 


1802 


2:1:1 


7 


Capt 


9 


10.0 


20 


7545 


7:3:2 


7 


Capt 


10 


10.0 


20 


4002 


4:2:1 


4 


ILt 


11 


10.0 


10 


2502 


3:2:1 


20 


LtCol 


12 


10.0 


10 


7557 


9:4:2 


26 


LtCol(FS) 


13 


10.0 


10 


7545 


7:3:2 


7 


Capt 


14 


10.0 


10 


0302 


1:1:1 


4 


ILt 


15 


5.0 


30 


4002 


4:2:1 


4 


ILt 


16 


5.0 


30 


0302 


1:1:1 


20 


LtCol 


17 


5.0 


30 


7204 


6:2:1 


26 


LtCol(FS) 


IS 


5.0 


30 


7563 


10 : 4 : 2 


7 


Capt 


19 


5.0 


20 


3060 


5:2:1 


7 


Capt 


20 


5.0 


20 


7545 


7:3:2 


20 


LtCol 


21 


5.0 


20 


1802 


2:1:1 


26 


LtCol(FS) 


22 


5.0 


20 


7563 


10 : 4 : 2 


4 


ILt 


23 


5.0 


10 


7204 


6:2:1 


20 


LtCol 


24 


5.0 


10 


4002 


4:2:1 


26 


LtCol(FS) 


25 


5.0 


10 


7523 


8:3:2 


4 


ILt 


26 


5.0 


10 


1802 


2:1:1 


7 


Capt 


27 


5.0 


5 


2502 


3:2:1 


7 


Capt 


28 


5.0 


5 


7557 


9:4:2 


20 


LtCol 


29 


5.0 


5 


3060 


5:2:1 


4 


ILt 


30 


5.0 


5 


7523 


8:3:2 


26 


LtCol(FS) 



(S ; L : M = Small MOS Group : Large MOS Group : Major MOS Group) 
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Table 4. TEST CASES FOR VECTOR METHOD 



No. 


To 


^0 


MOS 


YCS 


Grade 


Sendee 

Comp 


Comm 

Source 


2 


30.0 


30 


7523 


7 


Capt 


1 2 3 


(all) 


2 


30.0 


30 


7523 


7 


Capt 


(all) 


1 2 3 5 11 


3 


20.0 


30 


3060 


7 


Capt 


1 2 3 


(all) 


3 


20.0 


30 


3060 


7 


Capt 


(all) 


1 3 7 10 11 


6 


20.0 


20 


7557 


4 


ILt 


1 2 3 


(all) 


6 


20.0 


20 


7557 


4 


ILt 


(all) 


1 2 3 5 11 


10 


10.0 


20 


4002 


4 


ILt 


1 2 3 


(all) 


10 


10.0 


20 


4002 


4 


ILt 


(all) 


1 3 7 10 11 


11 


10.0 


10 


2502 


20 


LtCol 


1 2 3 


(all) 


11 


10.0 


10 


2502 


20 


LtCol 


(all) 


1 3 7 10 11 


20 


5.0 


20 


7545 


20 


LtCol 


1 2 3 


(all) 


20 


5.0 


20 


7545 


20 


LtCol 


(all) 


1 2 3 5 11 



Service Component: 

1 - regular 

2 - augmented regular 

3 - reserve 



Commissioning Sources used: 

1 - U.S. Naval Academy 

2 - Platoon Leader Class-Aviation 

3 - Platoon Leader Class-Ground 

5 - Aviation OlTicer Candidate 

7 - Officer Candidate Course-Ground 

10 - Enlisted Commissioning Program 

11 - N'ROTC-Scholarship 
Commissioning Sources not used: 

4 - Platoon Leader Class-Law 

6 - Marine Aviation Cadet 

8 - Officer Candidate Course- Law 

9 - Officer Candidate Course-Women 

12 - N'ROTC-Ground College 

13 - NROTC-Aviation College 

14 - XESEP 

15 - All Other 
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D. RESULTS 

For each test case, we apply the aggregation method to meet the threshold levels and 
then execute the estimation methods. While the ultimate use of these fhethods is to 
obtain an attrition rate estimate for the original cell, inspection of these estimates would 
be of httle value in evaluating and comparing the estimation methods. Thus the output 
from the program takes the form of the MOEs. 

The inclusion of the entire output from every test case would not only be cumber- 
some but would provide an inadequate means of comparing the methods. Therefore the 
output is summarized in Tables 5 through 7. They contain the output from the first five 
estimation methods only; the output from the vector test cases is presented later. Sam- 
ple output for the six methods is contained in Appendix C. 

The results summaries list the test case number, the cell inventory threshold. To, and 
the actual number of cells used, K. The level of expansion required to achieve these 
parameter levels is also listed. For example, test case one had a cell inventory threshold 
of 30.0, and 24 aggregated cells were obtained by expanding the small, large and major 
MOS group by YCSs four and five. The value for K is often different from the threshold 
number of cells, K^, listed in Table 3. When K is less than K^, maximum expansion oc- 
curred and the threshold was unattainable. When K is greater than K^, the expansion 
was the least amount possible to remain above the thresholds. From these test cases 
we can see that it is difficult to meet the threshold number of cells exactly. 

The results summaries then list the weighted average MOEs for each of the five es- 
timation methods. The first row within each test case contains the MAD values, the 
second row the chi square values and the third row the .VISE values. The maximum 
desired chi square value as given by Equation (61) is listed in parentheses, e.g., for test 
case one this value is 37.9. This affords easier comparison of the chi square values for 
the different methods. The values of MSE for both original scale methods are blank 
because MSE is a transformed scale MOE only. 

Before discussing the results in general, some additional comments about specific 
test cases are necessary. Test case four could not be executed by the Efron-.Morris 
method. This is because one of the aggregated cells contains zero attrition for all ten 
years of data. As a result, the iterative algorithm does not converge. 

Test case 12 was not possible because not even one aggregated cell meeting the cell 
inventory threshold was obtained with maximum expansion in major MOS group two. 
This extremely low inventoiy problem was generally true for all test cases involving the 
LtCol(FS)/26 YCS pair. Test cases 17, 21 and 24 obtained only three aggregated cells 
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with maximum expansion in major MOS group one. Test case 30 was also from major 
MOS group two, and obtained only one aggregated cell meeting the cell inventory 
threshold. As a result, test case 30 was changed to Major(FS)/I8 YCS so that results 
for these low thresholds could be obtained. 

As the thresholds became low (test cases 23-30) cells with inventory much larger 
than Tq were being obtained prior to any aggregation. This was especially true for the 
lLt/4 YCS and Capt/7 YCS pairs. To avoid masking the results of low inventory- 
thresholds by actually using large inventory cells, the service component/commissioning 
source parameters for test cases 26, 27 and 29 were changed. Rather than accepting all 
classifications of these parameters, only one classification for each was accepted. Thus 
test case 26 was executed with regular/USNA, test case 27 was executed with augmented 
regular/PLC-ground, and test case 29 was executed with regular/NROTC-scholarship. 

We now focus our attention on the results of these test cases with respect to the 
MOEs. The weighted average MAD figures vary little within each test case over the five 
methods. This suggests that the total deviation from the validation year empirical 
attrition rate was the same for all methods. However, this figure does not identify 
whether the deviations were overestimations or underestimations. The fraction of MAD 
from underage (not listed in the results summaries) was studied to gain more insight into 
this important consideration. For each of the 29 test cases (no results for test case 12), 
a weighted average fraction of .MAD from underage was computed for each of the five 
estimation methods (weighted by the number of cells just as the weighted averages for 
the MOEs). A weighted average of these 29 values was then computed. This overall 
weighted average indicates the tendency of the method to underestimate or 
overestimate-averages above 0.5 indicate a tendency to underestimate; averages below 
0.5 indicate a tendencv to overestimate. The author is unaware of anv information 
comparing the relative costs of underage and overage. Hence, as a default, we look for 
values of 0.5, which is a balance between overestimation and underestimation. The av- 
erages calculated were: TS 1 = 0.426; TS2 = 0.436; 051 = 0.481; 052 = 0.512; and 
E.M = 0.452. Although the MAD figures were generally the same for all methods, these 
averages indicate that the tendencies to overestimate or underestimate may not be the 
same. The original scale methods seem to have achieved more balance than the trans- 
formed scale methods. 

The chi square results were not entirely consistent across test cases nor across 
methods within a test case. These results are discussed first by comparisons between test 
cases; then by comparisons within test cases. 
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Of the first ten test cases, only three (Nos. 3, 5 and 8) had weighted average chi 
square values within the acceptable range. These three test cases expanded only into the 
large MOS group, whereas of the seven test cases which were unacceptable, all but one 
(No. 10) expanded into the major MOS group. Of the last 20 test cases, only four had 
chi square values outside the acceptable range (Nos. 13, 14, 15, and 18). Of these four, 
two expanded into the large MOS group, and two into the major MOS group. All the 
test cases with unacceptable values had either lLt/4 YCS or Capt/7 YCS pairs. They 
were fairly well spread across MOS groups. None of the test cases with lower thresholds 
(Nos. 19-30) had unacceptable chi square values. This suggests that lower thresholds, 
which result in less expansion, achieve more acceptable results with respect to this MOE. 

To investigate this claim further, different combinations of threshold levels for test 
cases seven and nine were used. The results are contained in Table 8. These results re- 
inforce the claim that lower thresholds are in fact better, since in both cases the chi 
square results improved as the thresholds, and therefore the level of expansion, were re- 
duced. When comparing these extra test cases, keep in mind that the chi square values 
for each set of threshold values should be compared to the acceptable range for that 
specific number of cells; comparisons across test cases with different values' of K are not 
valid. An important aspect of the argument for lower thresholds is that the thresholds 
must be considered jointly. For example, in test case seven, with a TJK pair of 10.0/6, 
the chi square values were nearly acceptable, whereas with 5.0/19 they were clearly un- 
acceptable. Thus we should be most aware of the value x Kq. 

We now turn our attention to comparing the chi square values within a test case. 
Several test cases had chi square values that varied significantly over the estimation 
methods (Nos. 4, 6 and 22). These three all had lLt/4 YCS pairs and were large MOS 
group four. In test cases four and six the chi square values for the transformed scale 
methods were not too much larger than the desired maximum; the values for the original 
scale methods were significantly larger than the desired maximum. In test case 22 the 
chi square values for the transformed scale methods were acceptable, however the values 
for the original scale methods were again significantly larger than the desired maximum. 

Several other test cases had varying chi square values, but to a lesser degree. In test 
case eight, only the transformed scale methods had acceptable chi square values, the 
original scale methods exceeded the desired maximum, although not by a significant 
margin. In test case 19, only the OS2 method exceeded the desired maximum. All 
methods for test cases 16, 25, and 27 were within the desired maximum, however the chi 
square values varied to a large degree over the five methods. Thus, using the chi square 
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MOE, it appears that the transformed scale methods were generally the same, and as a 
group outperformed the original scale methods. 

The .VISE was used only with the transformed scale methods, and thus no compar- 
ison with original scale methods can be made. The values for this MOE were generally 
equal between methods within a test case, and acceptable overall. 

The results for the vector method test cases are summarized in Table 9. The table 
lists the value Kx P (instead of K) because this is the number of estimates obtained and 
compared to the validation year empirical attrition rates with this method. 

Test case two with commissioning source as the vector component had to be modi- 
fied because only seven aggregated cells were obtained. As a result, K - P +2, and the 
vector method could not be conducted. Therefore, commissioning source three (Platoon 
Leader Class-ground) was deleted as a vector component and the test case run with only 
four commissioning sources (1, 2, 5 and 1 1). Test case 20 with service component as the 
vector component was infeasible. By starting with a low cell inventory threshold (5.0), 
when the cells were broken out into the vector components their inventory became ex- 
tremely low. As a result, when validating year five, two of the cells had zero inventory 
for all of the remaining nine years for service component three (reserve). Thus the value 
for XTB,j , given by Equation (45), becomes undefined and the method cannot be 
completed. 

The results of the vector method with respect to the MOEs is similar to the results 
of the previous five estimation methods. This method produced acceptable .MAD and 
MSE values, but its chi square values were fairly inconsistent. Test cases two and ten 
had unacceptable chi square values for both vector components; test case six had an 
unacceptable chi square value with service component as the vector component. Recall 
that these test cases also had unacceptable values with the vector component collapsed. 

No fair comparison with the first five estimation methods can be made using the 
summarized results. Obviously the .MAD and .MSE quantities will be larger since we are 
comparing three to five times as many estimates to empirical rates with the vector 
method. Thus a different evaluation technique must be used. 

The vector method is designed to take advantage of any correlation between the cells 
when broken out into vector components. To see if this is occurring we must look at 
the matrix S"‘ as given by Equation (54). In all of the vector test cases, this matrix was 
essentially diagonal, indicating little correlation between the vector components. In ad- 
dition, all of the eigenvalues, which become the elements of the diagonal matrix E, were 
less than {K — P — 2). Therefore, the eigenvalues were replaced by this quantity and the 
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Table 5. SUMMARY OF RESULTS (CASES 1-10) 



No. 


To 




Expansion 

Required 


MOE 


TSl 


TS2 


OSl 


OS2 


EM 








S:{4,5) 


MAD 


0.102 


0.101 


0.102 


0.102 


0.103 


1 


30.0 


24 


L:(4,5) 


X(37.9) 


98.79 


99.28 


100.99 


100.96 


101.23 








M:{4,5) 


MSE 


0.079 


0.078 






0.080 








S:{7) 


MAD 


0.091 


0.090 


0.092 


0.091 


0.091 


2 


30.0 


8 


L:(7) 


X(16.0) 


28.16 


27.68 


28.71 


28.68 


27.73 








M:{7) 


MSE 


0.077 


0.076 






0.076 








S:{l-3,6-19) 


MAD 


0.056 


0.056 


0.055 


0.055 


0.055 


3 


20.0 


31 


L:(6-8) 


X(46.7) 


37.49 


37.80 


39.58 


42.85 


38.38 










MSE 


0.062 


0.062 






0.062 








S:(l-5,8-19) 


MAD 


0.029 


0.029 


0.026 


0.024 




4 


20.0 


23 


L:{l-5,8-19) 


X(36.6) 


39.34 


47.36 


95.41 


139.14 










M;{ 1-5,8- 19) 


MSE 


0.035 


0.037 














S:{l-3,6-19) 


MAD 


0.048 


0.049 


0.047 


0.047 


0.048 


5 


20.0 


20 


L:(7) 


X(32.6) 


20.75 


20.83 


24.34 


24.20 


21.41 










MSE 


0.048 


0.049 






0.049 








S:(l-5.8-19) 


MAD 


0.029 


0.029 


0.027 


0.025 


0.029 


6 


20.0 


22 


L:( 1-5.8-19) 


X(35.3) 


40.42 


50.70 


90.93 


127.86 


65.20 








.M:{3-5) 


MSE 


0.037 


0.039 






0.040 








S:{4,5) 


MAD 


0.129 


0.130 


0.130 


0.130 


0.133 


7 


10.0 


32 


L:(4,5) 


X(48.0) 


111.67 


113.23 


112.64 


113.64 


114.28 








M:(4) 


MSE 


0.137 


0.137 






0.140 








S;{l-3,6-19) 


MAD 


0.039 


0.038 


0.039 


0.038 


0.039 


8 


10.0 


28 


L:{l-3,6-ll) 


X(43.0) 


34.37 


34.70 


45.16 


49.60 


35.54 










MSE 


0.039 


0.040 






0.040 








S:{7) 


MAD 


0.117 


0.115 


0.117 


0.117 


0.115 


9 


10.0 


14 


L;(7) 


X{24.6) 


38.85 


37.50 


39.59 


39.62 


40.26 








M:{7) 


MSE 


0.130 


0.127 






0.132 








S:(4.5) 


MAD 


0.128 


0.129 


0.129 


0.128 


0.127 


10 


10.0 


27 


L:(4,5) 


X(41.7) 


64.33 


65.92 


65.62 


66.61 


62.64 










MSE 


0.137 


0.137 






0.135 



TSl - Transformed scale, time dependent variance 
TS2 - Transformed scale, time independent variance 

051 - Original scale, time dependent variance 

052 - Original scale, time independent variance 
EM - Efron-Morris method 
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matrix E always had {K — P — 2) as its diagonal elements. Because these results indi- 
cated no worthwhile improvements over the first five methods, no further testing of the 
vector method was conducted. 

Finally, the performance of the heuristic aggregation algorithm listed in Appendix 
A was also evaluated. For each test case, the total inventory of cells below Tq was 
summed and this value divided by . The integer part of this number provides an upper 
bound on the number of aggregated cells that can be obtained by any aggregation 
technique. This upper bound was compared to the actual number of aggregated cells 
produced by the algorithm. The algorithm achieved the maximum in 71.4% (20 of 28) 
of the test cases. It achieved one less than the maximum in 21.4% (6 of 28) of the test 
cases, and two less than the maximum in 7.2% (2 of 28) of the test cases (only 28 test 
cases required aggregation: No. 12 was infeasible; No. 14 all cells were above ). This 
performance is acceptable for our application. 
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Table 6. SUMMARY OF RESULTS (CASES 11-20) 



No. 


T, 


K 


Expansion 

Required 


MOE 


TSl 


TS2 


OSl 


OS2 


EM 








S:(20-25) 


MAD 


0.109 


0.106 


0.108 


0.107 


0.107 


11 


10.0 


11 


L:(20-25) 


X(20.4) 


18.79 


19.75 


19.96 


20.53 


18.99 










MSE 


0.141 


0.142 






0.139 








S:(26) 


MAD 












12 


10.0 


0 


L:(26) 


X(O.O) 


















M;(26) 


MSE 


















S;(7) 


MAD 


0.117 


0.115 


0.117 


0.117 


0.115 


13 


10.0 


14 


L:{7) 


X(24.6) 


38.85 


37.50 


39.59 


39.62 


40.26 








M:(7) 


MSE 


0.130 


0.127 






0.132 








S:(4.5) 


MAD 


0.113 


0.110 


0.113 


0.113 


0.115 


14 


10.0 


10 


L:(4,5) 


X(18.9) 


57.80 


60.07 


57.71 


57.80 


59.36 








MSE 


0.102 


0.103 






0.104 








S:(4,5) 


MAD 


0.137 


0.140 


0.139 


0.139 


0.139 


15 


5.0 


32 


L:(4,5) 


X(48.0) 


69.96 


70.72 


70.27 


72.01 


68.15 








MSE 


0.156 


0.158 






0.155 








S:(20-25) 


MAD 


0.122 


0.122 


0.121 


0.121 


0.121 


16 


5.0 


32 


L:(20-25) 


X(48.0) 


39.41 


40.38 


44.35 


47.84 


42.27 








M:(20-22) 


MSE 


0.146 


0.146 






0.148 








S:(26) 


MAD 


0.166 


0.166 


0.168 


0.168 


0.168 


17 


5.0 


3 


L:(26) 


X(7.9) 


3.66 


3.57 


3.88 


3.78 


4.02 








M:(26) 


MSE 


0.0.189 


0.188 






0.202 








S:(6.7) 


MAD 


0.141 


0.141 


0.142 


0.141 


0.141 


18 


5.0 


29 


L:(6.7) 


X(44.2) 


79.04 


78.23 


78.61 


81.29 


78.84 








M;(7) 


MSE 


0.176 


0.176 






0.177 








S:(2,3,6-10) 


MAD 


0.085 


0.082 


0.080 


0.084 


0.080 


19 


5.0 


24 


X(37.9) 


28.95 


29.28 


30.91 


43.36 


35.44 










MSE 


0.122 


0.123 






0.132 








S;(20-25) 


MAD 


0.111 


0.113 


0.113 


0.106 


0.110 


20 


5.0 


19 


L;(20-25) 


X(31.3) 


17.29 


17.29 


18.92 


23.37 


18.28 








M:(20) 


MSE 


0.113 


0.114 






0.114 



TSl - Transformed scale, time dependent variance 
TS2 - Transformed scale, time independent variance 

051 - Original scale, time dependent variance 

052 - Original scale, time independent variance 
EM - Efron-Morris method 
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Table 7. SUMMARY OF RESULTS (CASES 21-30) 



No. 


T, 


K 


Expansion 

Required 


MOE 


TSl 


TS2 


OSl 


OS2 


EM 








S;(26) 


MAD 


0.167 


0.166 


0.168 


0.168 


0.168 


21 


5.0 


3 


L:(26) 


X(7.9) 


3.66 


3.57 


3.88 


3.78 


4.02 








M:(26) 


MSE 


0.189 


0.188 






0.202 








S:(l-5,8-19) 


MAD 


0.032 


0.042 


0.024 


0.022 


0.032 


22 


5.0 


22 


L:{3,4) 


X(35.3) 


23.56 


15.03 


63.10 


92.85 


22.56 








MSE 


0.059 


0.041 






0.056 








S:(20-25) 


MAD 


0.121 


0.121 


0.119 


0.121 


0.120 


23 


5.0 


9 


L:(20) 


X(17.5) 


8.12 


8.19 


8.65 


9.25 


8.17 








MSE 


0.141 


0.141 






0.137 








S:(26) 


MAD 


0.166 


0.166 


0.168 


0.168 


0.168 


24 


5.0 


3 


L:(26) 


X(7.9) 


3.66 


3.57 


3.88 


3.78 


4.02 








,M;(26) 


MSE 


0.189 


0.188 






0.202 








S;{ 1-6,8- 19) 


MAD 


0.050 


0.051 


0.020 


0.020 


0.048 


25 


5.0 


10 


L:(l-6) 


X(18.9) 


5.09 


4.78 


15.35 


16.54 


4.85 








MSE 


0.044 


0.037 






0.039 








S;(1-3,6-19) 


MAD 


0.138 


0.138 


0.138 


0.135 


0.135 


26 


5.0 


11 


L:(6,7) 


X(20.4) 


13.61 


13.65 


14.09 


16.10 


14.28 








MSE 


0.184 


0.186 






0.189 








S:{6-8) 


MAD 


0.070 


0.069 


0.055 


0.050 


0.068 


27 


5.0 


6 


X(12.9) 


4.64 


4.55 


6.88 


9.41 


4.60 










MSE 


0.075 


0.074 






0.074 








S:(20,21) 


MAD 


0.114 


0.116 


0.118 


0.113 


0.116 


28 


5.0 


7 




X(14.5) 


7.59 


7.67 


7.99 


8.74 


8.37 










MSE 


0.140 


0.143 






0.147 








S:(4,5) 


MAD 


0.170 


0.169 


0.168 


0.167 


0.169 


29 


5.0 


5 


X(11.3) 


6.36 


6.23 


7.64 


8.25 


8.06 










MSE 


0.211 


0.210 






0.281 








S:(l-6.8-19) 


MAD 


0.112 


0.110 


0.110 


0.104 


0.108 


30 


5.0 


6 


L:(17-19) 


X(12.9) 


6.15 


6.40 


7.31 


8.73 


6.42 










MSE 


0.136 


0.136 






0.135 



TSl - Transformed scale, time dependent variance 
TS2 - Transformed scale, time independent variance 

051 - Original scale, time dependent variance 

052 - Original scale, time independent variance 
EM - Efron- Morris method 
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Table 8. SUMMARY OF RESULTS (CASES 7 AND 9 EXPANDED) 



No. 


T, 


K 


Expansion 

Required 


MOE 


TSl 


TS2 


OSl 


OS2 


EM 


7-1 


10.0 


32 


S:(4,5) 

L:(4,5) 

M:(4) 


X(48.0) 


111.67 


113.23 


112.64 


113.64 


114.28 


7-2 


10.0 


17 


S:(4,5) 

L:(4) 


X(28.7) 


45.31 


45.99 


45.90 


46.58 


43.68 


7-3 


10.0 


6 


S:(4,5) 


X(12.9) 


13.08 


12.98 


13.41 


13^53 


13.03 


7-4 


5.0 


19 


S:(4,5) 

L:(4) 


X(31.3) 


49.28 


49.53 


49.23 


50.84 


46.91 


7-5 


5.0 


7 


S:(4,5) 


X(18.9) 


14.42 


14.27 


14.82 


14.93 


14.18 


7-6 


5.0 


4 


S:(4) 


X(9.7) 


9.20 


8.88 


9.22 


9.20 


9.49 


9-1 


10.0 


14 


S:(7) 

L:(7) 

M:(7) 


X(24.6) 


38.85 


37.50 


39.59 


39.62 


40.26 


9-2 


10.0 


4 


S:(7) 


X(9.7) 


10.35 


10.14 


10.58 


10.57 


10.11 


9-3 


5.0 


17 


S:(7) 

L:(7) 

M;(7) 


X(28.7) 


44.95 


42.87 


44.40 


46.03 


46.15 


9-4 


5.0 


5 


S;(7) 


X(11.3) 


11.91 


11.69 


12.16 


12.15 


11.64 



TSl - Transformed scale, time dependent variance 
TS2 - Transformed scale, time independent variance 

051 - Original scale, time dependent variance 

052 - Original scale, time independent variance 
EM - Efron-Morris method 
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Table 9. SUMMARY OF RESULTS (VECTOR METHOD) 



No. 


T, 


KP 


Vector 

Component 


MOE 


Vector 

Method 










MAD 


0.149 


2 


30.0 


24 


SC 


X(37.9) 


61.03 










MSE 


0.229 










MAD 


0.222 


2 


30.0 


28 


CS 


X(43.0) 


115.34 










MSE 


0.398 










MAD 


0.158 


3 


20.0 


93 


SC 


X(120.3) 


110.96 










MSE 


0.219 










MAD 


0.157 


3 


20.0 


155 


CS 


X( 190.2) 


132.37 










MSE 


0.201 










MAD 


0.062 


6 


20.0 


66 


SC 


X(89.0) 


109.19 










MSE 


0.074 










MAD 


0.083 


6 


20.0 


100 


CS 


X( 128.3) 


58.58 










MSE 


0.081 










MAD 


0.212 


10 


10.0 


81 


SC 


X( 106.5) 


251.63 










MSE 


0.368 










MAD 


0.321 


10 


10.0 


120 


CS 


X(151.0) 


707.85 










MSE 


0.625 










MAD 


0.193 


11 


10.0 


42 


SC 


X(60.3) 


44.86 










MSE 


0.214 










MAD 


0.222 


11 


10.0 


60 


CS 


X(81.9) 


53.55 










MSE 


0.225 










MAD 




20 


5.0 


72 


SC 


X(O.O) 

MSE 












MAD 


0.254 


20 


5.0 


100 


CS 


X( 128.3) 


55.92 










MSE 


0.245 



SC - service component 
CS - commissioning source 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The results indicate that the desired stability in estimating attrition rates for low in- 
ventor}’ cells has been achieved with the aggregation and estimation methods presented 
in this study. The use of “shrinkage” methods applied to well selected groups of cells 
allows for the achievement of quality estimates of attrition in the face of low inventor}’ 
numbers for the individual cells. 

None of the six estimation methods stood out as a clear favorite. The vector 
method did not provide any additional benefits using service component or commis- 
sioning source as vector components. Since it is a more complicated method and has 
the potential to become unsolvable with zero inventory vector components, it appears 
to be the least favorite. Perhaps more success would be obtained with alternative clas- 
sifications for the vector component. 

The Efron-Morris method also involves more computational effort than the first 
four empirical Bayes methods. Its performance was very much similar to the trans- 
formed scale, time dependent variance method since the only difference between them 
is the iterative algorithm used to determine the amount of shrinkage. The Efron-Morris 
method has the potential to become unsolvable when a cell has zero attrition for eveiy 
year-a distinct possibility when dealing with low inventory cells. This suggests that it 
is the least favorite of the first five methods. 

Of the remaining four methods, there seems to be only small difference between the 
time dependent variance and the time independent variance methods on the same scale. 
In test cases where the chi square values were marginal or unacceptable, the time de- 
pendent variance methods were usually better. In these same test cases, the transformed 
scale methods performed better than the original scale methods. Therefore, if one 
method was to be singled out as best, it would be the first method; transformed scale, 
time dependent variance. 

The tendency to overestimate or underestimate as shown by the weighted average 
fraction of .MAD from underage may also be a consideration when selecting a method. 
An analysis of this type must weigh the costs of overestimating versus the costs of 
underestimating, which generally are not the same. This type of an analysis is beyond 
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the scope of this study. In addition, further testing of the methods would be required 
to gain a more accurate estimate of this tendency. 

The threshold levels also seem to strongly influence the performance of the esti- 
mation methods. It appears that expansion past the large MOS group begins to detract 
from homogeneous attrition behavior. While further study would be required to identify 
optimal threshold levels, it is apparent that both thresholds should not exceed 20.0, and 
the value of Ti, x should not exceed 100. 

A method for dealing with cells whose inventory is much greater than must be 
developed. In some test cases, cells with inventory three or more times as large as Tq 
were obtained and used in the estimation process. This did not seem to afiect the results, 
as they were present in almost all test cases. These cells could be disaggregated into 
multiple cells with inventory closer to the threshold, although the effect of this has not 
been determined. 

B. RECOMMENDATIONS 

The proposed aggregation method should be implemented as a method of identifying 
additional cells to be used in the attrition rate estimation process. This method provides 
greater homogeneity of attrition behavior among cells over the current method. 

The empirical Bayes estimation methods developed in this study are recommended 
for use in estimating the attrition rates for low personnel inventory cells. 

Overall, the empirical Bayes estimation methods when combined with the proposed 
aggregation method have achieved the stability in attrition rate estimation that is re- 
quired to provide a foundation for manpower planning. 
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APPENDIX A. AGGREGATION ALGORITHMS 



A. HEURISTIC ALGORITHM 

The heuristic algorithm for aggregating cells is as follows: 

1. Given a set of cells, S, and the (time average) inventory of each cell, INV^, partition 
S into two subsets as follows: 

5, = { c : c e S ; INV, ^ 7„ } 

S, = ( c : c 6 S ; INV, < T , } 

2. Put the cells in Sj into the set of aggregated cells, K. 

3. Order the cells in Sj according to size of their inventory: 

INV, < INV, < . . . < INV„ n = 1 S, 1 

4. Start with c„, the cell in with the largest inventory. Find the smallest cell in S^, 
c , that when united with c„ the resulting total inventory will meet or exceed . 
Combine its data with c„, put c„ into K, and remove c from Sj (the modified set 
will now be referred to as SJ). Repeat the procedure with c„_j , and so forth. 

5. If no cell in SJ when combined with the current largest cell, c„_, , exceeds Tp, use 
the next largest cell, . and remove from SJ . This will create an aggregated 
cell that is still below threshold. Return to the procedure in Step 4 of trx'ing to find 
c' . If no such cell is contained in SJ, use and so forth. 

6. Continue this procedure until the sum of all the cells remaining in is less than 
Tp. These cells are sequentially added to the aggregated cells in K in Step 7. 

7. Add the largest cell in to the smallest cell in K, and update its average inventory. 
Add the next largest cell in to the current smallest cell in K, and update the in- 
ventor}'. Continue until all cells in SJ have been used. 

We now have | K \ aggregated cells which exceed the threshold, Tp , to use in the 
attrition rate estimation procedure. 

B. INTEGER LINEAR PROGRAM 

The formulation as an integer linear program is as follow's; 

Index Use 

c cell (before aggregation) 

a aggregated cell 

Given Data 

7.VF, average inventory of cell c 
Tp threshold cell inventory 
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Binarv' Variables 



1 indicates use cell c in aggregated cell a 
Z„ 1 indicates use cell a for aggregation 

Formulation 
MAX I Z, 

a 

subject to 

Z ^ 1 V c 
t INV, . > T,- Z, Wa 

i;„,z, e {0,1} 



(each cell used at most once) 
(aggregated cell must have size > Tj ) 
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APPENDIX B. COMPUTER PROGRAMS 






A. GENERAL 

A computer program written in FORTRAN is used to conduct the cross validation 
using the methods developed in this thesis. Although the program consists of 33 sub- 
routines, 6 function subroutines, and over 2000 lines of code, it can be easily summarized 
by breaking it into the two areas of the thesis: cell aggregation and estimation methods. 

The main program and aggregation subroutines (listed in paragraph B) read the in- 
put parameters and execute the expansion and aggregation methods discussed in Chap- 
ter II. An existing program written by Luis Uribe, an independent contractor under the 
direction of Professor Read, underwent extensive modification to fulfill these require- 
ments. The input parameters- To, , MOS, YCS, grade, service component(s), and 
commissioning source(s)-are read by Subroutine GETPAR either in the interactive 
mode via the terminal or by using MC87 EXEC (listed in paragraph E). Uribe uses an 
innovative method to estimate the amount of expansion required to meet the threshold 
parameters. This approach precludes the requirement to read the data base after each 
step in the expansion process which would be extremely computer time intensive. In- 
ventorv’ information is extracted from the data base and stored in a separate data file for 
each pay grade (a sample data file and the program used to create it are listed in para- 
graph F). The data file is accessed via the user's A-disk which is significantly faster than 
accessing the data base through MVS. Subroutine READET reads the appropriate data 
file for the specified grade and constructs a table of cells for those records that are in the 
same major MOS group as the user defined MOS, and meet the service component and 
commissioning source parameters. All YCSs are accepted, since the extent of expansion 
is not yet determined. Function NCEVAL screens this table using the current level of 
expansion and estimates the number of aggregated cells with average inventor>’ greater 
than or equal to To that will be obtained. If this number is less than K^, Subroutine 
EXPAND begins the expansion stages as described in paragraph Il.C.. After each in- 
crement of expansion, NCEVAL screens the table and estimates the number of aggre- 
gated cells that will be obtained. This loop through EXPAND and NCEVAL continues 
until the estimated number of aggregated cells meets the threshold, A^o- The estimated 
number of cells and the level of expansion are then displayed on the terminal screen. 
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The user may elect to go forth and read the data base to determine the actual number 
of aggregated cells obtained, or may elect to change the level of expansion. 

The level of expansion is changed through the variable AGGPCT. This variable 
estimates the effectiveness of the heuristic aggregation method listed in Appendix A. 
To estimate the number of aggregated cells that will be obtained, NCEVAL compares 
the cells which meet the expansion criteria to the minimum inventory threshold, T„. 
Those that are greater than will obviously produce one aggregated cell. The inventory 
of those that are less than is summed. The estimated number of cells is then the total 
of the number of cells greater than and AGGPCT times the sum of the cell inventory 
below To divided by Tq. The variable AGGPCT is initially set at 0.9, but can be inter- 
actively changed via the terminal. By increasing the value of AGGPCT we can decrease 
the level of expansion; by decreasing the value of AGGPCT we can increase the level 
of expansion. 

Once we decide to go forth and read the actual data base. Subroutine READER 
extracts records meeting the expansion criteria developed using EXPAND and 
NCEVAL and pools them into cells. Subroutine AGGREG aggregates these cells to 
meet the average inventory threshold, . The actual number of aggregated cells ob- 
tained is then compared to the threshold number of cells, Xo. Again, the user has the 
option of changing the level of expansion to obtain more or fewer cells, or continuing 
on to the estimation process. 

The first five estimation methods are called by SUBROUTINE MC87BZ (listed in 
paragraph C). The estimation methods are contained in separate subroutines: EBTSl, 
EBTS2, EBOSl, EBOS2 and EMTS (EB-empirical Bayes; EM-Efron-Morris; 
TS-transformed scale; OS-original scale; 1-time dependent variance; 2-time independent 
variance). The iterations required by the first four methods are conducted in Subroutine 
EBITER; the Efron-Morris iterations are conducted in Subroutine EMITER. The 
MOEs are then computed by Subroutines MSE and OSMOE. 

If the vector method is to be used. Subroutine BKDOWN then breaks the cells out 
by their vector components (a vector of length three for service component; a vector of 
length five for commissioning source). The vector estimation method is contained in 
Subroutine MC87V (listed in paragraph D). Since all of its computations are unique, 
this subroutine is self-contained with the exception of the transformation formula, which 
is contained in Function FTTV. 
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B. MAIN PROGRAM AND AGGREGATION SUBROUTINES 



--- PROGRAM TO CONDUCT AGGREGATION AND ESTIMATION METHODS 



C 



C 

C 



C 



C 



C 



--- PARAMETER MXY MUST BE UPDATED TO REFLECT EXACT NO. YEARS OF DATA 
--- PARAMETER MXP IS MAX LENGTH OF 3RD DIMENSION P-VECTOR 
--- PARAMETER MXK IS MAX NUMBER OF AGGREGATED CELLS (MAX NO) 
PARAMETER (MXX=600, MXY=10, MXP=6, MXK=50) 

PARAMETER (NMS=81, NG=14, NLG=6, NMG=4, NYB=4, NYE=18, NYEG=4) 

INTEGER ST1,ST2,LYR 
INTEGER SYCS(31), NYCS 

INTEGER SYCSG(31) ,SYCSL(31) ,SYCSM( 31) ,NYCSG,NYCSL,NYCSM 

INTEGER SMOS(30), NMOS 

INTEGER SVCMP(S), NSC 

INTEGER SCSRC(16),NCSR 

INTEGER SGRD 

INTEGER*2 M0SGR(2 ,NMS) , YCSB(NYE ,NYB ,NYEG) , VYC(NYE) 

INTEGER*2 LGRP(NG), MGRP(NLG) 

REAL INV(MXX,MXY),Y(MXX,MXY), SINV(MXX,MXY) ,SY(MXX,MXY) 

INTEGER DATA(MXY) 

--- ARRAYS FOR MC87BZ 

REAL XTB(MXX) .VXTB(MXX) ,XEB(MXX) ,A(MXX) 

--- ARRAYS FOR MC87V 

REAL XTBJI(MXP,MXK), DELTA(MXP,MXK) , X(MXP,MXK) 

REAL XVYR(MXP.MXK), VYRINV(MXP.MXK) , VYRY(MXP ,MXK) 

REAL BSTAR(MXP,MXP), S(MXP,MXP), GAMMA(MXP,MXP) 

REAL XBBJ(MXP), EVAL(MXP) 



INTEGER*2 PTRTBL(MXX, 2), INDX(MXX), MKG(MXX), RETTBL(MXX, 3) 
INTEGER*2 PTBL(MXX, 3), BKTBL(MXX,3) 

REAL AVINV(MXX), RETINV(MXX) 

DATA MKG/MXX*0/ 

--- ASSIGN MOS TO SMALL, LARGE AND MAJOR MOS GROUP 



DATA MOSGR/013,1, 


020,2 


, 027,2, 


, 038 


,2, 


039 


,2, 










005,3, 


007,3, 


049,3, 


052 


,3, 
















* 


074,4, 


079,4, 


085,4, 


101 


,4, 
















* 


016,5, 


060,5, 


064,5, 


076 


,5, 


111, 


5 , 


116, 


5 , 










132,6, 


134,6, 


135,6, 


139 


,6, 
















* 


143,7, 


147,7, 


150,7, 


153 


,7, 


154, 


7, 


155, 


7, 170,7; 


> 






149,8, 


151,8, 






















* 


160,9, 


161,9, 


164,9, 


166 


,9, 


167, 


9, 


168, 


9, 178,9; 






* 


173,10 


, 174,10 


, 175, 


10, 


176, 


,10, 


177 


,10, 


179, 


10, 


144, 


,10 




145,10 


, 165,10 


> 




















* 


001,11 


, 006,11 


, 012, 


11, 


015, 


,11, 


019 


,11, 


026, 


11, 


037, 


,11 


* 


048,11 


, 051,11 


, 059, 


11, 


070, 


,11, 


075 


,11, 


078, 


11, 


084, 


,11 


* 


087,11 


, 100,11 


, no. 


11, 


115, 


,11, 


131 


,11, 


138, 


11, 


217, 


,11 


* 


172,12 


, 187,12 


, 188, 


12, 


189, 


,12, 














* 


142,13 


, 146,13 


, 148, 


13, 


152, 


,13, 


156 


,13, 


163, 


13, 


169, 


,13 


* 


088,14 


/ 























DATA LGRP/1,1,4*2,3,3,4,4,5,5,5,6/ 

DATA MGRP/1,1,2,2,3,4/ 

--- CREATE YCS EXPANSION BOUNDS 

DATA YCSB/1,2,3,4,5,6, 8,9,10,11,12,13,14,15,16,17,18,19, 

* 7,17*0, 20,21,22,23,24,25,12*0, 26,17*0, 

* 1,2, 3, 4, 5, 8,9,10,11,12,13,14,15,16,17,18,19,1*0, 



MC800010 

MC800020 

MC800030 

MC800040 

MC800050 

MC800060 

MC800070 

MC800080 

MC800090 

MC800100 

MC800110 

MC800120 

MC800130 

MC800140 

MC800150 

MC800160 

MC800170 

MC800180 

MC800190 

MC800200 

MC800210 

MC800220 

MC800230 

MC800240 

MC800250 

MC800260 

MC800270 

MC800280 

MC800290 

MC800300 

MC800310 

MC800320 

MC800330 

MC800340 

MC800350 

MC800360 

MC800370 

MC800380 

MC800390 

MC800400 

MC800410 

MC800420 

MC800430 

MC800440 

MC800450 

MC800460 

MC800470 

MC800480 

MC800490 

MC800500 

MC800510 

MC800520 

MC800530 

MC800540 
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* 

* 

* 

* 

* 



6,7,16*0, 20,21,22,23,24,25,12*0, 26,17*0, 

1,2, 3, 4, 5, 6, 8,9,10,11,12,13,14,15,16,17,18,19, 
7,17*0, 20,21,22,23,24,25,12*0, 26,17*0, 

1,2,3, 6,7,8,9,10,11,12,13,14,15,16,17,18,19,1*0, 
4,5,16*0, 20,21,22,23,24,25,12*0, 26,17*0 / 



C --- 



2 

1 

C --- 



C - — 



C --- 



C 



INITIALIZE INVENTORY AND ATTRITION ARRAYS 
DO 1 1=1, MXX 
DO 2 J=1,MXY 
SINV(I,J)=0 
SY(I,J)=0 
INV(I,J)=0 
Y(I,J)=0 
CONTINUE 
CONTINUE 

DEFINE FILE FOR OUTPUT 

CALL EXCMSCFILEDEF 11 DISK MC87 OUTPUT A') 

FIRST/ LAST YEAR OF DATA ON TAPE. UPDATE WHEN NECESSARY 

ST1=77 

NYR=MXY 

LYR=ST1+NYR-1 

INITIAL VALUE FOR AGGREGATION ESTIMATION PERCENTAGE 

AGGPCT=0. 9 

ICYCLE=1 

GET INPUT PARAMETERS 

CALL GETPARC AIMIN,NO ,NMOS , SMOS ,NYCS , SYCS ,SGRD , 

^ NSC,SVCMP, NCSR,SCSRC, IGR,MOSGR,NMS , ISFLAG) 

MAJOR GROUP IS MG, LARGE GROUP LG, GROUP IGR, YCS BLOCK lY 
LG=LGRP(IGR) 



MC800 

MC800 

MC800; 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800 

MC800f 

MC800f 



C --- 



C --- 

10 



MC800f 

MG=MGRP(LG) MC8006 

WRITE(6,*) ' ’ MC800£ 

WRITE(6,*) ' GR,LG,MG=’ ,IGR,LG,MG MC800£ 

WRITE(6,*) ' ’ MC800S 

READ EVALUATION TABLE. SELECT ONLY RECS PASSING SELECTION CRITERIAMC8008 

CALL READET(RETTBL,RETINV,MXX,NRET, SGRD, NSC,SVCMP, NCSR,SCSRC, MC8008 
* MG,LGRP,MGRP, MOSGR,NMS) 

RC=0 
IGX=IGR 
LGX=0 
MGX=0 
NYCSG=1 

SYCSG(1)=SYCS(1) 

NYCSL=1 

SYCSL(1)=SYCS(1) 

NYCSM=1 

SYCSM(1)=SYCS(1) 

NCTOT=0 

NC=NCEVAL( AIMIN , IGX , LGX , MGX , NYCSG , SYCSG , RETTBL , RETINV , NRET , MXX , 

LGRP , MGRP , NMS , AGGPCT , I GR , LG ) 

DO WHILE NCTOT<NO & RC=0 (EXPAND AS LONG AS NO NOT MET) 

IF(NC .GE. NO) THEN 

WRITE(6,*) '$GG EVAL NC, SYCSG=' ,NC, ( SYCSG( II ), 11=1 , NYCSG) 

GO TO 60 
END IF 

IF( NYCSG. EQ. 1) THEN 

CALL GETVYCC SYCS ( 1 ) , LG , YCSB , NYE , NYB , NYEG , VYC ) 

WRITE(6,*) '== VYC=' ,(VYC(I),I=1,NYE) 



MC8008 

MC8009 

MC8009 

MC8009 

MC8009 

MC8009 

MC8009 

MC8009 

MC8009 

MC8009 

MC8009 

MC8010 

MC8010 

MC8010J 

MC8010;. 

MC8010-' 

MC8010."' 

MC8010( 

MCSOlO: 

MC8010{ 

MC8010' 

MC801K 
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no on on 



ENDIF MC801110 

CALL EXPAND(NyCSG,SYCSG,VYC,NYE,IGX,LGX,MGX,LG,HG,RC) MC 801120 

IF(IGX .EQ. 0) GO TO 20 MC801130 

NC=NCEVAL( AIMIN , IGX , LGX , MGX , NYCSG , SYCSG , RETTBL , RETINV , NRET , MXX , MC80 1 140 

* LGRP,MGRP,NMS,AGGPCT,IGR,LG) MC801150 

GO TO 10 MC801160 

20 NCTOT=NC MC801170 

WRITE(6,*) '$$G EVAL NC,SYCSG=' ,NCTOT,(SYCSG( II), 11=1, NYCSG) MC801180 

MC801190 

--- EXPAND TO LARGE MOS GROUP MC801200 

WRITE(6,*) ' ' MC801210 

WRITE(6,*) '= EXPANDING BY LARGE GROUP: ' ,LGX MC801220 

NC=NCEVAL( AIMIN , IGX , LGX , MGX , NYCSL , SYCSL , RETTBL , RETINV , NRET , MXX , MC 80 1 2 30 

* LGRP,MGRP,NMS,AGGPCT,IGR,LG) MC801240 

30 IF((NCTOT+NC) . GE. NO) THEN MC801250 

WRITE(6,*)' $LL EVAL NC,SYCSL=' , (NCTOT+NC) , ( SYCSL( II) , II=1 ,NYCSL)MC801260 
GO TO 60 MC801270 

ENDIF MC801280 

IF(NYCSL. EQ. 1) CALL GETVYC( SYCS( 1) ,LG,YCSB ,NYE ,NYB ,NYEG,VYC) MC801290 

CALL EXPAND(NYCSL,SYCSL,VYC,NYE,IGX,LGX,MGX,LG,MG,RC) MC801300 

IF(LGX .EQ. 0) GO TO 40 MC801310 

NC=NCEVAL( AIMIN, IGX, LGX, MGX, NYCSL, SYCSL, RETTBL, RETINV, NRET, MXX, MC801320 

* LGRP,MGRP,NMS,AGGPCT,IGR,LG) MC801330 

GO TO 30 MC801340 

40 NCTOT=NCTOT+NC MC801350 

WRITE(6,*) '$$L EVAL NC ,SYCSL=' ,NCT0T,( SYCSL( II) , 11=1 , NYCSL) MC801360 

MC801370 

--- EXPAND TO MAJOR MOS GROUP MC801380 

WRITE(6,*) ' ' MC801390 

WRITE(6,*) '= EXPANDING BY MAJOR GROUP: ' ,MGX MC801400 

NC=NCEVAL( AIM IN , IGX , LGX , MGX , NYCSM , SYCSM , RETTBL , RETINV , NRET , MXX , MC80 14 10 

* LGRP,MGRP,NMS,AGGPCT,IGR,LG) MC801420 

50 IF( (NCTOT+NC) . GE. NO .OR. RC . NE. 0) THEN MC801430 

WRITE(6,*)'$MM EVAL NC,SYCSM=' , (NC+NCTOT) , ( SYCSM( II) , 11=1 ,NYCSM)MC801440 

MC801450 
MC801460 
MC801470 
MC801480 
MC801490 
MC801500 
MC801510 
MC801520 
MC801530 
MC801540 



GO TO 60 
ENDIF 

IF( NYCSM. EQ. 1) CALL GETVYC( SYCS( 1) ,LG,YCSB,NYE,NYB ,NYEG, VYC) 
CALL EXPAND ( NYC SM , S YC SM , VYC , NYE , I GX , LGX , MGX , LG , MG , RC ) 

NC=NCEVAL( A IM IN , I GX , LGX , MGX , NYC SM , S YC SM , RETTB L , RETI NV , NRET , MXX , 
* LGRP,MGRP,NMS,AGGPCT,IGR,LG) 

GO TO 50 

--- EXPANSION FINISHED 
60 IF(RC .NE. 0) THEN 



WRITE(5, *)'*** REQUIRED NO MAY NOT BE MET: N0,NC=' ,N0,(NC+NCTOT)MC801550 



ENDIF 

C --- ALLOW USER TO CHANGE EXPANSION LEVEL 

WRITE(5,*) 'ESTIMATED NUMBER OF CELLS =' ,NC+NCTOT 
70 WRITE(5,*) 

WRITE (5,*) 'ENTER 1 TO CALL READER, 0 TO CHANGE EXPANSION' 
READ(5,*) NPICKl 
IF(NPICK1 .EQ. 1) THEN 
GO TO 80 
ELSE 

WRITE(5,*) 'AGGPCT IS CURRENTLY =' , AGGPCT 
WRITE(5,*) 'ENTER NEW VALUE FOR AGGPCT' 



MC801560 
MC801570 
MC801580 
MC801590 
MC801600 
MC801610 
MC 80 1620 
MC801630 
MC801640 
MC801650 
MCS01660 
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o o o o o o o 



80 



RE AD (5,*) AGGPCT 
GO TO 5 
ENDIF 

WRITE(5,*) 'CALLING READER' 



90 



- USER ELECTS TO READ THE DATA BASE - DETERMINE MOS EXPANSION LEVEL 
CALL GETMOS ( SMOS , NMOS , MGX , LGX , MG , LG , I GR , MOSGR , LGRP , MGRP , 

* NMS,NG,NLG) 

■ READ THE DATA BASE AND CREATE THE CELLS 

CALL READERC DATA , INV , Y , MXX , NMOS , NYCSG , NYCSL , NYC SM , NSC , NCSR , NYR , 

* SMOS , S YCSG , SYCSL , SYCSM , SGRD , SVCMP , SCSRC , NRC , PTRTBL , LGX , MGX , IGR , 

* LG , MGRP , LGRP , MOSGR , NMS , NG , NLG , ICYCLE , NPT , PTBL , I SFLAG , S INV , SY ) 

- PERFORM CELL AGGREGATION TO MEET INVENTORY THRESHOLD 
CALL AGGREGC INV , Y , MXX , NYR , SMOS , SYCSG , 

* NRC, NRCOLD, PTRTBL, INDX,AVINV,AIMIN,MKG) 

■ ALLOW USER TO CHANGE EXPANSION LEVEL 
WRITE(5,*) 'NUMBER OF CELLS =' ,NRC 
WRITE(5,*) 

WRITE(5,*) 'enter 1 TO CONTINUE, 0 TO CHANGE EXPANSION' 

READ(5,*) NPICK2 
IF(NPICK2 .EQ. 1) THEN 
GO TO 100 
ELSE 

WRITE(5,*) 'AGGPCT IS CURRENTLY =' , AGGPCT 
WRITE(5,*) 'ENTER NEW VALUE FOR AGGPCT' 

READ(5,*) AGGPCT 
ICYCLE=ICYCLE+1 
GO TO 5 
ENDIF 



100 



C -- 
C -- 



* USER ELECTS TO CONDUCT ESTIMATION 
CONTINUE 

WRITE( 11,201) 'EXPANSION INFORMATION: ' 

WRITE (11, 203) 'ACTUAL NO. OF CELLS USED= ' ,NRC 
WRITE( 11,202) 'MOS GROUP #',IGR,' YCS"S USED=' , 

* (SYCSG(I),I=1,NYCSG) 

IF(LGX .GT. 0) THEN 

WRITE( 11,204) 'large MOS GROUP ?'/',LG,' YCS"S USED=' , 

* (SYCSL(I),I=1,NYCSL) 

ELSE IF(MGX . GT. 0) THEN 

WRITE( 11, 204) 'LARGE MOS GROUP s'/', LG,' YCS"S USED=' , 

* (SYCSL(I),I=1,NYCSL) 

WRITE( 11,204) 'MAJOR MOS GROUP #',MG,' YCS"S USED=' , 

* (SYCSM(I),I=1,NYCSM) 

ENDIF 

■ PERFORM ALL BUT VECTOR ESTIMATION METHODS IN MC87BZ 
CALL MC87BZ(INV,Y,NRC,NYR,XTB,VXTB,XEB,A,MXX,MXY) 

■ VECTOR METHOD- -BREAK CELLS INTO VECTOR, CONDUCT ESTIMATION 
IFdSFLAG .GT. 0) THEN 

CALL BKDOWNC PTBL , NPT , PTRTBL , NRCOLD , INDX , MKG , MXX , MXY , 

* SINV,SY,INV,Y,BKTBL,NBK) 

CALL MC8 7 V( INV , Y , MXX , NYR , NRC , XTB JI , DELTA , X , XVYR , VYRINV , VYRY , 

* BSTAR , S , GAMMA , XBB J , EVAL , MXP , MXK , BKTBL, NBK , NSC , NCSR , I SFLAG) 
ENDIF 



201 F0RMAT(/1X,A) 



MC801,? 
MC80HS 
MC801i? 
MCSOllj 
MC801I 
MC801J 
MC801'.! 
MC80li 
MC80i: 
MC80i: 
MC80i: 
MC80i: 
MC80i: 
MC801J 
MC801J 
MC8015 
MC801£ 
MC801? 
MC801{ 
MC8018 
MC8018 
MC8018 
MC8018 
MC801 
MC801S 
MC801S 
MC8019 
MC8019 
MC8019 
MC8019! 
MC8019; 
MC8019 
MC8019 
MC8020 
MC8020. 
MC8020. 
MC8020. 
MC 8 O 20 ; 
MC 8 O 20 : 
MC8020) 
MC8020) 
MC8020J 
MC8020!! 
MC80211 
MC8021. 
MC802l;f 
MC802i;i 
MC8021<I 
MC8021i-i 
MC8021('i 
MC802i:| 
MC8021J 
MC8021'! 
MC8022C' 
MC8022] 
MC80221 
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202 


FORMATC IX, A, 12, A/ IX, 18(13)) 


MC802230 


203 


FORMATC IX, A, 12) 


MC802240 


204 


FORMATC 1X,A,I1,A/1X, 18(13)) 


MC802250 




END 


MC802260 


C 




MC802270 


***********^v************************->v**********************vr***********VcMc 802280 


c 




MC802290 




SUBROUTINE EXPAND( NYCSX , S YCSX , VYC , NYE , I GX , LGX , MGX , LG , MG , RC ) 


MC802300 


c --- 


EXPAND YCS IF FEASIBLE, ELSE EXPAND MOS TO LG/MG 


MC802310 




INTEGER SYCSX(31), NYCSX 


MC802320 




INTEGER*2 VYC (NYE) 


MC802330 


c --- 


FIND POSITION OF ORIGINALLY REQUESTED SYCS(l) 


MC802340 




IY=0 


MC802350 




DO 10 1=1, NYE 


MC802360 




IF(SYCSXd) .EQ. VYC(I)) IY=I 


MC802370 


10 


CONTINUE 


MC802380 




IF(IY. EQ. 0) GO TO 30 


MC802390 


c --- 


FIND NEAREST NON-ZERO YCS TO USE FOR EXPANSION 


MC802400 




DO 20 1=1, NYE 


MC802410 




J=IY-I 


MC802420 




IF(J. GE. 1) THEN 


MC802430 




IF(VYC(J).GT. 0) GO TO 50 


MC802440 




ENDIF 


MC802450 




J=IY+I 


MC802460 




IF(J, LE.NYE) THEN 


MC802470 




IF(VYC(J).GT. 0) GO TO 50 


MC802480 




ENDIF 


MC802490 


20 


CONTINUE 


MC802500 


30 


CONTINUE 


MC802510 


C --- 


NO MORE YCS EXPANSION POSSIBLE. SEE IF MOS EXP. FEASIBLE 


MC802520 




IF(IGX. GT. 0) THEN 


MC802530 


C --- 


EXPAND FROM GROUPS TO LARGE GROUP LGX 


MC802540 




IGX=0 


MC802550 




LGX=LG 


MC802560 




ELSE IF(LGX. GT. 0) THEN 


MC802570 


C --- 


EXPAND FROM LARGE GROUP LGX TO MAJOR GROUP MGX 


MC802580 




LGX=0 


MC802590 




MGX=MG 


MC802600 




ELSE 


MC802610 




RC=1 


MC802620 




ENDIF 


MC802630 




RETURN 


MC802640 


C 




MC802650 


C --- 


EXPAND WITH YCS IN POSITION J & CLEAR VYC(J) 


MC802660 


50 


CONTINUE 


MC802670 




NYCSX=NYCSX+1 


MC802680 




SYCSX(NYCSX)=VYC(J) 


MC802690 




VYC( J)=0 


MC802700 




END 


MC802710 


C 




MC802720 


*Vc**^Wc**^V***********Vr***^V*******'>Wc**'>V******^V**^V**^V******^V'>Wc***********'>V'>V 802730 


C 




MC802740 




FUNCTION NCEVALC AIMIN , IGX , LGX , MGX , NYCSX , SYCSX , RETTBL , RETINV , 


MC802750 




NRET , MXX , LGRP , MGRP , NMS , AGGPCT , IGR , LG) 


MC802760 



C --- COMPUTE ESTIMATED NO. CELLS TO BE OBTAINED WITH CURRENT EXPANSION MC802770 
INTEGER SYCSX(31) ,NYCSX MC802780 
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INTEGER*2 LGRP( 14) ,MGRP(6) MC802'. 

INTEGER*2 RETTBL(MXX, 3) MC802( 

REAL RETINV(MXX) MC802( 

NCEVAL=0 ^ MC802( 

IF(IGX. EQ. 0 .AND. LGX. EQ. 0 .AND. MGX. EQ. 0) RETURN MC802i 

TAINV=0. 0 MC802f 

DO 100 I=1,NRET MC802f 

C --- SCREEN ON YCS MC802E 

DO 10 J=1,NYCSX MC802f 

IF(RETTBL(I,2) . EQ. SYCSX(J)) GO TO 15 MC802£ 

10 CONTINUE MC8028 

GO TO 100 MC8029' 

C --- SCREEN ON MOS BY GROUP, L. GRP OR MG DEPENDING ON IGX,LGX,MGX MC802S 

15 CONTINUE MC802S 

M0S=RETTBL(I,1) MC8029 

IGP=RETTBL(I,3) MC802S' 

LGP=LGRP(IGP) MC802S 

IF(MGX .GT. 0) THEN MC8029: 

IF(MGRP(LGP) .EQ. MGX) THEN MC8029 

IF(LGP .NE. LG) GO TO 80 MC8029 

ENDIF MC8029 

ELSE IF(LGX . GT. 0) THEN MC8030 

IF(LGP .EQ. LGX) THEN MC8030 

IF(IGP .NE. IGR) GO TO 80 MC8030 

ENDIF MC8030: 

ELSE MC8030; 

IF(IGP .EQ. IGX) GO TO 80 MC8030 

ENDIF MC8030. 

GO TO 100 MC8030; 

80 CONTINUE MC8030: 

C --- ACCEPTED MC8030: 

IF(RETINVd) .GE. AIMIN) THEN MC80311 

NCEVAL=NCEVAL+1 MC80311 

ELSE MC8031 

TAINV=TAINV+RETINV(I) MC8031, 

ENDIF MC8031; 

100 CONTINUE MC8031i 

C --- FINAL ESTIMATE IS NCEVAL MC8031l> 

IF(AIMIN. GT. 0) NCEVAL=NCEVAL + AGGPCT^dAINV/AIMIN MC80311 

END MC8031I 

C MC8031; 



Vf***Vf VfVfVfVfVf *Vc***Vf Vf ******** ******Vf**VfVf**Vf VrVfVf g Q 3 2 C 



C MC8032. 

SUBROUTINE GETVYC(SYCS, LG, YCSB,NYE,NYB,NYEG, VYC) MC8032: 

INTEGER*2 YCSB(NYE,NYB,NYEG) , VYC(NYE) , LGEX(6) MC8032: 

INTEGER SYCS MC8032' 

DATA LGEX/4,4,1,2,4,3/ MC8032! 

C --- L INDICATES LAST DIMENSION IN YCS EXPANSION TABLE MC8032( 

L=LGEX(LG) MC8032: 

C --- FIND TO WHICH YCS BLOCK SYCS BELONGS AND MAKE COPY IN VYC MC80321 

DO 10 J=1,NYB MC8032' 

DO 20 1=1, NYE MC8033C 

IF(SYCS .EQ. YCSB(I,J,L)) THEN MC8033I 

DO 30 K=1,NYE MC8033i 

VYC(K)=YCSB(K,J,L) MC8033; 

30 CONTINUE MC8033-: 
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RETURN MC803350 

ENDIF MC803360 

20 CONTINUE MC803370 

10 CONTINUE MC803380 

WRITE(6,*) ******** YCS NOT FOUND IN YCSB TABLE YCS=' ,SYCS MC803390 

END MC803400 

C MC803410 



SUBROUTINE READET(RETTBL,RETINV,MXX,NRET, SGRD, NSC,SVCMP, 

* NCSR.SCSRC, MG,LGRP,MGRP, MOSGR.NMS) 

C --- READ TABLE WITH ALL EXISTING COMBINATIONS FOR SELECTION CRITERIA 
C --- ACCEPT RECS WITH MATCHING PG,MG,CS ,SVC. ACCEPT ALL YCS 
INTEGER SVCMP(5), NSC, SVC 
INTEGER SCSRC(16),NCSR, CS 
INTEGER SGRD, PG 
INTEGER MOS,YCS 

INTEGER*2 MOSGR( 2 ,NMS) , MGRP(*) ,LGRP(*) 

INTEGER*2 RETTBL(MXX, 3) 

REAL RETINV(MXX), AI 
NRET=0 

DO 10 1=1,999999 

READ(10+SGRD,100,END=999) PG,MOS,YCS,SVC,CS, NRECS,AI 
IF(PG .NE. SGRD) GO TO 10 
IGR=IGFIND(MOS, MOSGR,NMS) 

LG=LGRP(IGR) 

IF(MGRP(LG) .NE. MG) GO TO 10 
DO 20 J=1,NSC 

IF(SVC .EQ. SVCMP(J)) GO TO 21 

20 CONTINUE 
GO TO 10 

21 CONTINUE 

DO 30 J=1,NCSR 

IF(CS .EQ. SCSRC(J)) THEN 

CALL ACCEPTC MOS , YCS , I GR , RETTBL , MXX , NRET , RETINV , AI ) 

GO TO 10 
ENDIF 

30 CONTINUE 
C 

10 CONTINUE 
999 CONTINUE 

IFCNRET .GT. MXX) THEN 

WRITE(6,*) ***** ERROR - TOO MANY RECORDS IN RETTBL* 

STOP 

ENDIF 

100 F0RMAT(I2,I4,I3,I2,I3,I4,F7. 2) 

END 



MC803430 

MC803440 

MC803450 

MC803460 

MC803470 

MC803480 

MC803490 

MC803500 

MC803510 

MC803520 

MC803530 

MC803540 

MC803550 

MC803560 

MC803570 

MC803580 

MC803590 

MC803600 

MC803610 

MC803620 

MC803630 

MC803640 

MC803650 

MC803660 

MC803670 

MC803680 

MC803690 

MC803700 

MC803710 

MC803720 

MC803730 

MC803740 

MC803750 

MC803760 

MC803770 

MC803780 

MC803790 

MC803800 

MC803810 

MC803820 






C MC803840 

SUBROUTINE ACCEPT(MOS ,YCS, IGR, RETTBL, MXX , NRET, RETINV ,AI) MC803850 

C --- ACCEPT ENTRY. ACCUMULATE IF ALREADY SAME COMBINATION IS PRESENT MC803860 

INTEGER MOS, YCS MC803870 

INTEGER*2 RETTBL(MXX, 3) MC803880 

REAL RETINV(MXX), AI MC803890 

DO 10 1=1, NRET MC803900 
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IFCMOS.EQ. RETTBLCI,!) .AND. YCS. EQ. RETTBL( 1 ,2) ) THEN MC803S 

RETINV(I)=RETINV(I) + AI MC803S 

RETURN MC803S 

ENDIF MC803S 

10 CONTINUE MC803S' 

C --- NEW COMBINATION MC8039 

NRET=NRET+1 MC8039 

RETTBL(NRET,1)=M0S MC8039: 

RETTBL(NRET,2)=YCS MC8035 

RETTBL(NRET,3)=IGR MC8040 

RETINV(NRET)=AI MC8040 

END MC8040: 

C MC8040: 






C 



c 



c 



c 



c 



c 



c 



SUBROUTINE GETPAR( AIMIN , NO , NMOS , SMOS , NYCS , SYCS , SGRD , 

* NSC,SVCMP, NCSR,SCSRC, IGR,MOSGR,NMS, ISFLAG) 

--- GET SELECTION CRITERIA FROM USER AND VALIDATE 
INTEGER SYCS(31), NYCS 
INTEGER SM0S(20), NMOS 
INTEGER SVCMP(5), NSC 
INTEGER SCSRC(16),NCSR 
INTEGER SGRD 
INTEGER*2 MOSGR(2,NMS) 

WRITE(5,*) ' ENTER THRESHOLD MIN. FOR AVERAGE INVENTORY’ 

READ(5,*) AIMIN 

WRITE(5,*) ' ENTER THRESHOLD MIN. FOR NUMBER OF CELLS' 

READ(5,*) NO 

WRITE(5,*) ' THRESHOLDS TO USE AIMIN, N0=' , AIMIN, NO 

WRITE(5,*) ' ENTER MOS (ONLY 1 ACCEPTED)' 

NMOS=l 

READ(5,*) SMOS(l) 

WRITE(6,*) ' MOS SELECTED;', SMOS(l) 

IGR=IGFIND(SMOS(l), MOSGR,NMS) 

WRITE(6,*) ' GROUP TO USE:’, IGR 
IF(IGR. EQ. 0) THEN 

WRITE(5,*) '**** ERROR - INVALID MOS SELECTED: ' ,SMOS( 1) 
STOP 

ENDIF 

WRITE(5,*) ’ ENTER YCS (ONLY 1 ACCEPTED)' 

NYCS=1 

READ(5,*) SYCS(l) 

WRITE(6,*) ’ YCS SELECTED:’, SYCS(l) 

WRITE (5,*) ’ ENTER GRADE' 

RE AD (5,*) SGRD 

WRITE(6,*) ' grade SELECTED', SGRD 

WRITE(5,*) ' ENTER NO. OF SVC. COMPS & ARRAY (1-3, 4=1+2, 5=ALL)' 
READ(5,*) NSC, (SVCMP(I), 1=1, NSC) 

--- EXPAND 4 TO 1,2 AND 5 TO 1,2,3 
DO 10 1=1, NSC 

IF(SVCMPd). EQ. 4 .OR. SVCMP( I). EQ. 5) THEN 
NSC=SVCMP(I)-2 



MC8040; 

MC8040 

MC8040; 

MC8040 

MC8040 

MC80411 

MC8041! 

MC8041. 

MC80411 

MC8041. 

MC8041; 

MC804K 

MC8041' 

MC80411 

MC8041< 

MC8042( 

MC8042: 

MC8042: 

MC8042: 

MC8042^ 

MC8042! 

MC8042« 

MC8042: 

MC8042£> 

MC8042S! 

MC8043C 

MC80431 

MC80432 

MC80433; 

MC80434: 

MC80435; 

MC80438: 

MC80437; 

MC80438;j 

MC8O439! 

MC80440: 

MC80441' 

MC80442 

MC80443 

MC80444 

MC80445 

MC80446 
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15 



10 

11 



c --- 



DO 15 J=1,NSC 
SVCMP(J)=J 
CONTINUE 
GO TO 11 
ENDIF 
CONTINUE 
CONTINUE 
WRITE(6,*) 



SERVICE COMPONENTS SELECTED', (SVCMP(I), 1=1, NSC) 



WRITE(5,*) 



25 



ENTER NO. COMM. SOURCES AND ARRAY (1-15, 16=ALL) ' 
READ(5,*) NCSR, (SCSRC(I), I=1,NCSR) 

IF 16 IS SELECTED THEN EXPAND ARRAY TO COVER ALL 1-15 
DO 20 1=1, NCSR 
IF(SCSRCd) .EQ. 16) THEN 
NCSR=15 

DO 25 J=1,NCSR 
SCSRC(J)=J 
CONTINUE 



GO TO 
ENDIF 

20 CONTINUE 
26 CONTINUE 
WRITE(5,*) 



26 



COMM. SOURCES SELECTED:', (SCSRC(I), 1=1, NCSR) 



MC804470 

MC804480 

MC804490 

MC804500 

MC804510 

MC804520 

MC804530 

MC804540 

MC804550 

MC804560 

MC804570 

MC804580 

MC804590 

MC804600 

MC804610 

MC804620 

MC804630 

MC804640 

MC804650 

MC804660 

MC804670 

MC804680 

MC804690 



MC804700 

FLAG TO DETERMINE WHICH OF SVC OR CS WILL BE USED AS 3RD DIMENSIONMC8047 10 
WRITE(5,*) 'SELECT 3RD DIM. TO USE: 0=N0NE, 1=SVC, 2=C0MM. SOURCE ' MC804720 



C --- 



* 



READ (5,*) ISFLAG 

WRITE INPUT PARAMETER INFO TO OUTPUT FILE 
WRITE(11,101) 'TEST CASE INPUT PARAMETERS:' 

WRITE(11,102) 'INVENTORY THRESHOLD= ' ,AIMIN, 

'THRESHOLD NO. OF CELLS= ' ,N0 
'MOS= ' ,SM0S(1),'YCS= ' ,SYCS(1) , 'GRADE= ' ,SGRD 
'SERVICE COMPONENTS= ' , ( SVCMP( I ) , 1=1 ,NSC) 

'COMM SOURCES= ' , ( SCSRC( I) , I=1,NCSR) 

'3RD DIMENSION= ', ISFLAG 



MC804730 
MC804740 
MC804750 
MC804760 
MC804770 
MC804780 
MC804790 
MC804800 
MC804810 
MC804820 
MC804830 
MC804840 
MC804850 
MC804860 
MC804870 

C MC804880 



V^T^ITE(ll,103) 
WRITE(11,104) 
WRITE(11,104) 
WRITE (6,*) 



101 

102 

103 

104 



F0RMAT(1X,A) 

F0RMAT(1X,A,F4. 1,7X,A,I2) 
FORMATC 1X,A,I3,2(5X,A,I2)) 
FORMATC 1X,A, 15(13)) 

END 



C MC804900 

SUBROUTINE GETMOS( SMOS , NMOS , MGX , LGX , MG , LG , IGR , MOSGR , LGRP , MGRP , MC8049 10 

* NMS,NG,NLG) MC804920 

C --- BUILD SMOS ARRAY BASED UPON EXPANSION MC804930 

INTEGER SM0S(30) MC804940 

INTEGER*2 MOSGR( 2 ,NMS) , LGRP(NG), MGRP(NLG) MC804950 

NMOS=0 MC804960 

IF(MGX .GT. 0) THEN MC804970 

C --- HAVE EXPANDED TO MAJOR MOS GROUP MC804980 

DO 10 1=1, NMS MC804990 

IGP=M0SGR(2,I) MC805000 

LGP=LGRP(IGP) MC805010 

IF(MGRP(LGP) .EQ. MG) THEN MC805020 
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NM0S=NM0S+1 MC805( 

SM0S(NM0S)=M0SGR(1,I) MC805( 

ENDIF MC805C 

10 CONTINUE MC805C 

RETURN MC805C 

ELSE IF(LGX . GT. 0) THEN MC805C 

C --- HAVE EXPANDED TO LARGE MOS GROUP MC805C 

DO 20 1=1, NMS MC805] 

IGP=MOSGR(2,I) MC8051 

IF(LGRP(IGP) .EQ. LG) THEN MC8051 

NMOS=NMOS+l MC805I 

SMOS(NMOS)=MOSGR(l,I) MC8051 

ENDIF MC8051. 

20 CONTINUE MC8051: 

RETURN MC8051: 

ELSE MC8051! 

C --- HAVE EXPANDED TO SMALL MOS GROUP MC80511 

DO 30 1=1, NMS MC8052. 

IF(MOSGR(2,I) .EQ. IGR) THEN MC8052; 

NMOS=NMOS+l MC8052; 

SMOS(NMOS)=MOSGR(l,I) MC8052. 

ENDIF MC8052: 

30 CONTINUE MC80521 

RETURN MC8052i 

ENDIF MC8052’. 

END MC8052i 

C MC8052'. 



Vf'A’VrVrVf Vryr****Vf Vf Vf Vf Vf Vr Vf VnVVrVf ** Vf ***** Vf ** *:?r*':V*****'A’ Vf Vf Vf Vf VfVc 8 0 5 3 1 



C MC8053: 

FUNCTION IGFINDCMOS, MOSGR,NMS) MC8053^ 

C FIND LOCATION OF MATCHING MOS IN GROUP TABLE. RETURN GROUP NO MC8053: 
INTEGER*2 M0SGR(2,NMS) MC8053^ 

DO 10 1=1, NMS MC8053! 

IF(M0SGR(1,I) .EQ. MOS) THEN MC8053( 

IGFIND=M0SGR(2,I) MC8053: 

RETURN MC80536 

ENDIF MC8053S 

10 CONTINUE MC8054C 

IGFIND=0 MC80541 

END MC80542: 

C MC80543 



************************************************************************[^jQ 80544 



MC80545; 

SUBROUTINE READER( DATA , INV , Y , MXX , NMOS , NYCSG , NYCSL , NYCSM , NSC , NCSR , MC8054& 
* NYR , SMOS , SYCSG , SYCSL , SYCSM , SGRD , SVCMP , SCSRC , NRC , PTRTBL , LGX , MGX , MC80547 

IGR,LG,MGRP,LGRP,M0SGR,NMS,NG,NLG,ICYCLE,NPT,PTBL,ISFLAG,SINV,SY)MC80548i 



REAL INV(MXX,NYR),Y(MXX,NYR), SINV(MXX,NYR) , SY(MXX,NYR) 
INTEGER*2 PTRTBL(MXX, 2), PTBL(MXX,3) 

INTEGER SYCSG(*), SYCSL(*), SYCSM(*) 

INTEGER SMOS(*), NMOS 
INTEGER SVCHP(*), NSC 
INTEGER SCSRC(*),NCSR 
INTEGER SGRD 

INTEGER TYPE , YCS , PG , MOS , SEX , CS , EDLV , SVC , MOS 1 , MOS2 , RACE 
INTEGER DATA(NYR) 

CHARACTER*? CITLS 



MC80549-: 

MC80550; 

MC805511I 

MC80552 

MC80553 

MC80554 

MC80555 

MC80556' 

MC80557I 

MC80558I 
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INTEGER*2 M0SGR(2,NMS) ,LGRP(NG) ,MGRP(NLG) 




MC805590 


c 








MC805600 


c 


--- REWIND DATA FILE AND RESET INV,Y IF 


CYCLING 


THRU READER 


MC805610 




IFdCYCLE .GT. 1) THEN 






MC805620 




REWIND 1 






MC805630 




DO 6 1=1, MXX 






MC805640 




DO 5 J=1,NYR 






MC805650 




INV(I,J)=0. 0 






MC805660 




Y(I,J)=0. 0 






MC805670 




SINV(I,J)=0.0 






MC805680 




SY(I,J)=0. 0 






MC805690 




5 CONTINUE 






MC805700 




6 CONTINUE 






MC805710 




ENDIF 






MC805720 


c 


--- READ RECORD AND STORE IN MATRIX 






MC805730 




ICR=0 






MC805740 




NRC=0 






MC805750 




NPT=0 






MC805760 




ICNT=0 






MC805770 




IYN0=0 






MC805780 




IYR=0 






MC805790 


c 








MC805800 




1 READ(1,101,END=999) TYPE,YCS,PG,MOS 


,SEX,CS, 


EDLV,SVC,MOSl,MOS2, 


MC805810 




* RACE, CITLS, DATA 






MC805820 




ICR=ICR+1 






MC805830 


c 


--- CHECK IF RECORD MEETS SELECTION CRITERIA. 


OTHERWISE REJECT. 


MC805840 


c 


--- COLLECT TYPES 0=INVENTORY, AND 1-5 


ALL LOSSES 


MC805850 




IFCTYPE. GT. 5) GO TO 999 






MC805860 


c 








HC805870 


c 


--- SCREEN FOR GRADE 






MC805880 




IF(PG .NE. SGRD) GO TO 1 






MC805890 


c 








MC805900 


c 


--- SCREEN FOR MOS 






MC805910 




IGP=IGFIND(MOS,MOSGR,NMS) 






MC805920 




IF(IGP. EQ. 0) GO TO 1 






MC805930 




LGP=LGRP(IGP) 






MC805940 




IF(MGX .GT. 0) THEN 






HC805950 


c 


HAVE EXPANDED TO MAJOR MOS GROUP 






MC805960 




IF(LGP .EQ. LG) THEN 






MC805970 




DO 10 I=1,NYCSL 






MC805980 




IF(YCS .EQ. SYCSL(I)) THEN 






MC805990 




IY=I 






MC806000 




GO TO 60 






MC806010 




ENDIF 






MC806020 




10 CONTINUE 






MC806030 




GO TO 1 






MC806040 




ELSE IF(MGRP(LGP) . EQ. MGX) THEN 






MC806050 




DO 20 I=1,NYCSM 






MC806060 




IF(YCS .EQ. SYCSM(D) THEN 






MC806070 




IY=I 






MC806080 




GO TO 60 






MC806090 




ENDIF 






MC806100 




20 CONTINUE 






MC806110 




GO TO 1 






MC806120 




ELSE 






MC806130 




GO TO 1 






MC806140 
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END IF 

ELSE IF(LGX . GT. 0) THEN 
C --- HAVE EXPANDED TO LARGE MOS GROUP 
IF(IGP .EQ. IGR) THEN 
DO 30 I=1,NYCSG 

IF(YCS .EQ. SYCSG(I)) THEN 
IY=I 

GO TO 60 
END IF 

30 CONTINUE 

GO TO 1 

ELSE IF(LGP .EQ. LGX) THEN 
DO 40 I=1,NYCSL 

IF(YCS .EQ. SYCSL(I)) THEN 
IY=I 

GO TO 60 
ENDIF 

40 CONTINUE 

GO TO 1 
ELSE 
GO TO 1 
ENDIF 
ELSE 

C --- HAVE EXPANDED TO SMALL MOS GROUP 
IF(IGP .EQ. IGR) THEN 
DO 50 I=1,NYCSG 

IF(YCS .EQ. SYCSG(I)) THEN 
IY=I 

GO TO 60 
ENDIF 

50 CONTINUE 

GO TO 1 
ELSE 

GO TO 1 
ENDIF 
ENDIF 

60 CONTINUE 
C 

DO 70 I=1,NM0S 

IF(MOS .EQ. SMOS(I)) THEN 
IM=I 

GO TO 80 
ENDIF 

70 CONTINUE 

WRITE(6,*) '*** ERROR IN MOS SCREENING ***' ,MOS 
WRITE(6,*) 'NMOS,SMOS=' ,NMOS , ( SMOS( I ) , 1=1 ,NMOS) 
GO TO 1 

--- SCREEN FOR SERVICE COMPONENT 
80 CONTINUE 

DO 90 1=1, NSC 

IF(SVC .EQ. SVCMP(I)) THEN 
IS=I 

GO TO 100 

END IF 
90 CONTINUE 



MC8061 

MC8061 

MC8061 

MC8061 

MC8061 

MC8062 

MC8062 

MC8062 

MC8062. 

MC8062 

MC8062: 

MC8062; 

MC806Z 

MC8062: 

MC8062 

MC8063 

MC8063: 

MC8063 

MC8063: 

MC8063: 

MC8063 

MC8063: 

MC8063’ 

MC8063: 

MC8063’ 

MC8064. 

MC8064. 

MC8064.- 

MC8064. 

MC80644 

MC8064> 

MC8064* 

MC8064 

MC8064. 

MC8064i 

MC8065( 

MC8065 

MC8065J 

MC8065; 

MC8065‘ 

MC8065.' 

MC8065« 

MC8065: 

MC8065! 

MC80655 

MC8066( 

MC8066: 

MC8066: 

MC8066: 

MC8066-! 

MC8066! 

MC8066( 

MC8066: 

MC8066? 

MC8066S 

MC8067C 
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c 

c --- 

100 



110 

c 

120 



C 

C 

999 



C --- 



C 



GO TO 1 

SCREEN FOR COMMISSIONING SOURCE 
CONTINUE 
DO 110 I=1,NCSR 
IF(CS .EQ. SCSRC(I)) THEN 
IR=I 

GO TO 120 

END IF 
CONTINUE 
GO TO 1 

CONTINUE 

RECORD ACCEPTED - INSTALL IT IN INV,Y,SINV,SY, PTRTBL AND PTBL 
ICNT=ICNT+1 
IF(ISFLAG. EQ. 1) THEN 
IW=IS 

ELSE IFdSFLAG. EQ. 2) THEN 
IW=IR 

ELSE 

IW=-99 

ENDIF 

MINV=GINV( PTRTBL, MXX.NRC, IM,IY,-99) 

MV=GINV(PTBL, MXX,NPT, IM, lY, IW) 

IF(TYPE.EQ. 0) THEN 

CALL INSINV(PTRTBL,MXX,NYR,NRC,MINV,IM,IY,-99,INV,DATA) 
CALL INSINVCPTBL, MXX , NYR , NPT , MV , IM,IY, IW, SINV,DATA) 

ELSE 

CALL INSY( MXX, NYR, MINV,Y, DATA) 

CALL INSY(MXX,NYR,MV, SY,DATA) 

IYR=IYR+1 

IF(MINV. EQ. 0) THEN 

WRITE(6,*) ’*** ERROR IN DATA BASE. LOSS W/0 INV. REC. ' 
WRITE(6,122) 'Y**:', MOS,YCS,PG,EDLV, SVC, RACE, 

(DATA(IT), IT=1,NYR) 

IYNO=IYNO+l 

ENDIF 

ENDIF 

GO TO 1 



CONTINUE 
WRITE(6,*) ’ ' 

WRITE(6,*) 'TOTAL RECORDS READ =’,ICR 

WRITE(6,*) 'TOTAL INV. MOS/YCS COMBINATIONS=' ,NRC 

WRITE(6,*) 'TOTAL INV. MOS/YCS/IW COMBINATIONS=' ,NPT 

WRITE(6,*) 'TOTAL RECORDS ACCEPTED =',ICNT 

WRITE(6,*) 'TOTAL LOSS RECORDS ACCEPTED =',IYR 

WRITE(6,*) 'TOTAL LOSS RECORDS NOT MATCHED =',IYNO 
TERMINATE IF NO DATA COLLECTED 
IF(NRC .EQ. 0) THEN 



WRITE(6,*) '**** NO DATA MEETS SELECTION REQS' 
STOP 

ENDIF 



MC806710 

MC806720 

MC806730 

MC806740 

MC806750 

MC806760 

MC806770 

MC806780 

MC806790 

MC806800 

MC806810 

MC806820 

MC806830 

MC806840 

MC806850 

MC806860 

MC806870 

MC806880 

MC806890 

MC806900 

MC806910 

MC806920 

MC806930 

MC806940 

MC806950 

MC806960 

MC806970 

MC806980 

MC806990 

MC807000 

MC807010 

MC807020 

MC807030 

MC807040 

MC807050 

MC807060 

MC807070 

MC807080 

MC807090 

MC807100 

MC807110 

MC807120 

MC807130 

MC807140 

MC807150 

MC807160 

MC807170 

MC807180 

MC807190 

MC807200 

MC807210 

MC807220 

MC807230 

MC807240 

MC807250 

MC807260 
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101 F0RMAT(3I2,I3,I1,I2,2I1,2I3,I1,A7, IX, 1014) MC807; 

121 FORMAT(A8,13I6) MC807: 

122 FORMAT(A8,7I6, 5X, 1216) MC807; 

131 F0RMAT(I4, 216) MC807.' 

132 F0RMAT(I4, 316, 10F7. 2) MC807: 

END MC807: 

C MC807: 

C MC807- 

FUNCTION GINV(PTBL, MXX,NPT, IM,IY,IW) MC807: 

--- FIND LOCATION OF INVENTORY ENTRY FOR MOS , YCS , SVC/CS COMBINATIONS MC807- 

--- 3RD DIMENSION CHECKED ONLY IN CASE IW>0 MC807i 

INTEGER*2 PTBL(MXX, *) MC807: 

DO 10 1=1, NPT MC8074 

IF(PTBL(I, 1) .EQ. IM .AND. MC8074 

* PTBLd, 2) .EQ. lY ) THEN MC8074 

IFCIW.LT.O .OR. (IW.GT.O .AND. PTBL(I, 3).EQ. IW)) THEN MC8074 

GINV=I MC8074 

RETURN MC8074 

ENDIF MC8074 

ENDIF MC8074 

10 CONTINUE MC8074 

GINV=0 MC8074 

END MC8075: 

C MC8075 






C MC8075; 

SUBROUTINE INSINV(PT,MXX,NYR,N,K, IM, lY, IW, INV,DATA) MC8075: 

C --- ACCUMM INTO KTH ENTRY. INSTALL IN POINTER TABLE IF NOT PRESENT MC8075': 

REAL INV(MXX, NYR) MC8075; 

INTEGER*2 PT(MXX, *) MC8075: 

INTEGER DATA(NYR) MC8075^ 

IF(K .EQ. 0) THEN MC8075 

C ... add new entry MC8076- 

N=N+1 MC8076: 

IF(N .GT. MXX) THEN MC8076; 

WRITE(6,*) '*** error - TOO MANY INV. COMBINATIONS' ,N MC8076: 

STOP MC8O760 

ENDIF MC8076J 

K=N MC8076« 

PT(K, 1)=IM MC8076) 

PT(K, 2)=IY MC8076i 

IF(IW. GT. 0) PT(K, 3)=IW MC8076i 

ENDIF MC8077( 

DO 130 IT=1,NYR MC8077: 

INV(K,IT)=INV(K,IT) + . 25*FLOAT(DATA( IT)) MC8077; 

130 CONTINUE MC8077; 

END MC8077< 

C MC8077f 

Vf * ** Vf yryc*Vc**VcVr**yr**Vc****yr yr**Vc*VrV?** **** Vr* Vf** Vf Vf Vf 0 7 7 ( 

C MC8077; 

SUBROUTINE INSY(MXX, NYR, K,Y, DATA) MC8077{ 

C --- ACCUMM INTO KTH ENTRY FOR LOSS MC80775 

REAL Y(MXX, NYR) MC8078C 

INTEGER DATA(NYR) MC8078] 

IF(K .EQ. 0) RETURN MC80782 
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DO 10 IT=1,NYR MC807830 

Y(K,IT)=Y(K,IT) + DATA(IT) MC807840 

10 CONTINUE MC807850 

END MC807860 

C MC807870 



*************Vf *Vf******************Vf***********Vf********^f***********t^C s 0 7 8 8 0 



C 

SUBROUTINE AGGREG( INV , Y , MXX , NYR , SMOS , SYCSG , 

* NRC,NRCOLD,PTRTBL,INDX,AVINV, AIMIN,MKG) 

C --- COMP. AVERAGE INV. & SORT 

REAL INV(MXX, NYR), Y(MXX, NYR), AVINV(MXX) 

INTEGER*2 PTRTBL(MXX, 2), INDX(MXX) ,MKG(MXX) 

INTEGER SYCSGC*), SMOS(*) 

REAL*8 TINV,TY 

--- RESET MKG (NECESSARY WHEN CYCLING THRU AGGPCT VALUES) 

DO 10 1=1, MXX 

MKG(I)=0 
10 CONTINUE 
TINV=0 
TY=0 

DO 100 1=1, NRC 

--- FIX INV. ENTRIES LOWER THAN CORRESP. LOSSES & COMP. AVG INV. 
AI=0 

DO 201 J=1,NYR 
TINV=TINV+INV(I,J) 

TY= TY+ Y(I,J) 

IF( INV( I , J) . LT. Y( I , J) ) INV( I , J)=Y( I , J) 

AI=AI+INV(I,J) 

201 CONTINUE 

AVINV(I)=AI/NYR 
INDX(I)=I 
100 CONTINUE 

WRITE(6,*) '== TOTAL INV,Y=' ,TINV,TY 



--- SORT ASCENDING BY AVG INVENTORY 
CALL SORT2(AVINV,INDX,NRC) 

C 

NS 1=0 

C ... display TABLE IN SORT SEQUENCE 

CALL DSPTBLC INV , Y , AVINV , PTRTBL , INDX , AIM IN , NRC , MKG , MXX , NYR , 

* SYCSG, SMOS ) 

C 

DO 200 K=NRC,1,-1 
IF(AVINV(K) .GE. AIMIN) THEN 
C --- MARK AS MEMBER OF SET SO 

MKG(K)=32767 

ELSE 

C --- INITIAL COUNT OF MEMBERS OF SET SI 

NS1=K 
GO TO 202 

ENDIF 

200 CONTINUE 
202 CONTINUE 

C --- DO AGGREGATIONS WITHIN SET SI UNTIL NO MORE POSSIBLE (KF GE 0) 



MC807890 

MC807900 

MC807910 

MC807920 

MC807930 

MC807940 

MC807950 

MC807960 

MC807970 

MC807980 

MC807990 

MC808000 

MC808010 

MC808020 

MC808030 

MC808040 

MC808050 

MC808060 

MC808070 

MC808080 

MC808090 

MC808100 

MC808110 

MC808120 

MC808130 

MC808140 

MC808150 

MC808160 

MC808170 

MC808180 

MC808190 

MC808200 

MC808210 

MC808220 

MC808230 

MC808240 

MC808250 

MC808260 

MC808270 

MC808280 

MC808290 

MC808300 

MC808310 

MC808320 

MC808330 

MC808340 

MC808350 

MC808360 

MC808370 



65 



KF=-1 

C --- DO WHILE KF<0 

300 IF(KF. GE. 0) GO TO 310 

CALL AGG1(AVINV,INDX,MKG,NS1,INV,Y,MXX,NYR,AIMIN,KF) 

GO TO 300 
310 CONTINUE 

C --- DISPLAY TABLE AFTER 1ST AGGREGATION 

CALL DSPTBL(INV,Y,AVINV,PTRTBL,INDX,AIMIN,NRC,MKG,MXX,NYR, 

* SYCSG,SMOS ) 

IF(NSl.EQ.NRC) THEN 

WRITE(6,*) ****** SET SO EMPTY. NO CELLS ABOVE THRESHOLD' 
STOP 

END IF 

C --- DO AGGREGATIONS FROM SET SI INTO SET SO UNTIL NO MORE POSSIBLE 
KF=1 

C --- DO WHILE KF>0 

320 IF(KF. LE.O) GO TO 330 

CALL AGG2(AVINV,INDX,MKG,NS1,NRC,INV,Y,MXX,NYR, KF) 

GO TO 320 
330 CONTINUE 

C --- DISPLAY TABLE AFTER 2ND AGGREGATION 

CALL DSPTBL(INV,Y,AVINV,PTRTBL,INDX,AIMIN,NRC,MKG,MXX,NYR, 

* SYCSG,SMOS ) 

C --- MOVE VALUES GE AIMIN TO BEGINNING OF ARRAYS 

CALL CMPRS ( INV , Y , MXX , NYR , NRC , NRCOLD , AIMIN , AVINV) 

C --- DISPLAY TABLE AFTER MOVING VALUES. 

DO 400 K=1,NRC 

WRITE(6,122)K,AVINV(K), ( INV(K, J) , J=1 ,NYR) 

WRITE(6,123) ( Y(K,J),J=1,NYR) 

400 CONTINUE 

122 FORMAT(/I5,14X,F8. 3, 6X, 10F7. 2) 

123 FORMATC 33X, 10F7.2) 

END 

^ Vr Vr *sV Vr Vc Vf 

SUBROUTINE AGG1( AVINV, INDX,MKG,NS1, INV, Y, MXX, NYR, AIMIN, KF) 

C --- DO ONE PASS OF AGGREGATION 

REAL INV(MXX, NYR), Y(MXX, NYR), AVINV(MXX) 

INTEGER*2 INDX(MXX) ,MKG(MXX) 

KF=0 

CI=0 

DO 10 I=NS1,1,-1 
IF(MKG(I).EQ. 0) THEN 
IF(KF. EQ. 0) THEN 

C --- THIS WILL BE THE COLLECTING CELL 



MC808: 

MC808; 

MC8084 

MC808^ 

MC808^ 

MC808^ 

MC808^ 

MC8084 

MC808-! 

MC808-! 

MC808^ 

MC808^ 

MC808I 

MC808! 

MC808' 

MC8085 

MC808S 

MC8085 

MC8085 

MC8085 

MC808S 

MC8085 

MC8086 

MC8086 

MC8086 

MC8086 

MC8086 

MC8086 

MC8086 

MC8086 

MC8086 

MC8086 

MC8087 

MC8087 

MC8087 

MC8087 

MC8087; 

MC8087: 

MC8087; 

MC8087' 

MC8087! 

MC8087; 

MC8088= 

MC8088] 



KF=I MC8088i 

CI=AVINV(I) MC8088! 

ELSE MC8088! 

IF(CI+AVINV(I).LT. AIMIN) THEN MC8088; 

C --- ACCUM. WITH CELL KF TEMPORARILY. SET TEMP. POINTER -KF MC8088' 

CI=CI+AVINV(I) MC8088 

MKG(I)=-KF MC8088. 

ELSE MC8088' 

C --- FIND SMALLEST CELL TO ADD MC8089i 

CALL AGG1A( AVINV, MKG, I, Cl, AIMIN, KF, MXX) MC8089 

ENDIF MC8089: 



I 
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IF(CI.GE. AIMIN) THEN 


MC808930 


c 


— 


MAKE THIS AGGREGATION PERMANENT AND EXIT 


MC808940 






AVINV(KF)=CI 


MC808950 






CALL AGG1B(INDX,MKG,KF,INV,Y,NYR,MXX) 


MC808960 






NS1=NS1-1 


MC808970 






MKG(KF)=32767 


MC808980 






KF=-1 


MC808990 






RETURN 


MC809000 






ENDIF 


MC809010 






ENDIF 


MC809020 






ENDIF 


MC809030 




10 


CONTINUE 


MC809040 


c 






MC809050 






IFCKF.EQ. 0) RETURN 


MC809060 


c 


— 


CLEAR TEMPORARY POINTERS LEFT. THIS WAS AN UNSUCCESSFUL AGGREG. 


MC809070 






DO 20 1=1, NS 1 


MC809080 






IF(MKG(I). LT. 0) MKG(I)=0 


MC809090 




20 


CONTINUE 


MC809100 






END 


MC809110 


c 




r*********** 


MC809120 






SUBROUTINE AGG1A(AVINV,MKG, ILAST, Cl , AIMIN, KF,MXX) 


MC809130 


c 


— 


FIND SMALLEST CELL TO ADD AND SET TEMPORARY POINTER 


MC809140 






REAL AVINV(MXX) 


MC809150 






INTEGER*2 MKG(MXX) 


MC809160 






DO 10 1=1, ILAST 


MC809170 






IF(MKG(I).EQ. 0) THEN 


MC809180 






IF(CI+AVINV(I) .GE. AIMIN) THEN 


MC809190 






CI=CI+AVINV(I) 


MC809200 






MKG(I)=-KF 


MC809210 






RETURN 


MC809220 






ENDIF 


MC809230 






ENDIF 


MC809240 




10 


CONTINUE 


MC809250 






WRITE(6,*) '*** ERROR IN AGGIA. NO VALUE FOUND ***’ 


MC809260 






STOP 


MC809270 






END 


MC809280 


c 


^V'5V^Wc’}WrfV'5V^V’sWc"5V*’)V/V 


MC809290 






SUBROUTINE AGGIB ( INDX , MKG , KF , INV , Y , NYR , MXX ) 


MC809300 


c 


— 


MAKE AGGREGATION PERMANENT 


MC809310 






REAL INV(MXX, NYR), Y(MXX, NYR) 


MC809320 






INTEGER*2 INDX(MXX) ,MKG(MXX) 


MC809330 






K=INDX(KF) 


MC809340 






DO 10 I=1,KF-1 


MC809350 






IF(MKG(I) .LT. 0) THEN 


MC809360 






IF(MKG(I).NE. -KF) STOP 777 


MC809370 






MKG(I)=KF 


MC809380 






L=INDX(I) 


MC809390 






DO 20 J=1,NYR 


MC809400 






INV(K,J)=INV(K,J)+INV(L,J) 


MC809410 






Y(K,J)= Y(K,J)+ Y(L,J) 


MC809420 




20 


CONTINUE 


MC809430 






ENDIF 


MC809440 




10 


CONTINUE 


MC809450 






END 


MC809460 


c 




jV Vc Vc V? Vc 


MC809470 
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o o 



c --- 
c --- 



c --- 



SUBROUTINE AGG2( AVINV, INDX,MKG,NS1 ,NRC , INV, Y,MXX,NYR, KF) 

DO ONE PASS OF AGGREGATION FROM SET SI TO SET SO 

ON EACH PASS ONE ELEMENT OF SI IS TAKEN & ADDED TO SMALLEST OF SO 

REAL INV(MXX, NYR) , Y(MXX, NYR), AVINV(MXX) 

INTEGER*2 INDX(MXX) ,MKG(MXX) 

KF=0 

FIND ELEMENT OF SI (ONLY THOSE WITH POINTER MKG(I)=0) 

DO 10 1=1, NSl 

THEN 



10 

12 

C --- 



IF(MKG(I).EQ. 0) 

KF=I 

GO TO 12 

ENDIF 

CONTINUE 

CONTINUE 

IF KF STILL 0 THEN NO MORE ELEMENTS IN SI LEFT 
IF(KF.EQ. 0) RETURN 



--- FIND SMALLEST ELEMENT OF SO AND ADD TO IT. ONLY WITH MKG(I)=32767 
ISM=NRC 

SMALL=AVINV( ISM) 

DO 20 1=1, NRC 
IF(MKG(I). EQ. 32767) THEN 
IF( AVINV( I ) . LT. SMALL) THEN 
ISM=I 

SMALL=AVINV(I) 

ENDIF 
ENDIF 

20 CONTINUE 

--- JOIN ELEMENT KF TO ELEMENT ISM 
AVINV(ISM)=AVINV(ISM) + AVINV(KF) 

MKG(KF)=ISM 
L=INDX(KF) 

K=INDX( ISM) 

DO 30 J=1,NYR 
INV(K,J)=INV(K,J)+INV(L,J) 

Y(K,J)= Y(K,J)+ Y(L,J) 

CONTINUE 
END 

Q Vr * Vr* ** Vc”>V* * 

SUBROUTINE CMPRS (INV , Y , MXX , NYR , NRC , NRCOLD , AIMIN , AVINV ) 

REAL INV(MXX, NYR), Y(MXX, NYR), AVINV(MXX) 



30 



C --- 



20 



10 



COMPRESS INV,Y IN PLACE. 

NRCOLD=NRC 

NRC=0 

DO 10 1=1, NRCOLD 
AI=CAINV( INV, I , MXX, NYR) 
IF(AI .GE. AIMIN) THEN 
TRANSFER ACTIVE CELL 
NRC=NRC+1 
AVINV(NRC)=AI 
DO 20 J=1,NYR 
INV(NRC,J)=INV(I,J) 
Y(NRC,J)= Y(I,J) 
CONTINUE 
ENDIF 
CONTINUE 



MOVE ALL ROWS GE AIMIN TO TOP 



I ---> NRC 



MC809il 

MC8094I 

MC8095i 

MC8095' 

MC8095. 

MC8095 

MC8095: 

MC8095 

MC8095; 

MC8095: 

MC8095’: 

MC8095; 

MC8096; 

MC8096: 

MC8096; 

MC8096: 

MC8096; 

MC8096; 

MC8096-, 

MC8096i 

MC8096( 

MC8096'. 

MC8097( 

MC8097: 

MC8097,' 

MC8097: 

MC8097^ 

MC8097i 

MC8097( 

MC8097: 

MC8097J 

MC8097' 

MC8098C 

MC8098] 

MC8098: 

MC8098: 

MC8098^ 

MC8098S 

MC80986 

MC80987; 

MC80988: 

MC80989' 

MC80990: 

MC809911 

MC80992; 

MC80993: 

MC809944 

MC80995) 

MC80996 I 

MC80997', 

MC80998' 

MC809991 

MC81000I 

MC8100K 

MC81002( 

MC81003( 
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END 

Q ** ** ****** ***** 

FUNCTION CAINV(INV,I,MXX,NYR) 

REAL INVCMXX, NYR) 

C --- COMPUTE AVERAGE INVENTORY FOR ROW I 
CAINV=0 
DO 10 J=1,NYR 
CAINV=CAINV+INV( I , J) 

10 CONTINUE 

CAINV=CAINV/NYR 

END 

Q *************** 

SUBROUTINE DSPTBL( INV , Y , A VINV , PTRTBL , INDX , A IMIN , NRC , MKG , MXX , NYR , 
* SYCSG,SMOS) 

C --- DISPLAY TABLE IN SORT SEQUENCE 

REAL INVCMXX, NYR), Y(MXX, NYR), AVINV(MXX) 

INTEGER*2 PTRTBL(MXX, 2), INDX(MXX) ,MKG(MXX) 

INTEGER SYCSGC*) 

INTEGER SMOSC*) 

INTEGER IATT(2) 

CHARACTER*! STI 
WRITE(6,121) 

WRITEC6,*) 'INV. THRESHOLD MIN. VALUE=' ,AIMIN 



WRITE (6,*) 



I INDX AVG 



MKG 



INVENTORY/ LOSSES ' 



DO 200 K=1,NRC 
STI=' ' 

I=INDX(K) 

AI=AVINV(K) 

IFCAI .LT. AIMIN) STI=’$' 
lATTC 1 )=SMOS ( PTRTBLC 1,1)) 

IATT(2)=SYCSG(PTRTBL(I,2)) 

WRITE(6,122)K,I,AI,MKG(K),STI,(INV(I,J),J=1,NYR),(IATT(J),J=1,2) 
* PTRTBLC 1,1), PTRTBLC I, 2) 

WRITEC6,123) C YCI,J),J=1,NYR) 

200 CONTINUE 
C 

121 FORMATC///) 

122 F0RMATC/2I5,F8. 3,I9,1X,A2, 10F7.2, 5X, 615) 

123 FORMATC 30X, 10F7. 2) 

END 

Vc Vc Vc Vc r Vc Vc Vc Vc Vc c 

SUBROUTINE S0RT2CY, INDX, N) 

C --- INPLACE SORT USING SHELL ALGORITHM ******** 

C --- SORTS ON Y AND DOES SAME REORDERING ON INDEXES INDX 
REAL YCN),TEMP 
INTEGER GAP 

INTEGER*2 INDXCN), ITEMP 
LOGICAL EXCH 
C 

GAP=CN/2) 

5 IF C. NOT. CGAP. NE. 0)) GO TO 500 
10 CONTINUE 

EXCH=. TRUE. 

K=N-GAP 
DO 200 1=1, K 



MC810040 
MC810050 
MC810060 
MC810070 
MC810080 
MC810090 
MC810100 
MC810110 
MC810120 
MC810130 
MC810140 
MC810150 
MC810160 
MC810170 
MC810180 
MC810190 
MC810200 
MC810210 
MC810220 
MC810230 
MC8 10240 
MC810250 
MC810260 
MC810270 
MC810280 
MC810290 
MC810300 
MC810310 
MC810320 
MC8 10330 
MC8 10340 
MC8 10350 
,MC810360 
MC810370 
MC810380 
MC8 10390 
MC810400 
MC810410 
MC8 10420 
MC8 10430 
MC8 10440 
MC810450 
MC810460 
MC810470 
MC810480 
MC8 10490 
MC810500 
MC810510 
MC810520 
MC810530 
MC8 10540 
MC8 10550 
MC810560 
MC810570 
MC810580 
MC810590 
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KK=I+GAP 


MC8106 




IF(.NOT. (Y(I). GT. Y(KK))) GO TO 100 


MC8106 




TEMP=Y(I) 


MC8106 




Y(I)=Y(KK) 


MC8106 




Y(KK)=TEMP 


MC8106 




ITEMP=INDX(I) 


MC8106 




INDX(I)=INDX(KK) 


MC8106 




INDX(KK)=ITEMP 


MC8106' 




EXCH=. FALSE. 


MC8106: 


100 


CONTINUE 


MC8106: 


200 


CONTINUE 


MC8107; 




IF (.NOT. (EXCH)) GO TO 10 


MC8107i 




GAP=(GAP/2) 


MC8107; 




GO TO 5 


MC8IO7; 


500 


CONTINUE 


MC8107' 




RETURN 


MC8107' 




END 


MC8107^ 


c 




MC8107 


Vf9V*^V^V?Vyc7»V^V‘5V:V*^V7V* V?‘>V“V^V:iV^V7?c Vf *‘>V*>V‘>V*‘V^VVnV:iV:^:iV‘>V‘jVT?f ■:V‘>V‘>V^V7V7Sr*';?f *‘>V'>V‘>V*‘>V*>V‘>V 8107 


c 




MC8107' 




SUBROUTINE BKDOWN( PTBL , NPT , PTRTBL , NRC , INDX , MKG , MXX , MXY , 


MC8108I 




SINV,SY,INV,Y,BKTBL,NBK ) 


MC8108) 


c --- 


BREAKDOWN AGGREGATED VALUES BY THE 3RD DIMENSION SVC/CS 


MC8108: 




REAL INV(MXX,MXY),Y(MXX,MXY), SINV(MXX,MXY) ,SY(MXX,MXY) 


MC8108: 




INTEGER*2 PTRTBL(MXX, 2), INDX(MXX), MKG(MXX) 


MC8108^ 




INTEGER*2 PTBL(MXX, 3), BKTBL(MXX,3) 


MC8108! 




REAL*8 TINV,TY 


MC8108( 




NBK=0 


MC8108: 


c — - 


TRAVERSE MKG ARRAY AND BUILD BKTBL 


MC8108? 




DO 10 1=1, NRC 


MC8108S 




IF(MKG(I). NE. 32767) THEN 


MC8109C 




ICELL=MKG(I) 


MC8109] 




ELSE 


MC8109J 




ICELL=I 


MC8109C- 




ENDIF 


MC81094 




IX=INDX(I) 


MC8109S 




IM=PTRTBL(IX,1) 


MC81096 




IY=PTRTBL(IX,2) 


MC81097. 




C ALL BLDBK( ICE LL , IM , I Y , PTBL , NPT , MXX , BKTBL , NBK) 


MC81098 


10 


CONTINUE 


MC8IO99; 


c --- 


DISPLAY BKTBL PRIOR TO SORTING 


MC8IIOO: 




WRITE(6,101) (I,(BKTBL(I,J),J=1,3), 1=1, NBK) 


MC81101i 




CALL SORT3( BKTBL, NBK, MXX) 


MC8110Z 




WRITE(6,101) (I,(BKTBL(I,J),J=1,3), 1=1, NBK) 


MC8IIO3: 


c --- 


SUMMARIZE SINV,SY INTO INV,Y FOR MATCHING ENTRIES IN BKTBL 


MC81104( 




CALL SUMBK( BKTBL , NBK , MXX , S INV , S Y , INV , Y , MXY) 


MC8IIO5; 




WRITE(6,102) (I,(INV(I,J) ,J=1,MXY),(BKTBL(I,J),J=1,2), 1=1, NBK) 


MC81106: 




WRITE(6,102) (I,( Y(I,J),J=1,MXY),(BKTBL(I,J),J=1,2), 1=1, NBK) 


MC81107’ 


101 


F0RMAT(I4, 316) 


MC81108 


102 


F0RMAT(I4, 10F7. 2,10X,2I4) 


MC81109 


103 


F0RMAT(/I5,10F7. 2) 


MC81110 


104 


FORMATC 5X,10F7.2) 


MC81111' 




END 


MC81112' 




MC8III31 




SUBROUTINE BLDBK( I CELL , IM , I Y , PTBL , NPT , MXX , BKTBL , NBK) 


MC81114t 
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INTEGER*2 PTBL(MXX, 3), BKTBL(MXX,3) 

C --- RECORD ALL ENTRIES IN PTBL WITH MATCHING IM,IY IN BKTBL 
DO 10 1=1, NPT 

IF(PTBL(I,1). EQ. IM .AND. PTBL( 1 ,2). EQ. lY) THEN 
C --- INSTALL WITH CELL ID, IW & POINTER 

NBK=NBK+1 

BKTBL(NBK, 1)=ICELL 
BKTBL(NBK,2)=PTBL(I,3) 

BKTBL(NBK,3)=I 

ENDIF 

10 CONTINUE 
END 

Q iV^^V^V’jV*****^^**^^'** 

SUBROUTINE S0RT3(T,N,MXX) 

C --- INPLACE SORT USING SHELL ALGORITHM ******** 

C --- SORTS ON 1ST 2 COLS. OF T & DOES SAME REORDERING ON 3RIr'C0LUMN 
INTEGER*2 T(MXX,3), ITEMP 
INTEGER GAP 
LOGICAL EXCH 



5 

10 



GAP=(N/2) 

IF (GAP. EQ. 0) GO TO 500 
CONTINUE 

EXCH=. FALSE. 

K=N-GAP 
DO 200 1=1, K 
KK=I+GAP 
IF(T(I,1).GT. T(KK,1) 



OR. 



MC811150 

MC811160 

MC811170 

MC811180 

MC811190 

MC811200 

MC811210 

MC811220 

MC811230 

MC811240 

MC811250 

MC811260 

MC811270 

MC811280 

MC811290 

MC811300 

MC811310 

MC811320 

MC811330 

MC811340 

MC811350 

MC811360 

MC811370 

MC811380 

MC811390 

MC811400 

MC811410 

MC811420 



(T(I,1).EQ.T(KK,1) .AND. T(I,2). GT. T(KK,2)) ) THEN MC811430 



IT1=T(I,1) 

IT2=T(I,2) 

IT3=T(I,3) 

T(I,1)=T(KK,1) 

T(I,2)=T(KK,2) 

T(I,3)=T(KK,3) 

T(KK,1)=IT1 
T(KK,2)=IT2 
T(KK,3)=IT3 
EXCH=. TRUE. 

ENDIF 

200 CONTINUE 

IF (EXCH) GO TO 10 
GAP=(GAP/2) 

GO TO 5 
500 CONTINUE 
RETURN 
END 

SUBROUTINE SUMBK( BKTBL , NBK , MXX , S INV , SY , INV , Y , MXY) 

C --- CREATE AGGREGATED ARRAYS INV,Y FROM CELL & 3RD DIM. INFO. 
REAL INV(MXX,MXY),Y(MXX,MXY), SINV(MXX,MXY) , SY(MXX,MXY) 
INTEGER*2 BKTBL(MXX,3) 

REAL*8 TINV,TY 

IP=0 

Il=-1 

I2=-l 



MC811440 
MC811450 
MC811460 
MC811470 
MC811480 
MC811490 
MC811500 
MC811510 
MC811520 
MC811530 
MC811540 
MC811550 
MC811560 
MC811570 
MC811580 
MC811590 
MC811600 
MC811610 
MC811620 
MC811630 
IN BKTBLMC8 11640 
MC811650 
MC811660 
MC811670 
MC811680 
MC811690 
MC811700 
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TINV=0 


MC81i; 






TY=0 


MCBli: 






DO 10 1=1, NBK 


MC81i: 






IF(BKTBL(I,1).NE. 11 .OR. BKTBL( 1 ,2) . NE. 12) THEN 


MC81i: 


c 


— 


CHANGE OF CELL,IW IDENTIFIERS 


MC81i: 






IP=IP+1 


MC811'/ 






I1=BKTBL(I,1) 


MC811/ 






I2=BKTBL(I,2) 


MC81i; 






DO 15 J=1,MXY 


MC811*/ 






INV(IP,J)=0 


MC8118 






Y(IP,J)=0 


MC8118 




15 


CONTINUE 


MC8118 






BKTBL(IP,1)=I1 


MC811S 






BKTBL(IP,2)=I2 


MC8118 






ENDIF 


MC8118 


c 


— 


ACCUMULATE 


MC8118: 






I3=BKTBL(I,3) 


MC8118 






DO 20 J=1,MXY 


MC8118 






INV( IP , J)=INV( IP , J)+SINV( 13 , J) 


MC8118' 






Y(IP,J)= Y(IP,J)+ SY(I3,J) 


MC8119; 






TINV=TINV+SINV(I3,J) 


MC8119^ 






TY= TY+ SY(I3,J) 


MC8119I 




20 


CONTINUE 


MC8119) 




10 


CONTINUE 


MC8119. 






WRITE(6,*) ' TOTAL INV,Y AFTER BREAKDOWN-’ ,TINV,TY 


MC8119.I 


c 






MC8119I 






NBK=IP 


MC8119‘ 


c 


— 


FIX INV. ENTRIES LOWER THAN CORRESP. LOSSES 


MC8119) 






DO 40 1=1, NBK 


MC81195 






DO 30 J=1,MXY 


MC8120( 






IF( INV( I , J). LT. Y( I , J) ) INV( I , J)=Y( I , J) 


MC8120: 




30 


CONTINUE 


MC8120: 




40 


CONTINUE 


MC8120: 






END 


MC8120^ 
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on onjj-noooooooo 



C. ESTIMATION SUBROUTINES 



C 



SUBROUTINE MC87BZ(INV,Y,NRC,NYR,XTB,VXTB,XEB,A,MXX,MXY) 
--- CONDUCTS FIRST FIVE ESTIMATION METHODS 

REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) 

REAL XTB(MXX), VXTB(MXX), A(MXX) 

CALL EBTS1(INV,Y,NRC,NYR,XTB,VXTB,XEB,MXX,MXY) 

PRINT *, 'COMPLETED EBTSl' 

CALL EBTS2(INV,Y,NRC,NYR,XTB,VXTB,XEB,MXX,MXY) 

PRINT *, 'COMPLETED EBTS2' 

CALL EB0S1(INV,Y,NRC,NYR,XTB,VXTB,XEB,MXX,MXY) 

PRINT *, 'COMPLETED EBOSl' 

CALL EB0S2( INV , Y , NRC , NYR , XTB , VXTB , XEB , MXX , MXY) 

PRINT *, 'COMPLETED EB0S2' 

CALL EMTS(INV,Y,NRC,NYR,XTB,VXTB,XEB,A,MXX,MXY) 

PRINT 'COMPLETED EMTS' 

END 



MC800010 

MC800020 

MC800030 

MC800040 

MC800050 

MC800060 

MC800070 

MC800080 

MC800090 

MC800100 

MC800110 

MC800120 

MC800130 

MC800140 

MC800150 

MC800160 

MC800170 

MC800180 

MC800190 

MC800200 

MC800210 

MC800220 

MC800230 



*************>V***'>V****Vr**‘***********yr****^***Vr***Vr)V*****')V****')V*;V[^Q 800240 



SUBROUTINE EBTS 1 ( I NV , Y , NRC , NYR , XTB , VXTB , XEB , MXX , MXY ) 

--- TRANSFORMED SCALE, TIME INDEPENDENT VARIANCE METHOD 
REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) 

REAL XTB(MXX), VXTB(MXX) 

REAL MAXL,MINL,L,MAXCHI,MINCHI 
INTEGER T, VYR 

DATA AA/1. 6835/, B1/-.8934/, B2/.9881/ 

MAXL= -1000. 0 
MINL= 1000. 0 
SUML= 0. 0 
KLSUM=0 

MAXCHI= -1000.0 
MINCHI= 1000.0 
SUMCHI= 0.0 
KSUM=0 
SUMMAD=0. 0 
KPSUM=0 

WRITE(11,32)' ' 

WRITE(11,21)'EMP BAYES TRANS SCALE - TIME DEP VAR: ' 

VRITECl 1,21) 'MEAN ABSOLUTE DEVIATION (ORIG SCALE):' 

WRITEC 11, 28) 'FRACTION CELLS' , 'FRACTION MAD' 

WRITE( 1 1 , 29 ) ' VALID YR ' , ' K ' , ' WITH UNDERAGE ' , ' FROM UNDERAGE ' 
WRITE(11,30) 

--- LOOP THROUGH VALIDATION YEARS 
DO 280 VYR=1, NYR 
--- LOOP THROUGH CELLS 
DO 260 IN=1, NRC 



MC800250 
MC800260 
MC800270 
MC800280 
MC800290 
MC800300 
MC800310 
MC800320 
MC800330 
MC800340 
MC800350 
MC800360 
MC800370 
MC800380 
MC800390 
MC800400 
MC800410 
MC800420 
MC800430 
MC800440 
MC800450 
MC800460 
'MAD' MC800470 
MC800480 
MC800490 
MC800S00 
MC800510 
MC800520 
MC800530 
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c 





T=0 


MC8005 




SUMXT=0 


MC8005; 




SUMVAR=0 


MC8005' 


— 


LOOP THROUGH YEARS OF DATA TO COMPUTE XTB AND VAR(XTB) 


MC8005^ 




DO 200 IT=1, NYR 


MC8005; 




IF(IT .NE. VYR) THEN 


MC8005i 




IF(INV(IN,IT) .NE. 0) THEN 


MC8006( 




X=FTT(INV(IN,IT), Y(IN,IT)) 


MC8006' 




C=SQRT(0. 5+INVC IN, IT) ) 


MC8006; 




XX=X+C*(3. 141592654/2. 0) 


MC8006J 




XT=X/C 


MC8006^ 




T=T+1 


MC8006! 




SUMXT=SUMXT+XT 


MC8006( 




IF(XX .LT. 1.001) XX=1. 001 


MC8006: 




V ARX=AA* ( XX**B 1 ) * ( XX - 1 ) **B 2 


MC8006J 




IF(VARX .GT. 1.0) VARX=1.0 


MC8006< 




VARXT=VARX/(0. 5+INV( IN, IT) ) 


MC8007C 




SUMVAR=SUMVAR+VARXT 


MC8007J 




ENDIF 


MC8007: 




ENDIF 


MC8007: 


200 


CONTINUE 


MC80074 




XTB(IN)=SUMXT/T 


MC80075 




VXTB( IN)=SUMVAR/T**2 


MC80076 


260 


CONTINUE 


MC80077 

MC80078 


— 


CONDUCT ALGORITHM TO FIND XEB 


MC80079 




CALL EBITERC NRC ,XTB , VXTB , XEB , MXX , VYR) 


MC80080 

MC80081 


— 


COMPUTE MEAN SQUARED ERROR 


MC80082 




CALL MSEC INV,Y, NRC, NYR, VYR, XEB, L, MXX, MXY,KL) 


MC80083: 




IF(L .LT. MINL) THEN 


MC80084 




MINL=L 


MC80085^i 




MINLK=KL 


MC80086 




MINLYR=VYR 


MC80087 




ELSE IF(L .GT. MAXL) THEN 


MC80088 




MAXL=L 


MC80089 




MAXLK=KL 


MC80090; 




MAXLYR=VYR 


MC80091! 




ENDIF 


MC80092; 




SUML=SUML+L*KL 


MC80093: 




KLSUM=KLSUM+KL 


MC80094H 

MC80095b 


— 


INVERT XEB TO ORIGINAL SCALE 


MC80096t 




CALL INVERT(NRC,XEB,MXX) 


MC80097f 

MC80098t; 


— 


COMPUTE MAD AND CHI SQUARE 


MC80099C 




CALL OSMOE(INV,Y,NRC,NYR,VYR,XEB,CHI,K,MXX,MXY, 


MC80100( 




FCELLU,FMADU,PMAD,KP) 


MC8010K 




IFCCHI .LT. MINCHI) THEN 


MC80102( 




MINCHI=CHI 


MC80103( 




MNCHIK=K 


MC80104( 




MNCHYR=VYR 


MC80105( 




ELSE IFCCHI .GT. MAXCHI) THEN 


MC80106C 




MAXCHI=CHI 


MC80107C 




MXCHIK=K 


MC80108( 




MXCHYR=VYR 


MC80109( 
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C 

C - — 



19 

21 

25 

26 

27 

28 

29 

30 

31 

32 



ENDIF 

SUMCHI=SUMCHI+CHI*K 
KSUM=KSUM+K 
KPSUM=KPSUM+KP 
SUMMAD=SUMMAD+PMAD*KP 
WRITE( 1 1 , 3 1 ) VYR , KP , FCELLU , FMADU , PMAD 
CONTINUE 

WRITE RESULTS TO OUTPUT FILE 
AVGL=SUML/KLSUM 
AVGCHI=SUMCHI/KSUM 
AVGMAD=SUMMAD/KPSUM 
WRITE(11,19)’AVG MAD = 

WRITEC 11,21) 'CHI 
WRITEd 1,26) 'MIN 
WRITEC 11, 26) 'MAX 
WRITEC 11,26) 'AVG 
WRITEC 11,21) 'MEAN 



’ ,MNCHIK, 'VALID 
’,MXCHIK,' VALID 



YR = 
YR = 



MC801100 
MC801110 
MC801120 
MC801130 
MC801140 
MC801150 
MC801160 
MC801170 
MC801180 
MC801190 
MC801200 
MC801210 
MC801220 
MC801230 
,MNCHYRMC801240 
,MXCHYRMC801250 
MC801260 
MC801270 
MC801280 
MC801290 
MC801300 
MC801310 
MC801320 
MC801330 
MC801340 
MC801350 
MC801360 
MC801370 
MC801380 
MC801390 
MC801400 
MC801410 
MC801420 

C MC801430 

•>V^V**VV**Vf***Vr^VVc">V*^V*Vr********^VVc*****>V**^V**^V***^V********^V***Vt*’>V^V*''>V*^V^V******J^(350 1440 



,AVGMAD 

SQUARE CORIG SCALE): 

CHI = ' ,MINCHI, 'K = 

CHI = ' ,MAXCHI,’K = 

CHI = ' ,AVGCHI 
SQUARED ERROR C TRANS SCALE): 



,MINL, 'K = 
' ,MAXL, 'K = 
' ,AVGL 



,MINLK, 'VALID 
' ,MAXLK,' VALID 



YR = 
YR = 



WRITEC 11,25) 'MIN MSE = 

WRITEC 11,25) 'MAX MSE = 

WRITEC 11,27) 'AVG MSE = 

F0RMATC38X,A,F5. 3) 

F0RMATC/1X,A) 

FORMATC 1X,A,F6. 3 ,5X, A, 13 ,5X, A, 12) 

FORMATC 1X,A,F9. 3,5X,A,I3,5X,A,I2) 

FORMATC 1X,A,F6.3/) 

FORMATC17X,A,2X,A) 

FORMATC 1X,A,4X,A,3X,A,2X,A,3X,A) 

FORMATC 1X,8C'-’),2X,4C'-'),2X,14C’-'),2X,13C'-'),2X,5C’-')) 
FORMATC 1X,I5,I8,8X,F5. 3,10X,F5. 3,6X,F5.3) 

FORMATC IX, A) 

RETURN 

END 



,MINLYR 



,MAXLYR 



C 

C 



SUBROUTINE E BTS 2 C I NV , Y , NRC , NYR , XTB , VXTB , XE B , MXX , MXY ) 
TRANSFORMED SCALE, TIME INDEPENDENT VARIANCE METHOD 
REAL INVCMXX,MXY), YCMXX,MXY), XEBCMXX) 

REAL XTB C MXX ) , VXTB C MXX ) 

REAL MAXL,MINL,L,MAXCHI,MINCHI 

INTEGER T, VYR 

MAXL= -1000.0 

MINL= 1000.0 

SUML= 0.0 

KLSUM=0 

MAXCHI= -1000.0 
MINCHI= 1000.0 
SUMCHI= 0. 0 
KSUM=0 
SUMMAD=0. 0 
KPSUM=0 

WRITEC 11,21) 'EMP BAYES TRANS SCALE - TIME INDEP VAR:’ 
WRITEC 11,21) 'MEAN ABSOLUTE DEVIATION CORIG SCALE):' 

WRITEC 11,28) 'FRACTION CELLS ',' FRACTION MAD’ 

WRITEC 11,29) 'VALID YR' , 'K' , 'WITH UNDERAGE ',’ FROM UNDERAGE 



MC801450 
MC801460 
MC801470 
MC801480 
MC801490 
MC801500 
MC801510 
MC801520 
MC801530 
MC801540 
MC801550 
MC801560 
MC801570 
MC801580 
MC801590 
MC801600 
MC801610 
MC801620 
MC801630 
MC801640 
'MAD' MC801650 
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c 

c 

c 



c 



c 

c 

c 

c 



WRITE( 11,30) 

--- LOOP THROUGH VALIDATION YEARS 
DO 380 VYR=1, NYR 
--- LOOP THROUGH CELLS 
DO 360 1N=1, NRC 

T=0 

SUMXT=0 

SUMXT2=0 

LOOP THROUGH YEARS OF DATA TO COMPUTE XTB AND VAR(XTB) 
DO 300 1T=1, NYR 

1F(1T .NE. VYR) THEN 

1F(1NV(1N,1T) .NE. 0) THEN 
X=FTT(1NV(1N,1T), Y(1N,1T)) 

XT=X/SQRT(0. 5+lNV(lN,lT)) 

T=T+1 

SUMXT=SUMXT+XT 
SUMXT2=SUMXT2+XT**2 
END IF 
ENDIF 
300 CONTINUE 

XTB(1N)=SUMXT/T 

VXTB( 1N)=( (r^SUMXT2) -(SUMXT**2) )/( (T-1)*T**2) 

360 CONTINUE 

--- CONDUCT ALGORITHM TO FIND XEB 

CALL EB1TER(NRC,XTB,VXTB,XEB,MXX,VYR) 

--- COMPUTE MEAN SQUARED ERROR 

CALL MSEC 1NV,Y, NRC, NYR, VYR, XEB, L,MXX,MXY,KL) 

1F(L .LT. MINL) THEN 
M1NL=L 
M1NLK=KL 
M1NLYR=VYR 

ELSE 1F(L .GT. MAXL) THEN 
MAXL=L 
MAXLK=KL 
MAXLYR=VYR 
ENDIF 

SUML=SUML+L*KL 

KLSUM=KLSUM+KL 



--- INVERT XEB TO ORIGINAL SCALE 
CALL 1NVERT(NRC,XEB,MXX) 

--- COMPUTE MAD AND CHI SQUARE 

CALL 0SM0E(1NV,Y,NRC,NYR,VYR,XEB 
* FCELLU,FMADU,PMAD,KP) 

IFCCHl .LT. MlNCHl) THEN 
M1NCH1=CH1 
MNCH1K=K 
MNCHYR=VYR 

ELSE IFCCHl .GT. MAXCHl) THEN 
MAXCH1=CH1 
MXCH1K=K 
MXCHYR=VYR 



MC8016 
MC8016 
MC8016‘ 
MC8016^ 
MC8017: 
MC8017: 
MC8017: 
MC8017; 
MC8017: 
MC8017: 
MC8017: 
MC8017; 
MC8017: 
MC8017 
MC8018 
MC8018 
MC8018 
MC8018 
MC8018 
MC8018. 
MC8018 
MC8018 
MC8018' 
MC8018‘ 
MC8019' 
MC8019 
MC8019' 
MC8019. 
MC8019* 
MC8019; 
MC8019. 
MC8019 
MC8019- 
MC8019' 
MC 8020 ! 
MC8020 
MC8020; 
MC8020; 
MC8020-I 
MC8020; 
MC8020( 
MC80201 
MC8020} 
MC8020? 
MC802K 
MC8021; 

CH1,K,MXX,MXY, MC8021; 

MC8021; 
MC8021- 
MC8021i 
MC802K 
MC8021' 
MC802H 
MC8021< 
MC8022( 
MC8022: 
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c 



380 



19 

21 

25 

26 

27 

28 

29 

30 

31 



ENDIF 

SUMCHI=SUMCHI+CHI*K 
KSUM=KSUM+K 
KPSUM=KPSUM+KP 
SUMMAD=SUMMAD+PMAD*KP 
WRITE( 11,31) VYR , KP , FCELLU , FMADU , PMAD 
CONTINUE 



FILE 



WRITE OUTPUT TO 
AVGL=SUML/KLSUM 
AVGCHI=SUMCHI/KSUM 
AVGMAD=SUMMD/KPSUM 
WRITE(11,19)'AVG MAD = 
WRITE( 11,21) 'CHI 
WRITEd 1,26) 'MIN 
WRITE( 11,26) 'MAX 
WRITE( 11,26) 'AVG 



' ,AVGMAD 
(ORIG SCALE): 



' ,MNCHIK, 'VALID YR = 
' ,MXCHIK,' VALID YR = 



SQUARE 

CHI = ' ,MINCHI, 'K = 

CHI = ' .MAXCHI, 'K = 

CHI = ' ,AVGCHI 
WRITE( 11,21) 'MEAN SQUARED ERROR (TRANS SCALE):' 

WRITE( 11,25) 'MIN MSE = ',MINL,'K = ' ,MINLK, ' VALID YR = 

',MAXL,'K = ' ,MAXLK,' VALID YR = 
' ,AVGL 



WRITEd 1,25) 'MAX MSE = 

WRITEd 1,27) 'AVG MSE = 

F0RMAT(38X,A,F5. 3) 

F0RMAT(/1X,A) 

F0RMAT(1X,A,F6. 3 ,5X, A, 13 ,5X, A, 12) 

FORMATC 1X,A,F9. 3,5X,A,I3,5X,A,I2) 

F0RMAT(1X,A,F6. 3/) 

FORMAT(17X,A,2X,A) 

F0RMAT(1X,A,4X,A,3X,A,2X,A,3X,A) 

F0RMAT(1X,8( ' -' ) ,2X,4( ' - ' ) ,2X, 14( ' - ' ) , 2X, 13( ' - ' ) ,2X,5( ' - ' ) ) 
F0RMAT(1X,I5,I8,8X,F5. 3 , 10X,F5. 3 ,6X,F5. 3) 

RETURN 

END 



MC802220 
MC802230 
MC802240 
MC802250 
MC802260 
MC802270 
MC802280 
MC802290 
MC802300 
MC802310 
MC802320 
MC802330 
MC802340 
MC802350 
,MNCHYRMC802360 
,MXCHYRMC802370 
MC802380 
MC802390 
MC802400 
MC802410 
MC802420 
MC802430 
MC802440 
MC802450 
MC802460 
MC802470 
MC802480 
MC802490 
MC802500 
MC802510 
MC802520 
MC802530 



,MINLYR 

,MAXLYR 



MC802540 

y? y? V? * Vf VrVf Vf Vf Vf VfVf ******Vf Vf ■jVVrVf *** VfVf VfVf '>VV?*Vr'5VVf VfVf Vf 802550 



SUBROUTINE EBOS 1 (INV , Y , NRC , NYR , XTB , VXTB , XEB , MXX , MXY ) 

--- ORIGINAL SCALE, TIME DEPENDENT VARIANCE METHOD 
REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) 

REAL XTB (MXX), VXTB (MXX) 

REAL MAXCHI, MINCHI 
INTEGER T, VYR 
MAXCHI= -1000.0 
MINCHI= 1000. 0 
SUMCHI= 0.0 
KSUM=0 
SUMMAD=0. 0 
KPSUM=0 

WRITE(11,21)'EMP BAYES ORIG SCALE - TIME DEP VAR:' 
WRITE(11,21)'MEAN ABSOLUTE DEVIATION .(ORIG SCALE): ' 

WRITEd 1,28) 'FRACTION CELLS ',' FRACTION MAD’ 

WRITEdl, 29) 'VALID YR' , 'K' , 'WITH UNDERAGE ',' FROM UNDERAGE ',' MAD ' 
WRITE(11,30) 

--- LOOP THROUGH VALIDATION YEARS 
DO 480 VYR=1, NYR 
--- LOOP THROUGH CELLS 



MC802560 

MC802570 

MC802580 

MC802590 

MC802600 

MC802610 

MC802620 

MC802630 

MC802640 

MC802650 

MC802660 

MC802670 

MC802680 

MC802690 

MC802700 

MC802710 

MC802720 

MC802730 

MC802740 

MC802750 

MC802760 

MC802770 
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c --- 



DO 460 IN=1, NRC 

T=0 

SUMXT=0 
SUMVAR=0 

LOOP THROUGH YEARS OF DATA TO COMPUTE XTB AND VAR(XTB) 
DO 400 IT=1, NYR 

IF(IT .NE. VYR) THEN 

IF(INV(IN,IT) .NE. 0) THEN 
PHAT=Y( IN, IT)/INV( IN, IT) 

T=T+1 

SUMXT=SUMXT+PHAT 
IF(PHAT .GT. 0.0) THEN 

SUMVAR=SUMVAR+PHAT*( 1-PHAT)/INV( IN, IT) 

ELSE 

SUMVAR=SUMVAR+l/( INV( IN, IT)+1)**2 
END IF 
ENDIF 
ENDIF 
CONTINUE 
XTB(IN)=SUMXT/T 
VXTB( IN)=SUMVAR/T’'^*2 
460 CONTINUE 



400 



--- CONDUCT ALGORITHM TO FIND XEB 

CALL EBITER(NRC,XTB,VXTB,XEB,MXX,VYR) 



--- COMPUTE MAD AND CHI SQUARE 

CALL OSMOE(INV,Y,NRC,NYR,VYR,XEB,CHI,K,MXX,MXY, 
* FCELLU,FMADU,PMAD,KP) 

IF(CHI .LT. MINCHI) THEN 
MINCHI=CHI 
MNCHIK=K 
MNCHYR=VYR 

ELSE IF(CHI .GT. MAXCHI) THEN 
MAXCHI=CHI 
MXCHIK=K 
MXCHYR=VYR 
ENDIF 

SUMCHI=SUMCHI+CHI*K 
KSUM=KSUM+K 
KPSUM=KPSUM+KP 
SUMMAD=SUMMAD+PMAD*KP 
WRI TE ( 1 1 , 3 1 ) VYR , KP , FCE LLU , FMADU , PMAD 
480 CONTINUE 



WRITE OUTPUT TO FILE 
AVGCHI=SUMCHI/KSUM 
AVGMAD=SUMMAD/KPSUM 
WRITE(11,19)'AVG MAD = ' ,AVGMAD 
WRITE(11,21) 'CHI SQUARE (ORIG SCALE):' 
WRITE(11,26)'MIN CHI = ', MINCHI, 'K = ' 



,MNCHIK,' VALID YR 
WRITE(11,26)'MAX CHI = ', MAXCHI, 'K = ' ,MXCHIK, ' VALID YR 
WRITEC 11,27) 'AVG CHI = ' ,AVGCHI 



MC802' 

MG802i 

MC802J 

MC802( 

MC802q 

MC802t 

MC&02( 

MCS02i 

MC802!i 

MC802f 

MC802( 

MC802{ 

MC802S 

MC802S 

MC802S 

MC802S 

MC8029 

MC802S 

MC8029 

MC8029 

MC8029 

MC8029 

MC803G 

MC803C 

MC8030 

MC8030 

MC8030 

MC8030 

MC8030 

MC8030 

MC8030 

MC8030 

MC8031: 

MC8031: 

MC8031; 

MC8031; 

MC803U 

MC803i:l 

MC8031i| 

MC803l1 

MC8031. 

MC8031'.| 

MC8032(| 

MC8032; 

MC8032:| 

MC8032: 

MC8032^| 

MC8032I 

MC8032(| 

MC8032: 

' ,MNCHYRMC8032f 
' ,MXCHYRMC8032S 
MC8033C 



19 

21 

26 



F0RMAT(38X,A,F5. 3) 

F0RMAT(/1X,A) 

FORMATC 1X,A,F9. 3,5X,A,I3,5X,A,I2) 



MC8033] 

MC80332 

MC80333 
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27 FORMATC 1X,A,F9. 3/) MC803340 

28 FORMAT(17X,A,2X,A) MC803350 

29 F0RMAT(1X,A,4X,A,3X,A,2X,A,3X,A) MC803360 

30 FORMATC 1X,8('-^),2X,4(’-^),2X,14('-'),2X,13('-’),2X,5(’-')) MC803370 

31 F0RMAT(1X,I5,I8,8X,F5. 3,10X,F5. 3,6X,F5. 3) MC803380 

RETURN MC803390 

END MC803400 

C MC803410 



SUBROUTINE EB0S2 ( INV , Y , NRC , NYR , XTB , VXTB , XEB , MXX , MXY) 

C ... ORIGINAL SCALE, TIME INDEPENDENT VARIANCE METHOD 
REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) 

REAL XTB(MXX), VXTB(MXX) 

REAL MAXCHI, MINCHI 
INTEGER T, VYR 
MAXCHI= -1000.0 
MINCHI= 1000.0 
SUMCHI= 0. 0 
KSUM=0 
SUMMAD=0. 0 
KPSUM=0 

WRITEC 11,21) 'EMP BAYES ORIG SCALE - TIME INDEP VAR:' 

WRITEC 11,21) 'MEAN ABSOLUTE DEVIATION (ORIG SCALE):' 

WRITEC 11,28) 'FRACTION CELLS’ , 'FRACTION MAD' 

WRITEC 11,29) 'VALID YR' , ' K’ , ' WITH UNDERAGE ',' FROM UNDERAGE' 
WRITEC 11,30) 

--- LOOP THROUGH VALIDATION YEARS 
DO 580 VYR=1, NYR 
C --- LOOP THROUGH CELLS 
DO 560 IN=1, NRC 

T=0 

SUMXT=0 

SUMVAR=0 

SUMY=0 

SUMINV=0 

C --- LOOP THROUGH YEARS OF DATA TO COMPUTE XTB AND VAR(XTB) 

DO 500 IT=1, NYR 

IFCIT .NE. VYR) THEN 

IF(INV(IN,IT) .NE. 0) THEN 
PHAT=Y( IN , IT) /INVC IN , IT) 

SUMXT=SUMXT+PHAT 

SUMY=SUMY+Y(IN,IT) 

SUMINV=SUMINV+INV( IN, IT) 

T=T+1 

SUMVAR=SUMVAR+1. 0/INV(IN,IT) 

ENDIF 

ENDIF 

500 CONTINUE 

XTB(IN)=SUMXT/T 
IFCSUMY .GT. 0.0) THEN 
PTILDE=SUMY/SUMINV 

VXTBC IN)=(PTILDE*( 1-PTILDE)*SUMVAR)/T**2 
ELSE 

VXTBC IN)=SUMINV*SUMVAR/( ( ( 1+SUMINV)**2)*T’'^*2) 



MC803430 
MC803440 
MC803450 
MC803460 
MC803470 
MC803480 
MC803490 
MC803500 
MC803510 
MC803520 
MC803530 
MC803540 
MC803550 
MC803560 
MC803570 
MC803580 
, 'MAD' MC803590 
MC803600 
MC803610 
MC803620 
MC803630 
MC803640 
MC803650 
MC803660 
MC803670 
MC803680 
MC803690 
MC803700 
MC803710 
MC803720 
MC803730 
MC803740 
MC803750 
MC803760 
MC803770 
MC803780 
MC803790 
MC803800 
MC803810 
MC803820 
MC803830 
MC803840 
MC803850 
MC803860 
MC803870 
MC803880 
MC803890 
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ENDIF 






MC803; 


560 


CONTINUE 






MC803! 

MC803‘ 


— 


CONDUCT ALGORITHM TO FIND XEB 






MC803! 




CALL EBITER(NRC,XTB, VXTB, XEB, MXX, VYR) 






MC803* 

MC803‘ 


— 


COMPUTE MAD AND CHI SQUARE 






MC8035 




CALL OSMOE(INV,Y,NRC,NYR,VYR,XEB,CHI,K,MXX,MXY, 






MC803S 




FCELLU,FMADU,PMAD,KP) 






MC803< 




IF(CHI .LT. MINCHI) THEN 






MC8035 




MINCHI=CHI 






MC804C 




MNCHIK=K ■ 






MC804C 




MNCHYR=VYR 






MC804C 




ELSE IF(CHI .GT. MAXCHI) THEN 






MC804C 




MAXCHI=CHI 






MC804C 




MXCHIK=K 






MC804C 




MXCHYR=VYR 






MC804C 




ENDIF 






MC804C 




SUMCHI=SUMCHI+CHI*K 






MC804C 




KSUM=KSUM+K 






MC804C 




KPSUM=KPSUM+KP 






MC8041 




SUMMAD=SUMMAD+PMAD*KP 






MC804] 




WRITE ( 1 1 , 3 1 ) VYR , KP , FCELLU , FMADU , PMAD 






MC804] 


580 


CONTINUE 






MC8041 

MC8041 




WRITE OUTPUT TO FILE 
AVGCHI=SUMCHI/KSUM 
A VGMAD=SUMMAD / KP SUM 
WRITE(11,19)’AVG MAD = ' ,AVGMAD 
WRITE(11,21)’CHI SQUARE (ORIG SCALE):' 

WRITE( 11,26) 'MIN CHI = ', MINCHI, 'K = ' ,MNCHIK, ' 
WRITE( 11,26) 'MAX CHI = ', MAXCHI, 'K = ',MXCHIK,' 
WRITE(11,27)'AVG CHI = ' ,AVGCHI 






MC8041 

MC8041 

MC8041 

MC8041; 

MC8041 




VALID 


YR = ' 


,MNCHYRMC8042 




VALID 


YR = ' 


,MXCHYRMC8042 








MC8042 


19 


FORMAT(38X,A,F5. 3) 






MC8042 


21 


F0RMAT(/1X,A) 






MC8042 


26 


FORMATC 1X,A,F9. 3,5X,A,I3,5X,A,I2) 






MC8042 


27 


F0RMAT(1X,A,F9. 3/) 






MC8042 


28 


FORMAT(17X,A,2X,A) 






MC8042: 


29 


FORMATC 1X,A,4X,A,3X,A,2X,A,3X,A) 






MC8042; 


30 


FORMATC IX, 8C ' - ' ) , 2X, 4C ' - ' ) , 2X, 14C ' - ' ) , 2X, 13C ’ - ' 


),2X,5C'-')) 


MC8042; 


31 


FORMATC IX, 15, 18, 8X,F5. 3,10X,F5. 3,6X,F5. 3) 






MC80433 




RETURN 






MC80433 




END 






MC8043! 



C MC8043) 



C MC80433 

SUBROUTINE EMTS( INV, Y,NRC ,NYR,XTB , VXTB ,XEB , A,MXX,MXY) MC8043« 

C --- EFRON-MORRIS METHOD MC8043! 

REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) MC8043. 

REAL XTB(MXX), VXTB(MXX), A(MXX) MC80433 

REAL MAXL,MINL,L,MAXCHI,MINCHI MC80441 

INTEGER T, VYR MC8044, 

DATA AA/1. 6835/, B1/-.8934/, B2/. 9881/ MC8044; 

MAXL= -1000.0 MC8044: 

MINL= 1000.0 MC8044^ 

SUML= 0. 0 MC8044; 
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no on 



KLSUM=0 

MAXCHI= -1000.0 
MINCHI= 1000.0 
SUMCHI= 0. 0 
KSUM=0 
SUMMAD=0. 0 
KPSUM=0 

WRITE(11,21)'EFR0N-M0RRIS TRANS SCALE - TIME DEP VAR: ' 
WRITE( 11,21) 'MEAN ABSOLUTE DEVIATION (ORIG SCALE): ' 

WRITE( 11,28) 'FRACTION CELLS ',' FRACTION MAD' 
WRITE(11,29)'VALID YR' , ' K' , ' WITH UNDERAGE ’,' FROM UNDERAGE' 
WRITE(11,30) 

C 

C --- LOOP THROUGH VALIDATION YEARS 
DO 280 VYR=1, NYR 
C --- LOOP THROUGH CELLS 
DO 260 IN=1, NRC 

T=0 

SUMXT=0 

SUMVAR=0 

C ... LOOP THROUGH YEARS OF DATA TO COMPUTE XTB AND VAR(XTB) 

DO 200 IT=1, NYR 

IF(IT .NE. VYR) THEN 

IF(INV(IN,IT) .NE. 0) THEN 
X=FTT(INV(IN,IT), Y(IN,IT)) 

C=SQRT( 0. 5+INV( IN , IT) ) 

XX=X+C*(3. 141592654/2. 0) 

XT=X/C 

T=T+1 

SUMXT=SUMXT+XT 
IF(XX .LT. 1.001) XX=1. 001 
VARX=AA*( XX**B 1 ) *( XX- 1 ) **B2 
IF(VARX .GT. 1.0) VARX=1. 0 
VARXT=VARX/(0. 5+INV( IN, IT) ) 

SUMVAR=SUMVAR+VARXT 

ENDIF 

ENDIF 

200 CONTINUE 

XTB(IN)=SUMXT/T 
VXTB( IN)=SUMVAR/'F'*2 
260 CONTINUE 

--- CONDUCT ALGORITHM TO FIND XEB 

CALL EMITERC NRC , XTB , VXTB , XEB , A , MXX , VYR ) 

--- COMPUTE MEAN SQUARED ERROR 

CALL MSEC INV,Y, NRC, NYR, VYR, XEB, L, MXX, MXY,KL) 

IF(L .LT. MINL) THEN 
MINL=L 
MINLK=KL 
MINLYR=VYR 

ELSE IF(L .GT. MAXL) THEN 
MAXL=L 
MAXLK=KL 
MAXLYR=VYR 
ENDIF 



MC804460 
MC804470 
MC804480 
MC804490 
MC804500 
MC804510 
MC804520 
MC804530 
MC804540 
MC804550 
'MAD' MC804560 
MC804570 
MC804580 
MC804590 
MC804600 
MC804610 
MC804620 
MC804630 
MC804640 
MC804650 
MC804660 
MC804670 
MC804680 
MC804690 
MC804700 
MC804710 
MC804720 
MC804730 
MC804740 
MC804750 
MC804760 
MC804770 
MC804780 
MC804790 
MC804800 
MC804810 
MC804820 
MC804830 
MC804840 
MC804850 
MC804860 
MC804870 
MC804880 
MC804890 
MC804900 
MC804910 
MC804920 
MC804930 
MC804940 
MC804950 
MC804960 
MC804970 
MC804980 
MC804990 
MC805000 
MC805010 
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SUML=SUML+L*KL 

KLSUM=KLSUM+KL 



INVERT XEB TO ORIGINAL SCALE 
CALL INVERT(NRC,XEB,MXX) 



280 

C 

C --- 



19 

21 

25 

26 

27 

28 

29 

30 

31 



COMPUTE MAD AND CHI SQUARE 

CALL OSMOE(INV,Y,NRC,NYR,VYR,XEB,CHI,K,MXX,MXY, 
FCELLU,FMADU,PMAD,KP) 

IF(CHI .LT. MINCHI) THEN 
MINCHI=CHI 
MNCHIK=K 
MNCHYR=VYR 

ELSE IF(CHI .GT. MAXCHI) THEN 
MAXCHI=CHI 
MXCHIK=K 
MXCHYR=VYR 
ENDIF 

SUMCHI=SUMCHI+CHI*K 
KSUM=KSUM+K 
KPSUM=KPSUM+KP 
SUMMAD=SUMMAD+PMAD*KP 
WRITE ( 1 1 , 3 1 ) VYR , KP , FCELLU , FMADU , PMAD 
CONTINUE 



WRITE OUTPUT TO FILE 

AVGL=SUML/KLSUM 

AVGCHI=SUMCHI/KSUM 

AVGMAD=SUMMAD/KPSUM 

WRITE(11,19)'AVG MAD = ' ,AVGMAD 

WRITE( 11,21) 'CHI SQUARE (ORIG SCALE):’ 

WRITE( 11,26) 'MIN CHI = ', MINCHI, 'K = ’ ,MNCHIK, ' VALID YR = 



WRITE( 11,26) 'MAX CHI = ', MAXCHI, 'K = ’ ,MXCHIK, ' VALID YR = 
WRITE(11,26)'AVG CHI = ' ,AVGCHI 
WRITE( 11,21) 'MEAN SQUARED ERROR (TRANS SCALE):’ 

WRITE( 11,25) 'MIN MSE = ',MINL,'K = ' ,MINLK, ' VALID YR = ' ,MINLYR 

',MAXL,'K = ' ,MAXLK, 'VALID YR = ' ,MAXLYR 
' ,AVGL 



MC805C 

MC805C 

MC805C 

MC805C 

MC805C 

MC805C 

MC805C 

MC805C 

MC8051 

MC8051 

MC8051 

MC8051 

MC8051 

MC8051 

MC8051 

MC8051 

MC8051 

MC80511 

MC8052 

MC8052 

MC8052: 

MC8052;| 

MC8052-: 

MC8052!l 

MC8052t| 

MC8052' 

MC8052}| 

MC8052! 

MC8053( 

MC8053: 

MC8053: 

,MNCHYRMC8053: 



,MXCHYRMC8053^ 

MC8053f 



WRITE( 11,25) 'MAX MSE = 

WRITE( 11,27) 'AVG MSE = 

FORMAT(38X,A,F5. 3) 

F0RMAT(/1X,A) 

F0RMAT(1X,A,F6. 3,5X,A, I3,5X,A, 12) 

FORMATC 1X,A,F9. 3,5X,A,I3,5X,A,I2) 

FORMATC 1X,A,F6. 3) 

FORMAT(17X,A,2X,A) 

FORMATC 1X,A,4X,A,3X,A,2X,A,3X,A) 

FORMATC 1X,8('-'),2X,4('-'),2X,14(’-'),2X,13(’-'),2X,5(' 
FORMATC 1X,I5,I8,8X,F5. 3,10X,F5. 3,6X,F5. 3) 

RETURN 

END 



•’)) 



MC8053q 

MC8053) 

MC8053E 

MC80539 

MC8054C| 

MC80541 

MC80542 

MC80543; 

MC80544 

MC80545: 

MC80546 

MC80547 

MC80548: 

MC80549; 

MC80550 

MC805511 






c --- 



SUBROUTINE EB ITER( NRC , XTB , VXTB , XEB , MXX , VYR) 
ITERATIVE ALGORITHM TO SOLVE FOR XEB 
REAL XTB(MXX), VXTB(MXX), XEB(MXX) 

INTEGER VYR 



MC80553 

MC8O554I 

MC80555 

MC80556 

MC80557 
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A=0 

ITER=0 

100 CONTINUE 
ITER=ITER+1 

IF(ITER .GT. 100) PRINT *,'EBITER GT 100 ' 
A0=A 

SUMALK=0 

--- SUM THE ALPHAS 
DO 200 1=1, NRC 

SUMALK=SUMALK+l/( A+VXTB( I ) ) 

200 CONTINUE 

--- COMPUTE XBB 
XBB=0 

DO 300 1=1, NRC 

ALPHA=1/(A+VXTB(I)) 

GAMMA=ALPHA/SUMALK 
XBB=XBB+GAMMA*XTB( I ) 

300 CONTINUE 

--- UPDATE VALUE OF A 
SUMNUM=0 
SUMDEN=0 
DO 400 1=1, NRC 

ALPHA=1/(A+VXTB(I)) 

SUMNUM=SUMNUM+ALPHA*( XTB( I ) -XBB )**2 
SUMDEN=SUMDEN+( (XTB( I) -XBB)**2)*ALPHA**2 
400 CONTINUE 

A=A- ( NRC- 1 -SUMNUM) /SUMDEN 
IF(A .LE. 0) THEN 
A=0 

GO TO 500 
ENDIF 

IF(ABS(A-A0) .GT. 0.0001) GO TO 100 
500 CONTINUE 

--- ITERATIONS CONVERGED, COMPUTE XEB 
DO 600 1=1, NRC 

X=XTB(I) 

V=VXTB(I) 

XE B ( I ) = ( A*X ) / ( A+V ) + ( V*XB B ) / ( A+V ) 

600 CONTINUE 
RETURN 
END 



MC805580 

MC805590 

MC805600 

MC805610 

MC805620 

MC805630 

MC805640 

MC805650 

MC805660 

MC805670 

MC805680 

MC805690 

MC805700 

MC805710 

MC805720 

MC805730 

MC805740 

MC805750 

MC805760 

MC805770 

MC805780 

MC805790 

MC805800 

MC805810 

MC805820 

MC805830 

MC805840 

MC805850 

MC805860 

MC805870 

MC805880 

MC805890 

MC805900 

MC805910 

MC805920 

MC805930 

MC805940 

MC805950 

MC805960 

MC805970 

MC805980 

MC805990 

MC806000 

MC806010 

MC806020 

MC806030 
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C MC806050 

SUBROUTINE EMITER(NRC,XTB,VXTB,XEM,A,MXX,VYR) MC806060 

C --- ITERATIVE ALGORITHM TO SOLVE FOR XEB FOR EFRON-MORRIS METHOD MC806070 

REAL XTB(MXX), VXTB(MXX), XEM(MXX), A(MXX) MC806080 

INTEGER VYR MC806090 

DO 100 1=1, NRC MC806100 

A(I)=0 MC806110 
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100 CONTINUE 




MC806] 

MC806] 


--- SUM THE ALPHAS 




MC806] 


SUMALK=0 




MC8061 


DO 200 1=1, NEC 




MC8061 


SUMALK=SUMALK+1/(A( I)+VXTB( I) ) 




MC8061 


200 CONTINUE 




MC8061 

MC8061 


--- COMPUTE XHAT 




MC8062 


XHAT=0 




MC8062 


DO 300 1=1, NEC 




MC8062 


ALPHA=1/(A(I)+VXTB(D) 




MC8062 


GAMMA=ALPHA/ SUMALK 




MC8062 


XHAT=XHAT+GAMMA*XTB( I ) 




MC8062 


300 CONTINUE 




MC8062 


311 XHATP=XHAT 




MC8062 


1=1 




MC8062 


333 AP=A(I) 




MC8062 


S=(XTB(I)-XHAT)**2 




MC8063 

MC8063 


--- COMPUTE SN AND SD 




MC8063 


SN=0 




MC8063 


SD=0 




MC8063' 


DO 400 J=1,NEC 




MC8063: 


IF(J .NE. I) THEN 




MC8063 


DEN=(A(J)+VXTB(J))**2 




MC8063 


SN=SN+((XTB(J)-XHAT)**2-VXTB(J))/DEN 




MC8063) 


SD=SD+ 1/DEN 




MC80635 


ENDIF 




MC8064i 


400 CONTINUE 




MC8064. 

MC8064: 


--- NEWTON-EAPHSON ITEEATIONS TO SOLVE FOE A 




MC8064; 


444 AD=A(I)+VXTB(I) 




MC8064i 


GNUM=S - 3*VXTB ( I ) +SN*AD**2 




MC8064! 


GDEN=3+SD*AD**2 




MC8064( 


GPEM=2*AD*SN/GDEN-2*GNUM*AD*SD/GDEN**2 




MC8064: 


G=GNUM/GDEN 




MC8064( 


A( I)=A( I) -( (A( I) -G)/( 1-GPEM) ) 




MC8064« 


IF(A(I) .LE. 0.) THEN 




MC8065C 


A(I)=0. 0 




MC8065; 


1=1+1 




MC8065: 


IF(I .LE. NEC) THEN 




MC8065i 


GO TO 333 




MC80654 


ELSE 




MC8065f 


GO TO 555 




MC8065e 


ENDIF 




MC8065: 


ENDIF 




MC8065£i 


IF(ABS(A(I)-AP) .LE. 0.0001) THEN 




MC80659 


1=1+1 




MC8066C 


IF(I .LE. NEC) THEN 




MC80661 


GO TO 333 




MC80662 


ELSE 




MC80663 


GO TO 555 




MC80664 


ENDIF 




MC80665 


ELSE 




MC80666 


AP=A(I) 




MC80667' 

1 
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DO O O 



555 



600 



700 



800 



900 



950 



C 

c 

c --- 



100 



c 



GO TO 444 
END IF 

TEST FOR CONVERGENCE: ABS(S-SP) LT EPSILON 

SUMALK=0 

DO 600 J=1,NRC 

SUMALK=SUMALK+1/(A( J)+VXTB( J) ) 

CONTINUE , 

XHAT=0 

DO 700 J=1,NRC 

ALPHA=1/(A(J)+VXTB(J)) 

GAMMA=ALPHA / SUMALK 
XHAT=XHAT+GAMMA*XTB ( J ) 

CONTINUE 
DO 800 J=1,NRC 
S=(XTB(J)-XHAT)**2 
SP=( XTB ( J ) -XHATP) **2 
IF(ABS(S-SP) .GT. 0.0001) GO TO 311 
CONTINUE 

ITERATIONS CONVERGED, COMPUTE XEM 
DO 950 K=1,NRC 
SD=0 

DO 900 J=1,NRC 

IF(J .NE. K) SD=SD+1/(A(J)+VXTB(J))**2 
CONTINUE 
AD=A(K)+VXTB(K) 

DSTAR=3+(AD**2)*SD 

B=(l, -4. /DSTAR)*VXTB(K)/AD 

IF(B . GT. 1. 0) B=l. 0 

IF(B . LT. 0. 0) B=0. 0 

XEM(K)=XHAT+(1-B)*(XTB(K)-XHAT) 

CONTINUE 

END 



MC806680 

MC806690 

MC806700 

MC806710 

MC806720 

MC806730 

MC806740 

MC806750 

MC806760 

MC806770 

MC806780 

MC806790 

MC806800 

MC806810 

MC806820 

MC806830 

MC806840 

MC806850 

MC806860 

MC806870 

MC806880 

MC806890 

MC806900 

MC806910 

MC806920 

MC806930 

MC806940 

MC806950 

MC806960 

MC806970 

MC806980 

MC806990 

MC807000 

MC807010 

MC807020 






SUBROUTINE MSE( INV,Y,NRC,NYR,VYR,XEB,L,MXX,MXY,KL) 

COMPUTES MEAN SQUARED ERROR MOE 

REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) 

REAL L,MU 
INTEGER VYR 
SUMSE=0 
KL=0 

DO 100 1=1, NRC 

IF(INV(I,VYR) .GT. 0.0) THEN 
X=FTT( INV( I , VYR) , Y( I , VYR) ) 

MU=X/SQRT(0. 5+INV(I,VYR)) 

SUMSE=SUMSE+(XEB( I) -MU)**2 
KL=KL+1 
END IF 
CONTINUE 
L=SUMSE/KL 
RETURN 
END 



MC807040 

MC807050 

MC807060 

MC807070 

MC807080 

MC807090 

MC807100 

MC807110 

MC807120 

MC807130 

MC807140 

MC807150 

MC807160 

MC807170 

MC807180 

MC807190 

MC807200 

MC807210 

MC807220 

MC807230 
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c 

SUBROUTINE OSMOE( INV,Y,NRC,NYR,VYR,XEB,CHI ,K,MXX,MXY, 
* FCELLU,FMADU,PMAD,KP) 

C --- COMPUTES MAD AND CHI SQUARE MOES 

REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) 

INTEGER VYR 
CHI=0. 0 
K=0 

SUMPM0=0. 0 
SUMPMU=0. 0 
KPMO=0 
KPMU=0 

DO 100 1=1, NRC 

P=XEB(I) 

E=P*INV( I ,VYR) 

A=Y(I,VYR) 

C 

C --- COMPUTE MAD FOR THIS CELL 

IF(INV(I,VYR) .GT. 0.0) THEN 
PA=A/INV(I,VYR) 

IF(P .GT. PA) THEN 
SUMPMO=SUMPMO+( P-PA) 

KPMO=KPMO+l 
ELSE 

SUMPMU=SUMPMU+( PA-P) 

KPMU=KPMU+1 
END IF 
ENDIF 

COMPUTE CHI SQUARE FOR THIS CELL 
IF(E .NE. 0.0 .AND. P .NE. 1.0) THEN 
K=K+1 

CHI=CHI+((A-E)**2)/(E*(1-P)) 

ENDIF 

100 CONTINUE 
C 

C --- COMPUTE WEIGHTED AVERAGES 
KP=KPMO+KPMU 

FCELLU=REAL( KPMU) /REAL( KP) 

FMADU=SUMPMU/( SUMPMU+SUMPMO) 

PMAD=( SUMPMU+SUMPMO) /KP 
RETURN 
END 
C 



MC8072 

MC8072 

MC807J 

MC8072 

MC8072 

MC8072 

MC8072 

MC8072 

MC8073 

MC8073 

MC8073 

MC8073 

MC8073 

MC8073 

MC8073 

MC8074 

MC8074 

MC8074 

MC8074 

MC8074 

MC8074- 

MC8074 

MC8074. 

MC8074 

MC8074. 

MC8075 

MC8075' 

MC8075' 

MC8075i 

MC8075.; 

MC8075. 

MC8075' 

MC8075 

MC8075< 

MC8075! 

MC8076t 

MC8076; 

MC8076; 

MC8076J 

MC8076i 

MC8076J 

MC8076( 

MC8076; 

MC8076( 
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C MC8077C 

SUBROUTINE INVERT( NRC ,XEB ,MXX) MC8077] 

C ... INVERT XEB TO ORIGINAL SCALE MC80772 

REAL XEB(MXX) MC80772 

DO 100 1=1, NRC MC8077i 

P=0. 5*(1+SIN(XEB(I))) MC80772 

IF (P .LT. 0.0) THEN MC80776 

P=0. 0 MC80777 

ELSE IF (P .GT. 1.0) THEN MC80778 

P=l. 0 MC8077S 
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ENDIF 
XEB(I)=P 
100 CONTINUE 
RETURN 
END 



MC807800 

MC807810 

MC807820 

MC807830 

MC807840 

MC807850 
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C MC807870 

FUNCTION FTT(INV,Y) MC807880 

C --- CONDUCTS FREMAN-TUKEY TRANSFORM MC807890 

REAL INV,Y MC807900 

TEMP =-l. + 2.*Y/(1.+INV) MC807910 

TEMP1=-1. + 2.*(1.+Y)/(1.+INV) MC807920 

IF(ABS(TEMP). GT. 1 . OR. ABS ( TEMPI ). GT. 1) THEN MC807930 

WRITE(6,*) ’FIT ERROR INV,Y=’ ,INV,Y, TEMP, TEMPI MC807940 

FTT=1 MC807950 

RETURN MC807960 

ENDIF MC807970 

FTT=SQRT(. 5+INV)*. 5*(ASIN(TEMP) + ASIN(TEMPl)) MC807980 

END MC807990 
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D. VECTOR METHOD SUBROUTINE 



SUBROUTINE MC87V(INV,Y,MXX,NYR,NRC,XTBJI, DELTA, X,XVYR,VYRINV, MC800C 

* VYRY , BSTAR , S , GAMMA , XBB J , E VAL , MXP , MXK , BKTBL , NBK , NSC , NCSR , I SFLAG) MC800C 



c 


— 


VECTOR METHOD 


MC800C 






REAL INV(MXX,NYR), Y(MXX,NYR) 


MC800C 






REAL XTBJI(MXP,MXK), DELTA(MXP ,MXK) , X( MXP, MXK) 


MC8000 






REAL XVYR(MXP,MXK), VYRINV(MXP,MXK) , VYRY(MXP,MXK) 


MC8000 






REAL BSTAR(MXP,MXP), S(MXP,MXP), GAMMA(MXP,MXP) 


MC8000 






REAL XBBJ(MXP), EVAL(MXP) 


MC8000 






INTEGER*2 BKTBL(MXX,3) 


MC8000 


c 






MC8001 






REAL MAXL , MINL , L , MAXCHI , MINCHI , MO , MU , MAD 


MC8001 






INTEGER T, VYR, P 


MC8001 


c 






MC8001 






MAXL= -1000.0 


MCSOOl 






MINL= 1000. 0 


MC8001 






SUML= 0.0 


MCSOOl 






KPSUM=0 


MC8001 






MAXCHI= -1000.0 


MC8001 






MINCHI= 1000.0 


MC8001: 






SUMCHI= 0.0 


MC8002. 






KCSUM=0 


MC8002. 






WRITE( 11,32)' ' 


MC8002 






WRITE(11,21)'EMP BAYES TRANS SCALE - VECTOR CASE:' 


MC8002, 






IF (ISFLAG .EQ. 1) THEN 


MC8002.; 






P=NSC 


MC8002: 






WRITE (11, 21) 'VECTOR IS BY SERVICE COMPONENT' 


MC8002: 






ELSE 


MC8002 






P=NCSR 


MC8002i 






WRITEC 11, 21) 'VECTOR IS BY COMMISSIONING SOURCE' 


MC8002' 






ENDIF 


MC8003i 






WRITE(11,22)'K=' ,NRC, 'P=' ,P, 'KP= ' ,(NRC*P) 


MC8003 






WRITEC 11,21) 'MEAN ABSOLUTE DEVIATION (ORIG SCALE):' 


MC8003. 






WRITEC 11,28) 'FRACTION CELLS' , 'FRACTION MAD' 


MC8003; 






WRITEC 11,29) 'VALID YR' , 'KP' , 'WITH UNDERAGE ',' FROM UNDERAGE' 


, 'MAD'MC8003< 






WRITEC 11,30) 


MC80033 






K=NRC 


MC8003t; 






IF(K .LE. (P+2)) THEN 


MC8003‘ 






WRITEC 6,*)'*** error IN VECTOR CASE: P+2 GT K ***' 


MC8003J 






STOP 


MC80033 






ENDIF 


MCS004C 


c 






MC8004:- 


c 


— 


CONDUCT VALIDATION 


MC8004: 






DO 999 VYR=1, NYR 


MC8004' 






DO 90 J=1,MXP 


MC8004^ 






DO 80 1=1, MXK 


MC80043 






XTBJI(J,I)=0. 0 


MC8004( 






DELTAC J,I)=0. 0 


MC8004: 






XVYR(J,I)=0. 0 


MC8004f 




80 


CONTINUE 


MC8004S 




90 


CONTINUE 


MC6005( 






KMKG=BKTBL(1,1) 


MC8005] 






NRC=1 


MC80052 


c 


— 


LOOP THROUGH CELLS IN VECTOR FORM 


Mcsoos:- 
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on on 



DO 130 1=1, NBK 

IF(BKTBL(I,1) .NE. KMKG) NRC=NRC+1 
DO 100 J=1,P 

IF(BKTBL(I,2) . EQ. J) THEN 
JP=J 

GO TO 110 
END IF 

100 CONTINUE 

WRITE(6,*) '*** ERROR IN P VECTOR ASSIGNMENT ***' 

110 T=0 

SUMXT=0 

SUMVAR=0 

C --- LOOP THROUGH YEARS OF DATA TO SOLVE FOR XTB AND VAR(XTB) 

DO 120 IT=1, NYR 

IF(IT .NE. VYR) THEN 

IF(INV(I,IT) .NE. 0) THEN 
XIJ=FTTV( INV( I , IT) , Y(I,IT)) 

C=0. 5+INV(I,IT) 

XT=XIJ/SQRT(C) 

SUMXT=SUMXT+XT 

SUMVAR=SUMVAR+1/C 

T=T+1 

ENDIF 

ENDIF 

120 CONTINUE 

XTB JI ( JP , NRC )=SUMXT/T 

C --- STORE VARIANCE MATRIX IN DELTA MATRIX (TEMPORARY) 

DELTA( JP , NRC ) =SUMVAR/r^*2 

C ... get validation year ESTIMATE, INVENTORY AND ATTRITION INFO 
IF(INV(I,VYR) .GT. 0.0) THEN 
XIJ=FTTV(INV(I,VYR), Y(I,VYR)) 

XT=XIJ/SQRT(0. 5+INV(I,VYR)) 

XVYR(JP,NRC)=XT 

ENDIF 

VYRINVC JP , NRC )=INV( I , VYR) 

VYRY( JP , NRC )=Y( I , VYR) 

KMKG=BKTBL(I,1) 

130 CONTINUE 

IF(K .NE. NRC) THEN 

WRITE (6,*) ’*** ERROR IN VECTOR CASE: K NE NRC ***' 

ENDIF 



--- COMPUTE XBB SUB J 
DO 210 J=1,P 

SUMXTB=0. 0 
DO 200 1=1, K 

SUMXTB=SUMXTB+XTB JI ( J , I ) 

200 CONTINUE 

XBBJ(J)=SUMXTB/K 
210 CONTINUE 

--- COMPUTE X SUB JI MATRIX, MAKE A COPY IN DELTA MATRIX (TEMPORARY) 
DO 230 J=1,P 

DO 220 1=1, K 

X( J , I )=( XTB JI ( J , I ) -XBB J( J) )*SQRT( DELTA( J , I ) ) 

DELTA( J,I)=X(J,I) 



MC800540 

MC800550 

MC800560 

MC800570 

MC800580 

MC800590 

MC800600 

MC800610 

MC800620 

MC800630 

MC800640 

MC800650 

MC800660 

MC800670 

MC800680 

MC800690 

MC800700 

MC800710 

MC800720 

MC800730 

MC800740 

MC800750 

MC800760 

MC800770 

MC800780 

MC800790 

MC800800 

MC800810 

MC800820 

MC800830 

MC800840 

MC800850 

MC800860 

MC800870 

MC800880 

MC800890 

MC800900 

MC800910 

MC800920 

MC800930 

MC800940 

MC800950 

MC800960 

MC800970 

MC800980 

MC800990 

MC801000 

MC801010 

MC801020 

MC801030 

MC801040 

MC801050 

MC801060 

MC801070 

MC801080 

MC801090 
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cj o no on oo ooo oo 



220 CONTINUE 
230 CONTINUE 

--- COMPUTE S MATRIX 

CALL MXYTF(P, K,X,MXP,P,K, DELTA, MXP,P,P,S,MXP) 
---DO EIGENANALYSIS OF S 

--- PUT EIGENVALUES INTO EVAL, EIGENVECTORS INTO GAMMA 
CALL EVCSF( P,S,MXP, EVAL, GAMMA, MXP) 

--- CREATE ESTAR INVERSE 
KP2=K-P-2 
DO 240 J=1,P 

IF(EVAL(J) .LT. KP2) THEN 
EVAL(J)=KP2 
END IF 

EVAL(J)=1. 0/EVAL(J) 

240 CONTINUE 

DO 260 1=1, P 

DO 250 J=1,P 

IF(I ,EQ. J) THEN 
BSTAR(I,J)=EVAL(J) 

ELSE 

BSTAR(I,J)=0. 0 
ENDIF 

250 CONTINUE 
260 CONTINUE 

--- CREATE BSTAR = I - (K-P-2) S TILDE INVERSE 

CALL MRRRRCP, P, GAMMA, MXP, P,P, BSTAR, MXP, P,P,S, MXP) 
CALL MXYTF( P, P,S, MXP, P,P, GAMMA, MXP, P,P, BSTAR, MXP) 
DO 280 1=1, P 

DO 270 J=1,P 

BSTAR( I , J)=KP2*BSTAR( I , J) 

IF(I .EQ. J) THEN 

BSTAR(I,J)=1. 0-BSTAR(I,J) 

ELSE 

BSTAR(I,J)=0. 0-BSTAR(I,J) 

ENDIF 

270 CONTINUE 
280 CONTINUE 

--- COMPUTE DELTA SUB JI 
DO 300 J=1,P 

DO 290 1=1, K 

X(J,I)=XTBJI(J,I)-XBBJ(J) 

290 CONTINUE 
300 CONTINUE 

CALL MRRRRCP, P, BSTAR, MXP, P,K,X, MXP, P,K,XTBJI, MXP) 
DO 320 J=1,P 

DO 310 1=1, K 

DELTA(J,I)=XBBJ(J)+XTBJI(J,I) 

310 CONTINUE 
320 CONTINUE 

--- COMPUTE MSE 



MC8011! 

MC8011 

MC8011 

MC8011 

MC8011 

MC8011 

MC8011: 

MC8011.! 

MC801I 

MC801I 

MC801Z 

MC8012; 

MC8012; 

MC8012. 

MC8012: 

MC8012 

MC8012 

MC8012. 

MC8012 

MC8012*. 

MC8013! 

MC8013! 

MC80133 

MC8013; 

MC8013-; 

MC8013. 

MC8013( 

MC8013' 

MC8013! 

MC8013! 

MC8014( 

MC8014; 

MC8014: 

MC8014; 

MC8014i 

MC8014f 

MC8014( 

MC8014: 

MC8014f 

MC8014S 

MC80150 

MC8015] 

MC80152; 

MC80153 

MC80154 

MC8015S: 

MC80156 

MC80157 

MC80158 

MC80159 

MC80160 

MC80161 

MC80162 

MC80163 

MC80164 

MC80165 
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on no 



KP=0 

SUMSE=0. 0 
DO 340 J=1,P 

DO 330 1=1, K 

IF(INV(I,VYR) .GT. 0.0) THEN 

SUMSE=SUMSE+( DELTA( J , I ) -XVYR( J , I ) )**2 
KP=KP+1 
ENDIF 

330 CONTINUE 
340 CONTINUE 
L=SUMSE/KP 

IF(L .LT. MIND then 
MINL=L 
MINLKP=KP 
MINLYR=VYR 

ELSE IF(L .GT. MAXL) THEN 
MAXL=L 
MAXLKP=KP 
MAXLYR=VYR 
ENDIF 

SUML=SUML+L*KP 

KPSUM=KPSUM+KP 

--- INVERT DELTA SUB JI BACK TO ORIGINAL SCALE 
DO 360 J=1,P 

DO 350 1=1, K 

PHAT=0. 5*( 1+S IN( DELTA( J , I ) ) ) 

IF (PHAT .LT. 0.0) THEN 
PHAT=0. 0 

ELSE IF (PHAT . GT. 1.0) THEN 
PHAT=1. 0 
ENDIF 

DELTA(J,I)=PHAT 
350 CONTINUE 
360 CONTINUE 

--- COMPUTE CHI SQUARE AND MAD 
CHI=0. 0 
KCHI=0 
SUMPMO=0. 0 
SUMPMU=0. 0 
KPM0=0 
KPMU=0 

DO 410 J=1,P 

DO 400 1=1, K 

PHAT=DELTA(J,I) 

E=PHAT’WYRINV(J,I) 

A=VYRY(J,I) 

IF(VYRINV(J,I) .GT. 0.0) THEN 
PACT=A/VYRINV(J,I) 

IF(PHAT .GT. PACT) THEN 
SUMPMO=SUMPMO+( PHAT-PACT) 
KPMO=KPMO+l 
ELSE 

. SUMPMU=SUMPMU+( PACT- PHAT ) 
KPMU=KPMU+1 



MC801660 

MC801670 

MC801680 

MC801690 

MC801700 

MC801710 

MC801720 

MC801730 

MC801740 

MC801750 

MC801760 

MC801770 

MC801780 

MC801790 

MC801800 

MC801810 

MC801820 

MC801830 

MC801840 

MC801850 

MC801860 

MC801870 

MC801880 

MC801890 

MC801900 

MC801910 

MC801920 

MC801930 

MC801940 

MC801950 

MC801960 

MC801970 

MC801980 

MC801990 

MC802000 

MC802010 

MC802020 

MC802030 

MC802040 

MC802050 

MC802060 

MC802070 

MC802080 

MC802090 

MC802100 

MC802110 

MC802120 

MC802130 

MC802140 

MC802150 

MC802160 

MC802170 

MC802180 

MC802190 

MC802200 

MC802210 
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PHAT.NE. 1. 0) THEN 



GT. MAXCHI) THEN 



ENDIF 
ENDIF 

IF(E. NE. 0. 0 .AND 
KCHI=KCHI+1 

CHI=CHI+((A-E)**2)/(E*(1-PHAT)) 

ENDIF 

400 CONTINUE 
410 CONTINUE 

KMAD=KPMO+KPMU 
FCELLU=REAL( KPMU) /REAL( KMAD) 

FMADU=SUMPMU/ ( SUMPMU+SUMPMO) 

PMAD=( SUMPMU+SUMPMO) /KMAD 
IF(CHI .LT. MINCHI) THEN 
MINCHI=CHI 
MNCHIK=KCHI 
MNCHYR=VYR 
ELSE IF(CHI 
MAXCHI=CHI 
MXCHIK=KCHI 
MXCHYR=VYR 
ENDIF 

SUMCHI=SUMCHI+CHI*KCHI 
KCSUM=KCSUM+KCHI 

WRITE( 11,31) VYR,KMAD,FCELLU,FMADU,PMAD 
999 CONTINUE 

AVGL=SUML/KPSUM 
AVGCHI=SUMCHI /KCSUM 

WRITE( 11,21) 'CHI SQUARE (ORIG SCALE):’ 

WRITE( 11,26) 'MIN CHI = ', MINCHI, 'KP = ' ,MNCHIK, 

* 'VALID YR = ' ,MNCHYR 
WRITEC 11,26) 'MAX CHI = ', MAXCHI, 'KP = ' ,MXCHIK, 

* 'VALID YR = ’,MXCHYR 
WRITEC 11,26) 'AVG CHI = ' ,AVGCHI 
WRITEC 11,21) 'MEAN SQUARED ERROR CTRANS SCALE):’ 

WRITEC 11,25) 'MIN MSE = ’,MINL,'KP = ’ ,MINLKP, ' VALID YR = 
WRITEC 11, 25)' MAX MSE = ',MAXL,'KP 
WRITEC 11,27) 'AVG MSE = ' ,AVGL 

21 F0RMATC/1X,A) 

22 F0RMATC1X,3CA,I3,5X)) 

25 F0RMATC1X,A,F6. 3 ,5X,A, 13 ,5X,A, 12) 

26 F0RMATC1X,A,F9. 3,5X,A,I3,5X,A,I2) 

27 F0RMATC1X,A,F6. 3/) 

28 F0RMATC17X,A,2X,A) 

29 F0RMATC1X,A,3X,A,3X,A,2X,A,3X,A) 

30 F0RMATC1X,8C'-'),2X,4C'-'),2X,14C’-’),2X,13C’-'),2X,5C’-')) 

31 F0RMATC1X,I5,I8,8X,F5. 3,10X,F5. 3,6X,F5. 3) 

32 F0RMATC1X,A) 

WRITEC 6,*) 'COMPLETED VECTOR CASE' 

END 



MC8022 
MC8022 
MC8022 
MC8022 
MC8022 
MC8022 
MC8022 
MC8022 
MC8023 
MC8023 
MC8023 
MC8023 
MC8023 
MC8023 
MC8023 
MC8023 
MC8023 
MC8023 
MC8024 
MC8024 
MC8024 
MC8024 
MC8024^ 
MC8024 
MC8024 
MC8024 
MC8024 
MC8024 
MC8025 
MC8025 
MC8025 
MC8025 
MC8025 
MC8025 
,MINLYR MC8025 
’ ,MAXLKP, 'VALID YR = ' ,MAXLYR MC8025 

MC8025} 



_ ! 



MC8025S 
MC8026( 
MC8026 
MC8026 
MC8026 
MC8026^ 
MC8026 
MC8026( 
MC8026 
MC8026? 
MC8026 
MC8027C 

C MC8027 

C MC8027 

MC8027^ 
MC8027 
MC80276 
MC80277 



FUNCTION FTTVCINV,Y) 

CONDUCTS FREMAN-TUKEY TRANSFORM 
REAL INV,Y 

TEMP =-l. + 2.*Y/C1.+INV) 
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TEMP1=-1. + 2.*(1.+Y)/(1,+INV) MC802780 

IF(ABS(TEMP),GT. 1 .OR. ABS(TEMPl). GT. 1) THEN MC802790 

WRITE(6,*) 'FTT ERROR INV,Y=’ ,INV,Y, TEMP, TEMPI MC802800 

FTT=1 MC802810 

RETURN MC802820 

ENDIF MC802830 

FTTV=SQRT(. 5+INV)*. 5*(ASIN(TEMP) + ASIN(TEMPl)) MC802840 

END MC802850 



93 



E. EXEC PROGRAM 



CP LINK MVS 103 103 RR 
ACC 103 K 

*VnVVfVc*V«V********'^r **************** 



FIL CLEAR 

FIL 01 K DSN F0968 MCOR87 DATA (RECFM FB LRECL 69 BLOCK 17940 
DISK MC87 TEMP 

150 

(RECFM F LRECL 25 
(RECFM F LRECL 25 
(RECFM F LRECL 25 
(RECFM F LRECL 25 
(RECFM F LRECL 25 
(RECFM F LRECL 25 



FIL 02 
FIL 06 
FIL 
FIL 
FIL 
FIL 
FIL 
FIL 



25 

27 

29 



DISK 
DISK 

30 DISK 

31 DISK 

32 DISK 
&BEGSTACK 
30. 0 
30 
13 
4 
15 
3 1 
1 16 
1 

&END 

LOAD MC87 



&1 (RECFM FBA LRECL 
DISK MCLASS PG15 
MCLASS PG17 
MCLASS PG19 
MCLASS PG20 
MCLASS PG21 
MCLASS PG22 



2 3 



/* AVG INV THRESHOLD T */ 

/* NO. CELLS THRESHOLD K */ 

/* MOS (ONLY 1) */ 

/* YCS (ONLY 1) */ 

/* GRADE (ONLY 1) */ 

/* NO. SVC COMPS AND ARRAY( 1-REG.2-AUGREG, 3-RES ,4=1+2,5=ALL */ 
/* NO. COMM SRCS AND ARRAY(1-15, 16=ALL) 

/* 3RD DIMENSION (0=N0NE, 1=SVC, 2=CS) 



(START CLEAR 
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F. SAMPLE DATA FILE 



1 


2 


3 


4 


5 


6 


7 


15 


0 


1 


1 


15 


1 


0. 02 


15 


0 


3 


2 


9 


5 


0. 42 


15 


0 


3 


3 


9 


6 


0. 75 


15 


0 


3 


2 


10 


4 


0. 15 


15 


0 


3 


1 


11 


1 


0. 10 


15 


0 


3 


3 


10 


2 


0. 15 


15 


0 


3 


3 


3 


1 


0. 02 


15 


0 


3 


3 


4 


1 


0. 02 


15 


0 


4 


2 


9 


6 


0. 87 


15 


0 


4 


2 


10 


2 


0. 20 


15 


0 


4 


3 


9 


2 


0. 07 


15 


0 


4 


1 


11 


1 


0. 10 


15 


0 


4 


3 


4 


2 


0. 05 


15 


0 


4 


3 


10 


1 


0. 10 


15 


0 


5 


2 


9 


3 


0. 27 


15 


0 


5 


2 


10 


2 


0. 20 


15 


0 


5 


1 


11 


1 


0. 07 


15 


0 


6 


2 


9 


1 


0. 02 


15 


1 


2 


3 


15 


1 


0. 05 


15 


1 


3 


1 


1 


2 


0. 20 


15 


1 


3 


3 


15 


1 


0. 05 


15 


1 


3 


3 


10 


1 


0. 10 


15 


1 


3 


3 


7 


3 


0. 12 


15 


1 


3 


3 


3 


5 


0. 30 


15 


1 


3 


3 


5 


2 


0. 65 


15 


1 


3 


3 


2 


5 


0. 45 


15 


1 


3 


3 


9 


3 


0. 12 


15 


1 


3 


2 


15 


1 


0. 05 


15 


1 


3 


1 


11 


3 


0. 12 


15 


1 


3 


1 


10 


2 


0. 10 


15 


1 


3 


2 


7 


1 


0. 02 


15 


1 


3 


3 


6 


1 


0. 10 


15 


1 


4 


3 


5 


2 


0. 12 


15 


1 


4 


3 


2 


5 


0. 40 


15 


1 


4 


2 


3 


2 


0. 12 


15 


1 


4 


1 


1 


1 


0. 02 


15 


1 


4 


2 


9 


2 


0. 07 


15 


1 


4 


3 


9 


1 


0. 02 


15 


1 


4 


3 


7 


5 


0. 25 


15 


1 


4 


3 


10 


1 


0. 02 


15 


1 


4 


3 


3 


2 


0. 05 


15 


1 


4 


1 


10 


1 


0. 02 


15 


1 


4 


3 


12 


1 


0. 02 


15 


1 


4 


2 


12 


1 


0. 07 



(remaining entries omitted) 



Column descriptions: 

1 - grade 

2 - MOS 

3 - YCS 

4 - service component 



5 - commissioning source 

6 - number of records 

7 - total average inventory 



95 



C --- PROGRAM TO CREATE INVENTORY DATA FILE BY GRADE 
PARAMETER (MXX=20000, MXY=10) 

C --- CLASSIF. TABLE: GRADE, MOS, YCS, SVC, CS 
INTEGER*2 PTRTBL(MXX, 5), NRECS(MXX) 

REAL AINV(MXX) 

INTEGER TYPE , YCS , PG , MOS , SEX , CS , EDLV , SVC , MOS 1 , MOS2 , RACE 

INTEGER DATA(MXY), SPG 

CHARACTER*? CITLS 

DATA AINV/MXX*0./, NRECS/MXX*0/ 

C 

WRITE(5,*) 'ENTER PG' 

READ(5,*) SPG 

WRITE(5,*) 'PG TO USE=' ,SPG 

ICR=0 

NRC=0 

NG=0 

DO 10 1=1,999999 

READ(1,100,END=999) TYPE, YCS, PG, MOS, SEX, CS, EDLV, SVC, MOSl,MOS2, 

* RACE, CITLS, DATA 
ICR=ICR+1 

C --- CLASSIFY ALL RECORDS TYPE 0 
IF(TYPE.GT. 0) GO TO 999 
C --- ADD NEW RECORD TO TABLE 

IF(PG.EQ. SPG) CALL ADDTBL(PG,MOS , YCS ,SVC ,CS , DATA,MXY, PTRTBL, 

* MXX, NRC,AINV,NRECS) 

IF(PG. EQ. SPG) NG=NG+1 

IF(MOD( ICR, 5000). EQ. 0) WRITE(6,*) ' ICR,NRC=' , ICR,NRC 
10 CONTINUE 
C 

999 CONTINUE 

WRITE (6,*) ' ' 

WRITE(6,*) 'TOTAL RECORDS READ =',ICR 

WRITE(6,*) 'TOTAL RECORDS ACCEPTED =' ,NG 

WRITE(6,*) 'TOTAL INVENTORY COMBINATIONS =' ,NRC 

DO 20 1=1, NRC 

WRITE(2,101) (PTRTBL(I,J),J=1,5), NRECS( I) ,AINV( I) 

20 CONTINUE 

100 F0RMAT(3I2,I3,I1,I2,2I1,2I3,I1,A7, IX, 1014) 

101 F0RMAT(I2,I4,I3,I2,I3, 14, F7. 2) 

END 

C 

SUBROUTINE ADDTBL(PG,MOS,YCS,SVC,CS, DATA,MXY, PTRTBL, MXX, NRC, 

* AINV,NRECS) 

C --- SET INVENTORY POINTER FOR THIS ENTRY AND ACCUMULATE 
INTEGER*2 PTRTBL(MXX, 5), NRECS(MXX) 

REAL AINV(MXX) 

INTEGER YCS,PG,MOS,CS,SVC 
INTEGER DATA(MXY) 

MINV=GETINV( PTRTBL, MXX, NRC, PG,MOS,YCS,SVC,CS) 

IFCMINV .EQ. 0) THEN 
C --- NEW COMBINATION 

NRC=NRC+1 

IF (NRC .GT. MXX) THEN 

WRITE(6,*) '*** ERROR - TOO MANY INV. COMBINATIONS' , NRC 
STOP 
ENDIF 



MC8000 

MC8000 

MC8000 

MC8000' 

MC8000 

MC8000 

MC8000; 

MC8000. 

MC8000, 

MC8001; 

MC8001I 

MC80011 

MC80011 

MC8001i 

MC8001; 

MC8001I 

MC8001 

MC8001; 

MC8001'. 

MC8002t 

MC8002; 

MC8002: 

MC8002: 

MC8002< 

MC8002i 

MC8002( 

MC8002; 

MC8002i 

MC8002'' 

MC8003( 

MC8003; 

MC8003: 

MC8003: 

MC8003-: 

MC8003I 

MC8003( 

MC80031 

MC80038 

MC8003S' 

MC8004C 

MC80041 

MC80042 

MC80042: 

MC80044 

MC80045: 

MC80046 

MC80047 

MC80048: 

MC80049' 

MC80050, 

MC80051 

MC80052 

MC80053 

MC 8 0054-1 

MC80055| 

MC80056I 

I 
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MINV=NRC 


MC800570 


PTRTBL(MINV, 1)=PG 


MC800580 


PTRTBLCMINV, 2)=M0S 


MC800590 


PTRTBL(MINV, 3)=YCS 


MC800600 


PTRTBLCMINV, 4)=SVC 


MC800610 


PTRTBLCMINV, 5)=CS 


MC800620 


NRECSCMINV)=0 


MC800630 


ENDIF 


MC800640 


AI=0 


MC800650 


DO 110 IT=1,MXY 


MC800660 


AI=AI + FLOATCDATACIT)) 


MC800670 


CONTINUE 


MC800680 


AINVCMINV)=AINVCMINV) + , 25*AI/MXY 


MC800690 


NRECSCMINV)=NRECSCMINV) + 1 


MC800700 


END 


MC800710 

MC800720 


FUNCTION GETINVCPTRTBL, MXX,NRC, PG,MOS ,YCS , SVC ,CS) 


MC800730 


FIND LOCATION OF MATCHING INVENTORY ENTRY FOR A LOSS 


MC800740 


INTEGER*2 PTRTBLCMXX, 5) 


MC800750 


INTEGER YCS,PG,MOS,CS,SVC 


MC800760 


DO 10 1=1, NRC 

IFCPTRTBLCI, 1) .EQ. PG .AND. 


MC800770 


MC800780 


* PTRTBLCI, 2) .EQ. MOS .AND. 


MC800790 


* PTRTBLCI, 3) .EQ. YCS .AND. 


MC800800 


* PTRTBLCI, 4) .EQ. SVC .AND. 


MC800810 


* PTRTBLCI, 5) .EQ. CS ) THEN 


MC800820 


GETINV=I 


MC800830 


RETURN 


MC800840 


ENDIF 


MC800850 


1 CONTINUE 


MC800860 


GETINV=0 


MC800870 


END 


MC800880 
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APPENDIX C. SAMPLE OUTPUT 



4^ 



A. GENERAL 

This appendix contains sample output from the computer program. A sample out- 
put for test cases one through 30 which use the first five estimation methods is shown 
in paragraph B. A sample output for the vector test cases is shown in paragraph C. 
These examples show the output that is produced by the WRITE statements for file 
definition 11, e.g., WRITE(1 1,101). The program also contains several WRITE and 
PRINT statements that provide interactive information to the user via the terminal 
screen, e.g., WRITE(6,’’'), WRITE(5,*) and PRINT *. This interactive output is 
omitted. 



B. SAMPLE OUTPUT (TEST CASES 1-30) 



TEST CASE INPUT PARAMETERS: 

INVENTORY THRESHOLD= 30. 0 THRESHOLD NO. OF CELLS= 30 

M0S= 13 YCS= 4 GRADE= 15 
SERVICE C0MP0NENTS= 123 

COMM S0URCES= 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

EXPANSION INFORMATION: 

ACTUAL NO. OF CELLS USED= 24 
MOS GROUP # 1 YCS'S USED= 

4 5 

LARGE MOS GROUP #1 YCS'S USED= 

4 5 

MAJOR MOS GROUP #1 YCS'S USED= 

4 5 

EMP BAYES TRANS SCALE - TIME DEP VAR: 

MEAN ABSOLUTE DEVIATION (ORIG SCALE): 

FRACTION CELLS FRACTION MAD 



VALID YR 


K WITH UNDERAGE 


FROM UNDERAGE 


MAD 


1 


24 


0. 458 


0. 478 


0. 127 


2 


24 


0. 250 


0. 187 


0. 099 


3 


24 


0.542 


0. 441 


0. 098 


4 


24 


0. 333 


0. 375 


0. 069 


5 


24 


0. 417 


0. 352 


0. 072 


6 


24 


0. 125 


0. 053 


0. 082 


7 


24 


0. 208 


0. 138 


0. 099 


8 


24 


0. 417 


0. 472 


0. 077 


9 


24 


0. 833 


0. 943 


0. 181 


10 


24 


0. 833 


0. 952 


0. 113 








AVG MAD = 


0. 102 
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24 

24 



CHI SQUARE (ORIG SCALE): 
MIN CHI = 48.590 K 

MAX CHI = 329. 334 K 

AVG CHI = 98.791 



VALID YR = 8 

VALID YR = 9 



MEAN SQUARED ERROR (TRANS SCALE): 



MIN 


MSE = 


0. 033 


K = 


24 


VALID 


YR = 


4 


MAX 

AVG 


MSE = 
MSE = 


0. 205 
0. 079 


K = 


24 


VALID 


YR = 


9 



EMP BAYES 


TRANS 


SCALE - TIME INDEP VAR: 






MEAN ABSOLUTE DEVIATION (ORIG SCALE): 










FRACTION CELLS 


FRACTION 


MAD 




VALID YR 


K 


WITH UNDERAGE 


FROM UNDERAGE 


MAD 


1 


24 


0. 458 


0. 504 




0. 127 


2 


24 


0. 250 


0. 204 




0. 094 


3 


24 


0. 583 


0. 477 




0. 101 


4 


24 


0. 375 


0. 432 




0. 071 


5 


24 


0. 417 


0. 396 




0. 073 


6 


24 


0. 083 


0. 061 




0. 076 


7 


24 


0. 250 


0. 152 




0. 094 


8 


24 


0. 458 


0. 517 




0. 075 


9 


24 


0. 792 


0. 952 




0. 183 


10 


24 


0. 875 


0. 961 




0. 118 








AVG MAD = 


0. 101 


CHI SQUARE 


; (ORIG SCALE): 








MIN CHI = 


45 


.452 K = 24 


VALID 


YR = 


6 


MAX CHI = 


344.445 K = 24 


VALID 


YR = 


9 


AVG CHI = 


99. 284 








MEAN SQUARED ERROR (TRANS SCALE): 






MIN MSE = 


0. 036 K = 24 


VALID YR 


= 4 




MAX MSE = 


0. 213 K = 24 


VALID YR 


= 9 




AVG MSE = 


0. 078 








EMP BAYES 


ORIG 


SCALE - TIME DEP 


VAR: 






MEAN ABSOLUTE DEVIATION (ORIG SCALE): 










FRACTION CELLS 


FRACTION 


MAD 




VALID YR 


K 


WITH UNDERAGE 


FROM UNDERAGE 


MAD 


1 


24 


0. 458 


0. 478 




0. 127 


2 


24 


0. 250 


0. 195 




0. 097 


3 


24 


0. 542 


0. 461 




0. 093 


4 


24 


0. 333 


0. 397 




0. 070 


5 


24 


0. 417 


0. 375 




0. 074 


6 


24 


0. 125 


0. 057 




0. 079 


7 


24 


0. 208 


0. 161 




0. 102 


8 


24 


0. 417 


0. 488 




0. 079 


9 


24 


0. 833 


0. 944 




0. 185 


10 


24 


0. 833 


0. 952 




0. 117 



AVG MAD = 0. 102 
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CHI SQUARE (ORIG SCALE): 

MIN CHI = 49. 207 K = 24 VALID YR = 6 
MAX CHI = 340. 035 K = 24 VALID YR = 9 
AVG CHI = 100. 985 



EMP BAYES ORIG SCALE - TIME INDEP VAR: 

MEAN ABSOLUTE DEVIATION (ORIG SCALE): 

FRACTION CELLS FRACTION MAD 



VALID YR 


K WITH 


UNDERAGE 


FROM UNDERAGE 


MAD 


1 


24 


0. 458 




0. 479 




0. 127 


2 


24 


0. 250. 




0. 196 




0. 097 


3 


24 


0. 542 




0. 462 




0. 093 


4 


24 


0. 333 




0. 399 




0. 070 


5 


24 


0. 417 




0. 377 




0. 075 


6 


24 


0. 125 




0. 056 




0. 079 


7 


24 


0. 208 




0. 161 




0. 101 


8 


24 


0. 417 




0. 490 




0. 079 


9 


24 


0. 833 




0. 945 




0. 185 


10 


24 


0. 833 




0. 953 




0. 117 










AVG MAD = 


0. 102 


CHI SQUARE 


(ORIG SCALE): 










MIN CHI = 


48. 836 


K = 


24 


VALID 


YR = 


6 


MAX CHI = 


339. 835 


K = 


24 


VALID 


II 


9 


AVG CHI = 


100. 960 













EFRON-MORRIS TRANS SCALE - TIME DEP VAR; 

MEAN ABSOLUTE DEVIATION (ORIG SCALE): 

FRACTION CELLS FRACTION MAD 



VALID YR 


K WITH 


UNDERAGE 


FROM UNDERAGE 


MAD 


1 


24 


0. 458 




0. 475 






0. 126 


2 


24 


0. 250 




0. 172 






0. 095 


3 


24 


0. 583 




0. 451 






0. 102 


4 


24 


0. 375 




0. 365 






0. 076 


5 


24 


0. 333 




0. 348 






0. 076 


6 


24 


0. 083 




0. 046 






0. 082 


7 


24 


0. 208 




0. 135 






0. 098 


8 


24 


0. 458 




0. 469 






0. 076 


9 


24 


0. 708 




0. 938 






0. 185 


10 


24 


0. 875 




0. 960 






0. 115 










AVG MAD 


= 


0. 103 


CHI SQUARE 


(ORIG SCALE): 












MIN CHI = 


46. 301 


K = 


24 


VALID 


YR 


= 


8 


MAX CHI = 


340. 712 


K = 


24 


VALID 


YR 


= 


9 


AVG CHI = 


101. 231 














MEAN SQUARED ERROR (TRANS SCALE) 










MIN MSE = 


0. 042 


K = 24 




VALID YR 


= 


4 




MAX MSE = 


0. 211 


K = 24 




VALID YR 


= 


9 




AVG MSE = 


0. 080 















100 



C. SAMPLE OUTPUT (VECTOR TEST CASES) 



TEST CASE INPUT PARAMETERS: 

INVENTORY THRESHOLD= 30. 0 THRESHOLD NO. OF CELLS= 30 

MOS= 151 YCS= 7 GRADE= 17 
SERVICE COMPONENTS= 123 

COMM SOURCES= 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

EXPANSION INFORMATION: 

ACTUAL NO. OF CELLS USED= 8 
MOS GROUP # 8 YCS'S USED= 

7 

LARGE MOS GROUP #3 YCS'S USED= 

7 

MAJOR MOS GROUP #2 YCS'S USED= 

7 



EMP BAYES TRANS SCALE - VECTOR CASE: 

VECTOR IS BY SERVICE COMPONENT 
K= 8 P= 3 KP= 24 

MEAN ABSOLUTE DEVIATION (ORIG SCALE): 

FRACTION CELLS FRACTION MAD 



VALID YR 


KP WITH UNDERAGE FROM UNDERAGE 


MAD 




1 


24 


0. 375 




0. 219 




0. 186 




2 


24 


0. 458 




0. 375 




0. 179 




3 


24 


0. 458 




0.417 




0. 130 




4 


24 


0. 458 




0.422 




0. 127 




5 


24 


0. 542 




0. 708 




0. 146 




6 


24 


0. 292 




0. 310 




0. 142 




7 


24 


0. 250 




0. 193 




0. 126 




8 


24 


0. 375 




0.434 




0. 092 




9 


24 


0. 458 




0. 705 




0. 161 


10 


24 


0. 792 




0. 912 




0. 202 


CHI 


SQUARE 


(ORIG SCALE): 










MIN 


CHI = 


27. 827 


KP 


= 24 


VALID 


YR 


= 8 


MAX 


CHI = 


165. 694 


KP 


= 24 


VALID 


YR 


= 10 


AVG 


CHI = 


61. 025 












MEAN SQUARED ERROR 


(TRANS SCALE): 








MIN 


MSE = 


0. 089 


KP = 


24 


VALID YR 


= 


8 


MAX 


MSE = 


0. 483 


KP = 


24 


VALID YR 




10 


AVG 


MSE = 


0. 229 
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